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Abstract 

Studies of meson production at threshold in the hadron-hadron interaction began in the fifties 
when sufficient energies of accelerated protons were available. A strong interdependence between 
developments in accelerator physics, detector performance and theoretical understanding led to a 
unique vivid field of physics. Early experiments performed with bubble chambers revealed already 
typical ingredients of threshold studies, which were superseded by more complete meson produc- 
tion investigations at the nucleon beam facihties TRIUMF, LAMPF, PSI, LEAR and SATURNE. 
Currently, with the advent of the new cooler rings as lUCF, CELSIUS and COSY the field is 
entering a new domain of precision and the next step of further progress. 

The analysis of this new data in the short range limit permits a more fundamental consideration 
and a quantitative comparison of the production processes for different mesons in the few-body final 
states. The interpretation of the data take advantage of the fact that production reactions close-to- 
threshold are characterized by only a few degrees of freedom between a well defined combination 
of initial and exit channels. Deviations from predictions of phase-space controlled one-meson- 
exchange models are indications of new and exciting physics. Precision data on differential cross 
sections, isospin and spin observables — partly but by no means adequately available — are 
presently turning up on the horizon. There is work for the next years and excitement of the 
physics expected. Here we try to give a brief and at the same time comprehensive overview of this 
field of hadronic threshold production studies. 



PACS numbers: 

13.60.Hb, 13.60.Le, 13.75. Cs, 13.75.-n, 13.85.Lg, 13.85.Ni, 13.85.Rm, 
14.40.Cs, 25.40.-h, 25.40.Ve, 29.20.Dh 



*This review article will be published in "Prog. Part. Nucl. Phys." Vol. 49, issue 1 (2002). 
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8.7 Epilogue 



1 Introduction 



The complexity of the hadronic structure is one of the actual and topical fields of physics concerning 
the microscopic scale. Different, complementary as well as competing experimental approaches are used 
to attack the challenge of the hadronic spectroscopy. 

A detailed understanding of the strong interaction dynamics on the symmetries of QCD is a fun- 
damental aim of hadron physics. At high energies the QCD serves as a reliable theory whereas at low 
energies the non-perturbative QCD phenomena must be investigated. Here calculations are performed 
in an effective field theory of hadrons represented by the chiral effective Lagrangian, a theory of Gold- 
stone bosons coupled to the octet of baryons and vector mesons. The QCD Lagrange function is written 
as: 

LqCD = q {il^d^ — M) q + Lgiuon + Lquark-gluon- (1) 

In principle, q describes the column vector of the quarks -u, rf, s, c, 6, t with each of them represented by 
a spinor, identical to the Dirac equation for spin-| electrons. At energies in the few GeV range only the 
first three quarks are of any relevance, M represents the diagonal of the matrix for the quark current 
masses. 

Data of high quality and precision on hadronic processes at low and intermediate momenta arc 
necessary in order to verify the systematic low energy expansion of the Chiral Perturbation Theory 
(xPT), which has already enforced an important insight into the structure and dynamics of nucleons 
and mesons [1,2]. It is known that the current masses of the lightest three quarks are significantly 
smaller than the typical hadronic energy scale, represented by the proton mass of m^, IGeV/c^. 
This is the reason why in first order the quark masses of the light quarks are neglected in theoretical 
considerations and the spin-| fields can be separated into two independent left- and right-handed 
parts, refiected by the chiral symmetry. The QCD Lagrangian separates into two identical images: 

LqCdIQ-, = Lqcd fee/t; + Lqcd [bright, Qright]- (2) 

Left-handed quarks do not communicate with right-handed quarks and vice versa. There is no 
Lorentz transformation which can change the handedness of a mass-less quark or particle. This sym- 
metry of the Lagrange function is not the symmetry of the spectrum of the particles since otherwise 
each hadron should have a partner of equal mass but inverted parity. Consequently there would be a 
second proton with negative parity which does not exist. There are no parity-doublets in the spectrum 
of the strong interacting particles and therefore the symmetry is spontaneously broken. Due to this 
breaking the mass-less so-called Goldstone bosons, the pseudoscalar tt, K and rj mesons appear. The 
restriction of the interaction of these particles among each other and with other hadrons can be used 
to analyze consequences of the chiral symmetry and its breaking in the framework of an effective field 
theory. 

In reality the lightest hadrons — the Goldstone bosons — are by no means mass-less, which can 
be understood by an explicit breaking of the chiral symmetry including a term for the quark masses. 
For each massive Fermion there exists a Lorentz transformation which transforms a left-handed to a 
right-handed field. The effect of this exphcit breaking now can be treated perturbatively since the 
masses of the light quarks are much smaller than the typical hadronic mass. 

Predicting the low-energy properties of nuclear and particle physics, xPT plays an important role 
in understanding hadron physics in the non-perturbative regime. Therefore it is very important to 
know the properties of the mesons, their structure and interaction in the hadronic environment. 

The physics program at the medium energy hadron accelerators was and is focusing on studies of 
the production and the decay of light mesons and baryon resonances and the conservation or violation 
of symmetries. 
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Following this sequence we will present and discuss production of mesons and meson pairs at thresh- 
old with the questions depicted by the interaction view of figure 1 [3] where in the nucleon-nucleon 
{NN) scattering a meson X is created in a one-boson-exchange model. For the particular case the 
questions have to be answered: how is the distortion of the incident NN waves (ISI) included, which 
mesons contribute to the exchange process, is there an intermediate baryon resonance, how significant 
are rescattering contributions of the exchange mesons and what is a reasonable treatment of the A^A^ 
and NX final state interactions (FSl)? Especially for the A^A^-interaction the FSl is crucial because 
of the nearby poles in the S-wave amplitudes corresponding to the deuteron bound state in the ^Si 
channel or the ^Sq virtual state [3]. These poles and the phase-space factors tend to determine much 
of the energy dependence of the total cross section for meson production. Furthermore, in any region 
where these poles dominate, it is possible to link quantitatively meson production in cases where the 
two nucleons emerge separately or as a bound deuteron state. 



NN^NNX 




Figure 1: Diagrammatic view of the one-boson-exchange meson production process. Produced meson 
= X, nucleon momentum = p, meson momentum = q, relative A^A^ momentum = 2 k [3]. 

Even at threshold the reaction mechanism of the basic process for the interrelation between real 
pions and virtual exchange meson currents as for instance in the pion production (A^A^ NNtt) is still 
not fully understood. Though first data suggested that s-wave pions were produced in a heavy-meson- 
exchange process, later measurements resulted in an interpretation of an interference among transition 
amplitudes as Ss, Ps, Pp, Sd and Ds, where the capital letter indicates the A^A^ final state wave and the 
small letter the angular momentum between the two nucleons and the meson produced. In addition, it 
has been concluded that for higher partial waves vr exchange rather than a heavier meson exchange is 
more significant. 

Especially the pion-nucleon interaction went through several stages of increasing sensitivity for tests 
of nuclear theory. As indicated above more than a decade after the initial work [4, 5] only in the sixties 
explicit calculations on cross sections of the pp ppn^, pmr^ transitions were possible [6, 7] which, 
however, failed the total cross section [8-10] by a factor of five. After further theoretical developments, 
see e.g. reference [11], only recently detailed calculations were provided within the Jiilich meson exchange 
model [12-14] including transitions in the threshold region beyond = and predicting analyzing 
powers and spin correlation coefficients making use of i) the basic diagrams, ii) final state interactions, 
iii) off-shell effects, iv) the exchange of heavier mesons and v) influences of intermediate excitation of 
the A resonance. On the other hand, the Jiilich model does not account for quark degrees of freedom 
i.e. calculations within the framework of the xPT which, in fact, predict cross sections too small [15] 
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compared to experiments. 

The heavier the meson produced in nucleon-nucleon scattering the larger the momentum transfers and 
thus in such processes the short-range parts of the production operators are tested. 

Double meson production processes with both mesons being either identical (e.g. 7r°7r° or K^K'^) or 
different (e.g. nr]) are in principle similar to the single ones, however, the possibly associated baryon 
resonances as intermediate states differ significantly. For instance, the two-pion production very likely 
is dominated by either the Pn A^*(1440) resonance — via an effective a exchange giving information of 
a scalar meson excitation of the Pn A^*(1440) resonance — or the simultaneously excited P33 A(1232) 
resonances. Due to selection rules here the choice of definite quantum numbers in the initial and final 
state might help to sort out the reaction mechanism. 

In any case, a large amount of the knowledge about the tt-tt interaction has been obtained by the 
ttN tittN reaction, whereas data employing proton beams just start to come. 

As long as only S-wave processes are involved, which is the privilege of threshold production studies, 
scattering length and effective range approaches are used to describe the interaction sign and strength. 
In case of strong attractive interactions a distinction between the final state scattering and the forma- 
tion of a baryonic resonance leading to bound or quasi-bound states can not be made uniquely. Such 
investigation will be presented. 

Isospin violation or charge symmetry breaking processes are a topical and interesting field of threshold 
production physics. The essential contribution to isospin breaking is the possible ir^ — rj mixing. There- 
fore experiments around the rj threshold as e.g. dd — > a7r'^{r]) or pd — > ^He7r^{^H n'^) should identify 
the TT^ — r] mixing angle which in turn might give hints to the mass difference of the up and down quarks. 

1.1 Aspects of threshold production 

Threshold production experiments are characterized by excess energies which are small compared to the 
produced masses ^. Consequently, in fixed target experiments the momenta of the final state particles 
transverse with respect to the direction of the incoming beam are small compared to the longitudinal 
components. Thus, ejectiles are confined to a narrow forward cone in the laboratory system around the 
centre-of-mass movement and — close-to-threshold — an experimental acceptance covering the full 
solid angle is feasible with comparatively small dedicated detector arrangements. 

Small relative momenta in the final state effectively limit the number of partial waves contribut- 
ing, simplifying the theoretical interpretation of experimental results. It should be noted, that already 
three-body final states require — in principle — a three-body Faddeev like approach, which has not 
been accomplished so far [16, 17]. However, as first described by Watson [18], with two strongly inter- 
acting particles in the final state, the energy dependence of the total cross section close-to-threshold 
is essentially determined by the (three-body) phase-space and the energy dependence of the on-shell 
final state interaction (FSI). Due to small relative velocities FSI effects are inherent to the experimental 
observables. Thus, the interpretation of data in terms of reaction dynamics requires a correct treatment 
of both initial ^ (ISI) and final state interactions [21] (for a review see [22]). On the other hand, FSI 
effects might provide access to low-energy scattering parameters, which are difficult to obtain otherwise 
in case of unstable particles. 

^With meson production up to masses of w 1 GeV/c^ discussed in this article, we will principally refer to experimental 
results and theoretical investigations covering excess energies up to 140 MeV (the 7r-mass) in excess of the respective 
thresholds. 

^At high relative momenta of the colliding nucleons necessary for meson production the on-shell nucleon-nucleon 
interaction exhibits a rather weak energy dependence and might well be approximated by an overall reduction of the cross 
section in magnitude [19,20]. However, for produced masses in the IGeV/c^ regime, even for proton-proton scattering 
no reliable model exists to allow a consistent evaluation of the initial state interaction [20]. 
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Meson production at threshold imphes high momentum transfers Ap, given by 



Ap = -y (mjv + niB + mx)^ -4m2^ (3) 

for the reaction type NN — > NBX with iV, B and X denoting a nucleon, baryon or meson in the 
initial and final state and m^v, niB and m^ as the respective masses. With momentum transfers in 
the range of 0.37GeV/c to l.lOGeV/c for 7r° and meson production, threshold production probes 
with corresponding distances between 0.53 fm and 0.18 fm the short range part of the nucleon-nucleon 
interaction. Consequently, the energy dependence of the primary production amplitude is expected to 
be weak, motivating Watson's approach [18]. 

Theoretical studies have been carried out mainly within the framework of hadronic meson exchange 
models, i.e. with baryons and mesons as effective degrees of freedom (for a recent review see [20]). Chiral 
perturbation theory has been applied for the description of data on tt production close-to-threshold [23]. 
However, in view of the characteristic distances mentioned above, QCD inspired models [24, 25] with 
constituent quarks and gluons or instantons as relevant degrees of freedom, might turn out to be 
appropriate. High quality exclusive data on close-to-threshold meson production in the energy range of 
non-perturbative strong interaction physics will be crucial for exploring the boundary between effective 
meson exchange models and (so far phenomenological) approaches based on quark-gluon degrees of 
freedom. 



1.2 Hadronic and electromagnetic probes 

It seems that the X-P^ is of limited value once strangeness is involved in the hadronic systems and 
non-perturbative coupled-channel considerations beyond xPT are on the market [26,27]. High quality 
data on the meson {7r,f],f]',LV,K~^,K~) production at threshold have been and arc being produced at 
the hadronic beam cooler rings using the baryon number B = 2 systems (essentially pp and pn). At 
least equally important are studies with hadronic interactions in the B = 1 sector. However, most 
measurements with tt- or X-beams used bubble chamber techniques which naturally suffer from very 
poor statistics. 

Accurate investigations of differential cross sections for meson production from the initial tt plus 
proton or K plus proton systems arc needed at several momenta to study the s- and p-wave dynamics 
and their interplay. There are several indications that as intermediate states both the Sn resonance 
and the quasi-bound state — followed by a final state interaction leading to nucleon and meson 
— are essential and could be investigated best via the B = 1 sector. 

Until now, most of the data for threshold meson production were observed with proton beams. The 

use of antiprotons has contributed considerably to the knowledge of the meson spectroscopy up to serious 
candidates for glueballs and exotic hybrids but is still suffering from the lack of unique interpretations. 

A comparison of the meson production in the fundamental process to productions on heavier nucleons 
would give constraints to eliminate uncertainties in the basic interactions. As outlined by W. Weise 
and U.-G. MeiBner [28], precise measurements of such processes set the necessary constraints for the 
effective Lagrangian of low-energy QCD with strange quarks. Its detailed knowledge has impacts on 
several other important issues, such as the strange quark content of the nucleon. 
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Figure 2: Proton-proton cross sections for different reactions of the light meson production in the 
threshold region. For comparison also the total and elastic cross sections of proton-proton scattering 
are shown. Relevant references to the data will be given throughout the present article at the appropriate 
sections. 



Though in the present article we will restrict ourselves to hadron induced reactions, in general 
induced with p-, c?-, vr- and p-beams, investigations on the structure of the nucleon as seen with 
leptons should at least be mentioned [29]. The use of leptonic interactions via photon or electroweak 
boson exchanges has the attraction of combining a weakly interacting particle of well controlled energy 
momentum transfer with the complex multi particle system containing charged constituents. Photo- 
production provides a crucial complementary information compared to hadronic processes, since they 
are characterized by the electromagnetic probe as a virtual vector meson or an ss fluctuation and as 
such are a powerful source of spin-1 states. The well proven perturbative methods allow us to define a 
hierarchy of operators which, at least in principle, can be extracted separately. Perturbative methods 
finally allow to overcome one of the limitations of the electroweak interaction, that is the access to gluon 
distributions. In fact the strong interaction in the limit of hard scales is becoming a probe for the soft 
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region. All this can be complemented by the addition of the spin and flavour degrees of freedom [29]. 

Cross sections are very small in the threshold region down to eight orders of magnitude compared to 
the total yield. Due to the rapid growth of the phase-space volume the total cross section of the meson 
production reactions increases by orders of magnitudes in a few MeV range of excess energy. Such 
studies have been made possible only due to the low emittance and small momentum spread proton 
beams available at the storage ring facilities. Figure 2 demonstrates threshold production measurements 
of light mesons in the proton-proton interaction. 

Especially hadron accelerators offer the possibility to excite nuleonic resonances as A^*'s and A^*^'s 
with the same quantum numbers as those of the proton via the pa N*a reaction, where the a-channel 

— with its spin and isospin zero — serves as a spin-isospin filter. 

This scalar-isoscalar tool is complementary to the vector excitation of nuclei with electromagnetic 
probes. Recent results from different electron beam facihties as ELSA and MAMI in Germany and 
NIKHEF in The Netherlands certainly deserve a detailed and circumstantial presentation. Due to the 
present limitation to hadronic interactions we can not give appropriate credit to this field of physics and 
refer to [29,30]. At the Continuous Electron Beam Accelerator (CEBAF) of the Thomas Jefferson Na- 
tional Accelerator Facility (TJNAF) an electron beam is used for exciting particular nuclear resonances 
and observing their individual decay modes. 

1.3 Hadronic interactions 

The interaction of hadrons — being the reflection of the strong force acting between their constituents 

— delivers indirect information about their structure and the strong interaction itself. In the frame of 
the optical potential model the hadronic interaction can be expressed in terms of phase-shifts, which in 
turn are described by the scattering length and effective range parameters. These are quite well known 
for the low-energy nucleon-nucleon interaction [31, 32] yet they are still poorly established in case of 
meson-nucleon or even meson-meson interactions. This is partly due to the absorption of mesons when 
scattering on a baryon. To account for this effect the scattering length becomes a complex quantity 
where the imaginary part — e.g. in case of the nucleon-r/ interaction — describes the rjN — > ttN 
and r]N — > multi-n N processes. Moreover, the short life-time of all neutral ground state mesons 
prohibits their utilization as secondary beams and therefore the study of their interaction with hadrons 
is accessible only via observations of their influence on the cross section of the reactions in which they 
were produced (A'"A'" — > NN Meson). When created close to the kinematical threshold with the relative 
kinetic energy being in the order of a few MeV, the final state particles remain much longer in the range 
of the strong interaction than the typical life-time of A^* or A baryon resonances with 10^^'^ s and hence 
they can easily experience a mutual interaction before escaping the area of an infiuence of the hadronic 
force. This — as introduced by Watson [18] — final state interaction (FSI) may significantly modify 
both the original distributions of relative momenta of the outgoing reaction products and the magnitude 
of the production cross section. Though it is easy to understand that the FSI changes the distributions 
of the differential cross sections it is rather difficult to cope for the influences on the magnitude of the 
total reaction rate, since one tends to separate the primary production from the final state interaction 
in space as well as in time [18]. Considering the primary production as a separate process it is well 
worth trying to understand the phenomenon qualitatively. If there were no final state interactions the 
total cross section would be fully determined by the kinematically available phase-space volume, Vps, 
where each interval is populated with a probability governed by the primary production amplitude only: 



The approximation in the equation results from the assumption that |Mp = constant in a few MeV 
range above the production threshold [33-35]. F denotes the flux factor of the colliding particles. 




const. ■ K 



(4) 
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In the classical picture we might imagine that the reaction particles are created together with 
their appropriate force field and when escaping the interaction region they acquire a potential energy 
increasing or decreasing their kinetic energy depending whether the interaction is repulsive or attractive. 
For an attractive interaction they could be created also in those phase-space partitions which are not 
available for non-interacting particles and subsequently be "pulled down" to the energetically allowed 
regions by final state interaction. The temporary growth of the primary production phase-space would 
than increase the reaction rate. Contrary, in case of a repulsive interaction the particles must be 
produced in the lower phase-space volume since leaving the interaction area they will acquire additional 
kinetic energy. The reduction of the total cross section, for example in case of the repulsive Coulomb 
force, can easily be understood when considering the production in a coordinate space of point-like 
objects. Here — in contrast to non-interactive particles — a strongly repulsive object can not be 
produced at appropriately small distances since their later acceleration would lead to the breaking of the 
energy conservation and thus the space available to primary production is reduced. One can also argue 
relativistically that the primary mechanism creates the particles off the mass shell and subsequently 
they are lifted onto the mass shell by the final state interaction. The solid line in figure 3 (left) depicts 
the boundary of the Dalitz plot in case of the pp — > pprj reaction calculated at the total centre-of-mass 
energy a/s = 2433.8 MeV exceeding by 10 MeV the threshold energy. The area surrounded by that 
curve is a direct measure of the kinematically available phase-space volume. The dotted line shows the 




2.19 



2.19 2.2 2.21 2.22 2.23 2.24 2.25 

(mp^)' [ GeV V ] 



2.26 
^ 2.25 

^ 2.24 
CD 

o 

^ 2.23 

CM 

!^ 

E 2.22 
2.21 
2.2 



■ □ □ n □ □□□□ 

□ □□□□□ DDl^ 

□ □□□□□□□□I 
□□□□□□□□□ 

□□□□□□ 



□ □ 

□ □ 



□ □ □ □ 



□ □ 

□ □ 

□ □ 



□ □ n □ □ 



□ □□□□□□ 



□ □ 

□ □ 

□ □ 

□ □ 

□ □ □ □ □ 

□ □ □ □ □ 

□ □ □ □ □ 



□ 
□ 



□ □ □ 

□ □ □ □□ 

□ □ □ □ □ 
□ 
□ 
□ 
□ 
□ 



□ □ 

□ □ 

□ □ 

□ □ 

□ ■ 



□ □□□□□□ 



2.2 2.21 2.22 2.23 2.24 2.25 2.26 

(nip^)^ [ GeVV ] 



Figure 3: (left) The solid line indicates the kinematically available area for the pprj system with a total 
centre-of-mass energy y/s = 2433.8 MeV. The dotted fine depicts the range assuming a reduction of 
the proton and ?7-meson masses by 2 MeV. The phase-space volume results from an integral over the 
closed area: Vps = / / dmp_^^dmp^^. (The powers of 2n are skipped here and will be included into 
the fiux factor F according to the convention introduced by Byckling and Kajantie [36]). 
(right) Distribution of the phase-space for the pprj system modified by the proton-proton interaction 
and calculated for an excess energy of Q = 16 MeV. The area of squares are proportional to the number 
of entries. The largest square of the figure corresponds to 260 events. 



corresponding plot at the moment of the primary creation if the mass of each particle was reduced 
by 2 MeV, demonstrating that now the available phase-space grows significantly. Indeed, as shall be 
inferred from the experimental results presented in subsequent sections, at excess energies of a few MeV 
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above threshold, the mutual interaction among the outgoing particles enhances drastically — by more 
than an order of magnitude — the total cross section and modifies appreciably the occupation of the 
phase-space. Figure 3 (right) indicates the phase-space distribution expected for the pprj system at 
an excess energy of Q = ^/s — 2 m^ — m^ = 16 MeV, assuming a homogeneous primary production 
and taking into account the S-wave interaction between the protons. The proton-proton FSI modifies 
the homogeneous Dalitz plot distribution of "non-interacting particles" , enhancing its population at a 
region where the protons have small relative momenta. The interaction of the proton-77 system would 
manifest itself at low invariant masses m^^ corresponding to small relative momenta between the proton 
and the 77 meson. Such effects observed in the experiments are presented in section 4. 

At the vicinity of the threshold, relative angular momenta larger than I = Oh will not play any role 
due to the short range of the strong interaction and small relative momenta of the produced particles. 
It can be inferred from parity and angular momentum conservation laws that the production of an 
N-N-Meson system (for pseudoscalar or vector mesons) in relative S-waves may only occur if the 
nucleons collide in P-wave. Thus at threshold the transition P ^ Ss is the only possible one, with 
capital letters denoting the angular momentum between nucleons and the small letter for the meson 
angular momentum with respect to the pair of nucleons. An interaction of the nucleons in the entrance 
channel infiuences the production process appreciably [21], similarly as that described above between the 
outgoing particles. For example, for 77 production in the pp pprj reaction the initial state interaction 
(ISI) reduces the total cross section by about a factor of 3-5 [19,21] due to the repulsive proton-proton 
^Po-wave potential. This factor keeps constant in the range of a few tens of MeV [19] and hence, 
does not infiuence the energy dependence of the total cross section of the meson production, which 
remains predominantly determined by the final state interaction. In particular, the dominant S-wave 
nucleon-nucleon final state interaction is by far stronger than any of the low-energy meson-nucleon 
ones, with the only exception of the proton-i^'~ FSI (see section 3). The effect of the nucleon-nucleon 
FSI diminishes with increasing excess energy since it significantly infiuences only that partition of the 
phase-space at which the nucleons have small relative momenta. Whereas this fraction keeps constant, 
the full volume of the phase-space grows rapidly: An increase of the excess energy from Q = 0.5 MeV 
to Q = 30 MeV corresponds to a growth of Vpg by more than three orders of magnitude [37]. As a result 
the S-wave nucleon-nucleon FSI is of less importance for higher excess energies where it affects a small 
fraction of the available phase-space volume only. A more quantitative discussion about the influence 
of the ISI and FSI effects at close-to-threshold production cross sections will be presented in section 3. 



2 Basics of free and quasi— free production 
2.1 Free NN scattering 

Investigations of the production of mesons and their interactions with nucleons are based on measure- 
ments determining the total and differential production cross sections and their dependence on the 
energy of the interacting nucleons. Therefore, to enable a quantitative discussion on the mechanisms 
leading to the transformation of the energy-of-motion of nucleons into matter in the form of mesons 
let us recall the formula of the reaction cross section. In case of the NN — > A^A^ Meson process — with 
the four-momenta of the colliding nucleons denoted by and and with n = 3 particles in the exit 
channel — this reads: 

^ = 4 / dV,,\M,t,-.i23\' i{d% ■ 8{r, = m^) . e(P^) • 5^(P„ + P6 - X^P,-) • |M„,_.i23p, (5) 

where |Mab-»i23p denotes the square of the Lorentz-invariant spin averaged matrix element describing 
the probability to create two nucleons and a meson with four-momenta of Pj = (Ej,Pj) and i — l..n, 
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respectively. The energy and momentum conservation as well as the on-shellness of the created particles 
is ensured by the Dirac-(5 and the Heaviside-O functions. The formula holds also for n 7^ 3. 

The total cross section is then defined as an integral of the probabilities to populate a given phase- 
space interval over the whole kinematically available range — determined by energy and momentum 
conservation — normalized to the flux factor F of the colliding nucleons. In the case of non-interacting 
final state particles the matrix element close to threshold |Ma5^i23p, is nearly constant [33-35] and 
hence the allowed phase-space volume, Vps, is the decisive quantity which governs the growth of the 
cross section with increasing excess energy Q. The latter — defined as the total kinetic energy — is 
shared among the outgoing particles in the reaction centre-of-mass frame ^: 

n 

Vs-^m, (6) 

i=l 

where s = |Pa + Pfep = | denotes the square of total centre-of-mass energy. Exactly at 

threshold, where the particles' kinetic energy in centrc-of-mass system is equal to zero (Q = OMeV), 
the total reaction energy ^/s amounts to: ^/sth = ^"=1 nij. Before writing explicitly the formula for Vps 
let us introduce the kinematical triangle function A defined as [36]: 

A(x, y, z) = + + 2;^ - 2xy — 2yz — 2zx, (7) 

which enables the expression of many useful kinematical variables in a very compact and Lorentz- 
invariant form. In particular, the momenta of particle i and j in their rest system equal: 



p: = p; = ^ o ^ — > (8) 



/Sij 



where Sij = |Pj + P^p stands for the square of the invariant mass of the ij system considered as one 
quasi-particle. The above relation gives an expression for the flux factor F from equation (5) in terms 
of the colliding masses of the nucleons and the total energy s only, namely: 



F = 4V~s {2nf^-' p: = 2 {2n)" ^X{s,ml,ml), (9) 

where we have chosen the convention introduced by Byckling and Kajantie [36] and included the (2n) 
factors for the phase-space (27r)^" and for the matrix element (27r)~^ into the definition of F. It is 
important to note that at threshold, for an excess energy range of a few tens of MeV, the small 
fractional changes of total energy = ; — -, -7^) causes weak variations of the flux factor and 

o \ mi + m2 + ma + Q / 

influences only slightly the shape of the energy dependence of the total cross section. 

For a two-particle flnal state ab — > 12 (for instance for the reaction pn drj) the phase-space 
integral defined in equation (5) reduces to Vps := / dVps = Pi = f^A/A(s, mf,m2) and the total cross 
section for such reactions (when neglecting variations due to the dynamical effects (|Mab^i2| = const.)) 
should increase linearly with the centre-of-mass momentum of the produced meson in the vicinity of 
the threshold. The total cross section can be expressed analytically as a function of the masses of the 
particles participating in the reaction and the square of total reaction energy s: 



Vps const 1 v/A(s.iii2.iii^) 

(^ab^i2 = const ■ = — = — — =. 10) 

F 167rsp* 167rs^A(s, m2,m^) 



^Traditionally by centre-of-mass system we understand a frame in which the momenta of all particles add to zero, 
called sometimes more explicitly centre-of-momentum frame. 
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Near threshold, at a given excess energy Q = y^ — mi — m2, the emission of the reaction products in the 
centre-of-mass frame is isotropic and the whole dynamics of the process manifests itself in the absolute 
value of the transition matrix element |Mafe^i2|. The underlying production mechanisms can also be 
extracted from the deviations of the total cross section energy dependence following the prediction 
of relation (10). A visualization of possible differences in the dynamics for the production of various 
mesons can be realized by comparing the total cross section of the studied reactions at the same value 
of the phase-space volume normalized to the flux factor (Vps/F). It is important to recognize that 
the comparison at the same centre-of-mass meson momentum (pg) gives direct information about the 
amplitude differences only if the production of mesons with the same masses is concerned (for example 
pp —>■ dTT"*" and nn — >• dTT~), yet the reactions pn — >• drj and pn — >• (i7r° have by far different Vps at the 
corresponding P2. 

However, if one is interested in the decomposition of the production amplitude according to the 
angular momenta of the final state particles then indeed the appropriate variable for the comparison 
even for particles with different masses is the meson momentum in the reaction centre-of-mass frame. 
Correspondingly for the more than two-body final state, the adequate variable is the maximum meson 
momentum, since it is directly connected to the maximum angular momentum by the interaction range. 
Considering for a three-body exit channel {ab — > 123) that the meson (m^) possesses maximum mo- 
mentum when the remaining two particles are at rest relative to each other and employing definition (8) 
one obtains: 



_ VA(s,(ml+m2)^mi) 
^max 2Vs ■ ^ ' 

Contrary to the two-body final state, at a fixed excess energy the dynamics of the meson production 
associated with two or more particles reflects itself not only in the absolute value of the square of the 
matrix element but also in distributions of variables determining the final state kinematics. Usually, in 
non-relativistic calculations of the total cross section, one takes the Jacobi momenta, choosing as inde- 
pendent variables the 'q'-meson momentum in the reaction centre-of-mass frame and 'k'-momentum 
of either nucleon in the rest frame of the nucleon-nucleon subsystem. By means of the A function they 
can be expressed as: 

yXMpE, and ^^VKsi2MM) (12) 
^ 2Vi 2^^ ' ^ ' 



with Si2 = |Pi -I- P2|^ denoting the square of the invariant mass of the nucleon-nucleon subsystem. In 
a non-relativistic approximation the expression of the total cross section defined by formula (5) for a 
meson production reaction in nucleon-nucleon interactions of the type NN NN Meson simplifies 
to: 



a oc 



/ kq^\Mab^i23\^dq. (13) 
Jo 



Denoting by L and I the relative angular momentum of the nucleon-nucleon pair and of the meson 
relative to the NN system, respectively and approximating the final state particles by non-distorted 
plane waves (case of non-interacting objects), whose radial parts ipi{q,r) are given by the spherical 
Bessel functions: 

^Kq,r)«jKqr)'^°(^|^, (14) 
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an expansion of the transition amplitude iM^^^p for the LI partial wave around qr = leads to 

iM^^l^ocq^'k^^ (15) 

Applying the above proportionality in equation (13) and solving the integral, yields the partial cross 
section 

« C^l^aT^' « v'm'^'^^\ (16) 

where tjm = Qmax/i^M with niM denoting the mass of the created meson Thus — at threshold — for 
the Ss partial wave the cross section for the NN NN Meson reaction should increase with energy of 
the fourth power of rjM- The dimensionlcss parameter rjM was introduced by Rosenfeld [5] as a variable 
allowing for the qualitative estimations of the final state partial waves involved in pion production. 
He argued that if the phenomenon of pion production takes place at a characteristic distance R from 
the centre of the collision, with R in the order of h/Mc, then the angular momentum of the produced 
meson is equal to / = Rq = hq/Mc. Hence, rjM denotes the classically calculated angular maximum 
momentum of the meson relative to the centre-of-reaction. The same arguments one finds in the work 
of Gell-Mann and Watson [4], where the authors do not expect the range of interaction to be greater 
than the Compton wave-length of the produced meson (h/Mc). However, it is rather a momentum 
transfer Ap between the colliding nucleons which determines the distance to the centre of the collision 
R ^ h/Ap at which the production occurs. Based on the data indications, it is emphasized in the 
original article of Rosenfeld [5] that R is slightly less than h/{2Mc), which is numerically close to the 
value of h/Ap. Directly at threshold, where all ejectiles are at rest in the centre-of-mass frame, Ap is 
equal to the centre-of-mass momentum of the interacting nucleons and hence, exploring equation (8) 
it can be expressed as: 



V K^th,T^l, ^ forpp^ppX ^ l_ _ ^n2 

2 a/s^ 



Aptft = z— = ^^^= \ nipuix + —. (17) 



Though the present considerations arc constrained to the spin averaged production only, it is worth 
noting that very close-to-threshold — due to the conservation laws and the Pauli excluding principle ^ 
— for many reactions there is only one possible angular momentum and spin orientation for the incoming 
and outgoing particles. 



Table 1: Partial wave transitions for the pp — > pp Meson and nn nn Meson reactions at threshold 



type 


meson 


spin and parity 


transition 


pseudoscalar 

vector 
scalar 


71,7], 7]' 
P,UJ,(f) 

do, fo 


0" 
1- 

0+ 


=^Po ^ 'Sos 

3Pi ^ %S 

'So ^ 'Sos 



For example the production of neutral mesons with negative parity — as pseudoscalar or vector 
mesons — may proceed in the proton-proton collision near threshold only via the transition between 

"'In former works [4,5] dealing only with pions this parameter is denoted by r], here in order to avoid ambiguities with 
the abbreviation for the eta-meson, we introduce an additional suffix M. 

^The Pauli principle for the nucleon-nucleon system implies that (— 1)L+S+T _ where L, S and T denote angular 
momentum, spin and isospin of the nucleon pair, respectively. For example if the nucleon-nucleon wave function is 
symmetric in space (L = 0) and spin (S = 1) then it must be antisymmetric in isospin-space (T = 0) to be totally 
antisymmetric. 
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and ^SqS partial waves (see table 1). Therefore only the coUision of protons with relative angular 
momentum equal to Ih and a parallel spin orientation — which flips during the reaction — may lead to 
the production of negative parity neutral mesons. In the case of the production of neutral scalar mesons, 
the protons or neutrons must collide with anti-parallel spin orientations which remains unchanged after 
the reaction. These simple considerations imply that in the close to threshold measurements with the 
polarized beam and target one should see a drastic effect in the reaction yield depending whether the 
polarization of reacting protons is parallel or anti-parallel. In the ideal case, having the product of beam 
and target polarization equal to unity, the close to threshold 7r°, 77 or rj' production can be realized only 
in case of the parallel spin orientation of the reacting protons and for the anti-parallel polarizations 
zero yield should be observed. Indeed, in reality, strong differences in the yield for various combinations 
of beam and target polarization have been determined in the pioneering measurements of the reaction 
pp^ ppn^ at the lUCF facility [38]. 



Due to both the strong nucleon-nucleon low-energy interaction in the So state and the mcson- 
nucleon forces in the exit channel, the assumption of the non-interactive plane waves leading to equa- 
tion (16) failed when confronted with experimental data. However, these deviations give the possibility 
for determining the still poorly known nucleon-meson interactions. 

15 




Figure 4: Monte-Carlo Simulations: (a) Phase-space distribution for the pp — > ppr] reaction at Q = 
16 MeV modified by the proton-proton final state interaction, (b) The dotted line shows the projection 
of the pure phase-space density distribution onto the Spp axis and the solid curve presents its modification 
by the proton-proton FSI. (c) The square of the scattering amplitude for the pp pp elastic scattering 
as a function of the proton-proton invariant mass in the range 2m.p < y^s^ < 2mp + Q, calculated 
according to formula (34) of section 3. (d) Phase-space density distribution modified by the proton-77 
interaction, with scattering length equal to a^^ = 0.7fm + i 0.3 fm. The proton-r^ scattering amplitude 
has been calculated according to the equation (37). A detailed descussion of the nucleon-nucleon and 
nucleon-meson interaction will be presented in section 3.1. 

The interaction between particles depends on their relative momenta. Consequently, for investiga- 
tions of final state interactions more instructive coordinates, than the q and k momenta, are the squared 
invariant masses of the two-body subsystems [39]. These are the coordinates of the Dalitz plot. In the 
original paper [40] Dalitz has proposed a representation for the energy partitions of three bodies in an 
equilateral triangle whose sides arc the axes of the ccntrc-of-mass energies. He took advantage of the 
fact that the sum of distances from a point within the triangle to its sides is a constant equal to the 
height. Therefore, the height of the triangle measures the total energy = + E2 + Eg and interior 
points — fulfilling four-momentum conservation — represent energy partitions. For a constant ^/s due 
to the energy conservation, without loosing any information, it is enough to consider the projection on 
any of the Ej E^ planes. The linear relation between E* and s^-^ (sjk = s + mf — 2 ^/sE*) allows to use 
Sjk Sifc or Ej Ej coordinates equivalently, with the following relation between the phase-space intervals: 
dE* dE* = ^ dsjk dski- According to the Kilian's geometrical representation [39] a Dalitz lies on a plane 
in the three dimensional space (si2, S13, S23) orthogonal to the space diagonal. The plane including a 
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Dalitz plot corresponding to a given total energy s is then given by the following scalar product [39] : 



(1, 1, l)(Si2,Si3,S23) =Si2+Si3+S23 = s + + + m|. (18) 

The second equality of equation (18) means that there are only two independent invariant masses of 
the three subsystems and therefore a projection onto any of the Sij Sjk planes still comprises the whole 
principally accessible information about the final state interaction of the three-particle system. In the 
case of no dynamics whatsoever and the absence of any final state interaction the occupation of the 
Dalitz plot would be fully homogeneous since each phase-space interval would be equally probable. The 
final state interaction would then appear as a structure on that area. Figure 4a shows — for the example 
of the pp — > ppr) reaction — how the uniformly populated phase-space density is modified by the S-wave 
(^So) interaction between outgoing protons. An enhancement in the range corresponding to the low 
relative momentum between protons is clearly visible. A steep decrease of the occupation density with 
increasing invariant mass of the proton-proton subsystem is even better seen in figure 4b presenting 
the projection of the phase-space distribution onto the Spp axis indicated by an arrow in figure 4a. This 
is a direct reflection of the shape of the proton-proton (^Sq) partial wave amplitude shown in figure 4c. 
Figure 4d shows the Dalitz plot distribution simulated when switching off the proton-proton interaction 
but accounting for the interaction between the ?7-meson and the proton. Due to the lower strength of 
this interaction the expected deviations from the uniform distributions are by about two orders of 
magnitude smaller, but one recognizes a slight enhancement of the density in the range of low invariant 
masses of proton-r^ subsystems. However, due to weak variations of the proton-ry scattering amplitude 
the enhancement originating from the T^-meson interaction with one proton is not separated from the 
r/-meson interaction with the second proton. Therefore an overlapping of broad structures occurs. It 
is observed that the occupation density grows slowly with increasing Spp opposite to the effects caused 
by the S-wave proton-proton interaction, yet similar to the modifications expected for the P-wave 
one [41]. From the above example it is obvious that only in high statistics experiments signals from the 
meson-nucleon interaction can appear over the overwhelming nucleon-nucleon final state interaction. 
It is worth noting, however, that the Dalitz plot does not reflect any possible correlations between the 
entrance and exit channel [39]. 

The Dalitz plot representation allows also for a simple interpretation of the kinematically available 
phase-space volume as an area of that plot. Namely, equation (5) becomes: 

^^"E^ir / ^^12 / rfS23 |M„b^i23|^, (19) 



F4s 

(mi+m2)2 sS™(si2) 



where the limits of integrations deflning the boundaries of the Dalitz plot can be expressed as [36]: 

sS"''(si2) = m^ + m^ - ^ I (si2 - s + m2)(si2 + m^ - m^) - ^A(si2, s, m2)A(si2, m^, m2)| , 

sS'" (^12) = + ml - ^ I (si2 - s + m^) (si2 + m^ - m^) + ^A(si2, s, mi)A(si2, m^, mf)| . (20) 

Thus, the phase-space volume kinematically available for the three-body flnal state can be written by 
means of only one integral: 

Vps^ J d^ps^-^ J ^^12 J ds23 = ^ J ~^\/ ^^(^12, s, m|)A(si2, m|, m^), (21) 

(mi+m2)2 s5J*"(si2) (mi+m2)2 
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whose solution leads, in general, to elliptic functions [36]. However, in the nonrelativistic approximation 
it has the following closed form: 



^mi m.2 ms 



(22) 



ps 



2 (mi + m2 + ms) 2 



where the substitution of the non-relativistic relation between 77^3 and Q 



m| + 2 m| (mi + m2) 




2 (mi + m.2) 



gives the Ss partial cross section of equation (16). On the basis of the formula (22) the kinematically 
available phase-space volume (V^s) can be as easily calculated as the excess energy Q. Close-to-threshold 

— in the range of a few tens of MeV — the non-relativistic approximation differs only by a few per 
cent from the full solution given in equation (21), which in fact with an up-to-date computer can be 
solved numerically with little effort. Therefore, in the following chapters we will describe the data as a 
function of Vps as well as of Q or 77^, if it is found to be appropriate. 

The choice of the proper observables for picking up the dynamical effects in the case of a four- 
particle exit channel is by far more complicated. It was for example discussed by Chodrow [42] who 
has given the extension of the Dalitz representation to the four-body final states. He has proven that 
the two-body correlation in the case where there are two identical particles in the four-body final state 

— say ij — can be directly visible in the density distribution of events on the Fj Ej-plane, where Ej 
denotes the energy of particle i and Fj is defined as: 



2.2 Scattering inside the deuteron 

Extending meson production experiments from the free nucleon-nucleon interaction to interactions of 
light nuclei is interesting for many reasons. First of all it offers the possibility to study more complex 
production processes with more than only one target nucleon involved. Furthermore, in the absence 
of free neutron targets, it gives the possibility to study the meson creation in the proton-neutron 
interaction at proton-beam facilities, complementary to a neutron-induced tt meson production with 
use of high-resolution (energy spread less than 1 MeV) secondary neutron beams [43-45] investigated 
at TRIUMF [46]. The technique based on quasi-free proton-nuclear reactions, with a beam of protons 
reacting with a nucleus in the target is utilized at the CELSIUS and COSY accelerators [47-51]. First 
experiments of 7r° meson production in the proton-neutron reaction with simultaneous tagging of the 
spectator proton have recently been carried out at CELSIUS and resulted in a resolution of the excess 
energy equal to a = 1.8 MeV [47]. In the quasi- free approximation (fig. 5) it is assumed that the 
bombarding proton interacts exclusively with one nucleon in the target nucleus and that the other 
nucleons affect the reaction by providing a momentum distribution to the struck constituent only. This 
assumption is justified if the kinetic energy of a projectile is large compared to the binding energy of the 
hit nucleus. In fact, as noticed by Slobodrian [52], also the scattering of protons on a hydrogen target, 
where the protons are bound by molecular forces, may serve as an extreme example of the quasi-free 
reaction. In that case, although the hydrogen atoms rotate or vibrate in the molecule, their velocities 
and binding forces are totally negligible with respect to the velocity and nuclear forces operating the 
scattering of the relativistic proton [52] . 




(23) 
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Figure 5: Spectator model for a particle production reaction via pd — > PspX. 
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Figure 6: (a) Momentum and (b) kinetic energy distribution of the nucleons in the deuteron, generated 
according to an analytical parametrization of the deuteron wave function [53, 54] calculated from the 
Paris potential [55]. 

The deuteron is also relatively weakly bound with a binding energy of 2.2 MeV, which is by 

far smaller — more than two orders of magnitude in the case of pion and already more than three 
orders of magnitude in the case of meson production — than the kinetic energy of the bombarding 
protons needed for the creation of mesons in the proton-neutron interaction. However, it has to be 
considered that even the low binding energy of ^ 2.2 MeV results in large Fermi momenta of the 
nucleons which can not be neglected. The momentum and kinetic energy spectra of the nucleons in the 
deuteron are shown in figure 6. 

In the considered approximation the Fermi motion of the nucleons is the only influence of the 
internucleon forces on the proton-neutron reaction and hence the struck neutron is treated as a free 
particle 

p f '>T' \ — P ^ Meson , . 

\ P J Psp 
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in the sense that the matrix element for quasi-free meson production off a bound neutron is identical to 
that for the free pn — > pn Meson reaction. The proton from the deuteron is considered as a spectator 
which does not interact with the bombarding particle, but rather escapes untouched and hits the detec- 
tors carrying the Fermi momentum possessed at the moment of the collision. From the measurement of 
the momentum vector of the spectator proton p^^ one can infer the momentum vector of the neutron 
p„ = —Psp the time of the reaction and hence calculate the excess energy Q for each event, provided 
that the beam momentum is known. As an example, a distribution of the excess energy in a quasi-free 
pn — >■ pnrj' reaction is presented in figure 7a [51]. Due to the large centre-of-mass velocity (/3 ~ 0.75) 
with respect to the colliding nucleons, a few MeV wide spectrum of the neutron kinetic energy inside a 
deuteron (fig. 6b) is spread by more than a factor of 30. 
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Figure 7: (a) Distribution of the excess energy Q for the pnr]' system originating from the reaction 
pd PspP^'V' calculated with a proton beam momentum of 3.350 GeV/c and the neutron momentum 
smeared out according to the Fermi distribution shown in figure 6a. (b) Spectrum of the ofi'-shell mass 
of the interacting neutron, as calculated under the assumption of the impulse approximation. 

Therefore, especially in the case of near-threshold measurements, where the cross section grows 
rapidly with increasing excess energy (see e.g. fig. 2) the total ccntre-of-mass energy ^/s has to be 
determined on an event-by-event level. For this purpose, in experiments, the spectator protons are 
usually registered by silicon pad- or /x-strip detectors [47,51,56] which determine their kinetic energy 
(Tsp) and polar emission angle (9). Thus, it is useful to express s as function of these variables: 



s = |Pp + P„|' = So - 2 Tsp (md + Ep) + 2 p^ ^T% + 2mpTsp cos(e) 



(25) 



with So denoting the centre-of-mass energy, assuming a vanishing Fermi motion. Measuring both the 
energy and the emission angle of the spectator protons it is possible to study the energy dependence of 
a meson production cross section from data taken at only one fixed beam momentum. 

It must be noted, however, that in the framework of the impulse approximation, illustrated in 
figure 5, the measured spectator proton is a physical particle, yet the reacting neutron is ofi^ its mass 
shell, where the explicit expression for its four-momentum vector in the rest frame of the deuteron, 
reads: 



(nid-mp-T,p, -p,p), 



(26) 
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with Tgp and p^^ denoting the kinetic energy and the momentum vector of a spectator proton, respec- 
tively. The mass spectrum of the interacting neutron (m^ = l^nP) resulting from the distribution of 
Fermi momentum is shown in figure 7b. It can be seen that the maximum of this spectrum differs only 
by about 3McV/c^ from the free neutron mass (m„ = 939.57 MeV/c^), however on the average it is off 
by about 9MeV/c^. In the frame of the discussed approximation, the struck neutron is never on its 
mass shell and the minimum deviation from the real mass occurs for vanishing Fermi-momentum and 
— as can be inferred from equation (26) — is equal to the binding energy = m^ — m„ — m^. Measure- 
ments performed at the CELSIUS and TRIUMF accelerators for the pp — > ppr] [50] and pp dn^ [57] 
reactions, respectively, have shown that within the statistical errors there is no difference between the 
total cross section of the free and quasi-free processes. In figure 8a the production of the rj meson in 
free proton-proton collisions is compared to the production inside a deuteron and in the overlapping 
regions the data agree within the statistical errors. These observations allow to anticipate that indeed 
the assumption of the identity for the transition matrix element for the meson production off free and 
quasi-free nucleons bound in the deuteron is correct. In case of the meson production off the deuteron, 
one can also geometrically justify the assumption of the quasi-free scattering since the average distance 
between the proton and the neutron is in the order of ^ 3 fm. Of course, the other nucleon may scatter 
the incoming proton and the outgoing meson. However, these nuclear phenomena are rather of minor 
importance in case of the production on the neutron bound in the deuteron, but should be taken into 
account for derivations of total cross sections from experimental data. 

The reduction of the beam flux on a neutron, due to the presence of the proton, referred to as 
a shadow effect, decreases for example the total cross section by about 4.5 % [59] for the ry-meson 
production. Similarly, the reduction of the total cross section due to the reabsorption of the outgoing 
T] meson on the spectator proton was found to be only about 3% [59]. The appraisals were performed 
according to a formula derived in reference [60] which shows that the cross section for the deuteron 
reduces by a factor of: 



R=l- a^^' < r-2 > /47r (27) 



compared to the free nucleon cross sections. Here cr^'j^' denotes the rjN inelastic cross section and 
< r~^ > stands for the average of the inverse square nucleon separation in the deuteron taking the 
nucleon size into account [59]. The latter effect for the production of mesons like (tt, uj, t]', (p) is expected 
to be much smaller, since (as shall be discussed in section 3) the s-wave interaction of the 77-meson 
with nucleons is by far stronger than for any of the mentioned ones. 



^The matter radius of the deuteron amounts to ~ 2fni [58]. 
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Figure 8: (a) Total cross sections for the pprj reaction as a function of excess energy measured for 
the scattering of protons in vacuum (open symbols) [61-66] and inside a deuteron (filled symbols) [50]. 
(b) Close-to-threshold total cross sections for the reactions pn —>■ dn^ (squares [44], diamonds [67]) and 
pd —>■ pdn^ (circles [68], triangles [69]). The solid curve illustrates the predictions for the total cross 
section of the pd pdn^ reaction assuming that the process is utterly due to the quasi-free pn —>■ dir^ 
reaction and taking into account the final state interaction between proton and deuteron. Dashed lines 
shows the result when the final state interaction was neglected. In the calculations the cross section of 
the pn — > dn'^ was parametrized as presented by the dotted line. 



Due to the relatively large distance between the nucleons in the deuteron one could argue that the 
quasi-free reactions will dominate the production of mesons even with a deuteron in the final state. 
The conviction that the production of mesons in nucleon-deuteron collisions proceeds predominantly 
via an elementary scattering of nucleons is also based on the comparative investigations of the 7r°- 
meson production in proton-neutron and proton-deuteron collisions [70]. The total cross section for 
the pd — s> pdn^ reaction close to the kinematical threshold is by more than two orders of magnitude lower 
compared to the elementary pn — > dn^ process (see figure 8b). This observation has been successfully 
explained assuming that the production of the vr*^ meson in the pd pdn^ reaction is proceeding 
entirely due to the quasi-free pn dir^ reaction [70] for tjm larger than 0.7 (dashed line). In the 
threshold region the result of the calculation (solid line) describes exactly the experimental data when 
the proton-deuteron FSI is additionally taken into account. The large difference in magnitude between 
the threshold cross section for the pn —>■ dn^ and the pd — > pdn^ reaction can be explained since 
the incident proton momentum needed to create a vr" meson via the pd pdn^ process is much 
below the kinematical threshold for the production in elementary proton-neutron collisions. Therefore 
only nucleons possessing large Fermi momenta can contribute to the vr" production and this appears 
with a small probability (compare fig. 6a). The other possible elementary processes pp — >■ ppir^ and 
pn pniT^ can be neglected since they are expected to lead predominantly to a four-body exit channel. 
This explanation could serve as a heuristic example, however recent data [69] on the cross sections of 
dp dpn^ show that the supposed spectator protons to the elementary np dir^ reaction are observed 
in phase-space regions corresponding to zero Fermi momenta and hence the quasi-free production could 
be significant only at energies close to the NN — > NNn threshold {rjM {pd — > pdn^) ~ 0.83). In this 
experiment the deuteron beam was impinging onto the hydrogen target, and hence the spectator protons 
were boosted in the laboratory permitting registration also of the low-energy part of the spectators' 
spectral distribution. 
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A natural extension of the meson production in quasi-free proton-neutron collisions with either a 
deuteron used as the target or as the beam is the combination of both methods into a double quasi-free 
production occurring during the collisions of deuterons as depicted in figure 9. 



Figure 9: Schematic depiction of the dou- 
ble quasi-free nn — > nnX reaction. Dur- 
ing the coUision of deuterons (left hand 
side of the figure, with the total momen- 
tum (solid arrow) resulting from the sum 
of the beam momentum (dotted arrow) 
plus the Fermi momentum (short arrow)) 
a double quasi-free neutron-neutron re- 
action may lead to the creation of mesons 
(small gray circle). The spectator pro- 
tons (black circles) leave the reaction re- 
gion with their initial momentum plus the 
Fermi momentum, which they possessed 
at the moment of the reaction. Neu- 
trons are plotted as open circles. Due 
to the large relative momenta between 
spectators and the outgoing neutrons (~ 
IGeV/c close to the threshold for the 
77 meson production) a distortion of the 
nnX system by the accompanying pro- 
tons can be neglected. 

Utilizing this technique [71] and registering both the slow and fast spectator protons could allow 
for the study of meson production in neutron-neutron collisions with a very precise determination of 
the excess energy which depends on the accuracy for the registration of the momentum vectors of the 
spectators. 



3 Phenomenology of the initial and final state interaction 

The simplest method to study the interactions between meson and nucleon is to scatter a beam of 
mesons off hydrogen or deuteron targets and to measure the behaviour of the total and differential 
cross section in the clastic scattering of the considered particles. Figure 10a presents the total cross 
section for K^p and K^p elastic scattering as a function of the centre-of-mass excess energy. When 
comparing the data to calculations [72] — including the changes of phase-space integral (eq. (10)) and 
Coulomb interaction in the initial and final states — one observes a huge enhancement for the K~p 
cross section with decreasing excess energy and a slight suppression in the case of K^p scattering. This 
observation may be attributed to the slight repulsion due to the kaon-proton hadronic interaction and 
the significantly larger scattering caused by the strong interaction between the K~ and proton due to 
the vicinity of the A(1405) hyperon resonance. The effect has a direct influence on the effective mass 
splitting of kaons immersed in dense nuclear matter [73,74]. 
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Figure 10: (a) Cross section for the X^-proton [75] and X^-proton [76] elastic scattering. Solid and 
dashed lines represent the changes of the phase-space integral (eq. (10)) modified by the initial and 
final state Coulomb interaction {crxp^Kp = const ■ C^/ (IGvrs)). Both curves are normalized to points 
of large excess energies, (b) Total cross section of the pp ppK^K~ reaction [77,78]. Statistical and 
systematical errors are separated by dashes. The solid line indicating calculations of reference [79] is 
described in the text, (c) Decomposition of the total cross section of the pp — > pp7r° reaction into Ss, Ps 
and Pp final state angular momenta. The dashed and solid lines represent the r]\j and ryf^ dependence 
of Ps and Pp partial cross sections, respectively. The remainder is indicated as the dotted line. The 
Sp partial wave is forbidden by the conservation laws and the Pauli excluding principle. Note that at 
r]M — 1 the Pp and Ps partial waves seem to dominate. However, the analysis of the differential cross 
section measured at CELSIUS at tim — 0.449 [80, 81] showed that also a d-wave pion production — 
due to the interference between Ss and Sd states — constitutes 7 % of the total cross section, when a 
meson-exchange model is assumed. The figure has been adapted from [82]. 



A Coulomb repulsion was taken into account by multiplying the expression (10) by the Coulomb 
penetration factor for both exit and entrance channels. determines the ratio of the probability 
of finding two particles close together to the probability of finding two uncharged particles together, all 
other things being equal [83] and can be expressed as [84] : 

where r]c is the relativistic Coulomb parameter, which for the collision of particles i,j reads: 



Qj q j a 



= qj q,- a 



m? 



m; 



A(sij,m2,m2) 



with the fine structure constant a, the relative velocity v of the colliding particles and with q^, q^ 
denoting their charges ^. The comparison of the K~y> scattering data with the phase space integral 
presented in figure 10a aims only to illustrate the qualitative influence of the meson-nucleon dynamics 
on the total cross section energy dependence. It is worth noting, however, that at present the data can 
be qualitatively well described in the frame of the effective chiral Lagrangian with a coupled channel 
potential [26,87,88]. 

''For collisions at an angular momentum / larger than h the of equation (28) needs to be multiplied by a factor 
ofnLi(l + (Wr^)') [85,86]. 
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Another method permitting such investigations is the production of a meson in the nucleon-nucleon 
interaction close to the kinematical threshold where the outgoing particles possess small relative ve- 
locities. The experimental observables are the excitation functions and the double differential cross 
section. A possible deviation from the expectation based on the assumption of homogeneous phase- 
space abundance corrected for the known nucleon-nucleon interaction delivers information about the 
meson-nucleon forces. For example, figure 10b presents the first data points on the close-to-threshold 
production of the K~ meson via the pp ppK^K~ reaction. The superimposed solid line depicts the 
results of calculations [79] taking into account the changes of the production amplitude as deduced from 
the K^p and K^p elastic scattering shown in figure 10a, but neglecting the influence of the dominant 
proton-proton interaction! On the other hand the proton-proton FSl can not be neglected in case of 
three-body final states, e.g. pp pp-n^ [8,89], pp ppr] [61,63], or pp pprj' [37,90], since it has 
strong infiuence on the total cross section energy dependence enhancing it by more than an order of 
magnitude for excess energies below Q ~ 15McV as shall be shown below. Thus it is surprising that 
in spite of its ncglcction one can describe the data of the pp — > ppK^K^ reaction. The origin of that 
effect will be investigated experimentally in the near future [91]. At present one can only speculate 
whether it is due to the partial compensation of the pp and K'p hadronic interaction or maybe due to 
the additional degree of freedom given by the four-body final state [92] . 




Figure 11: (a) The ''^So and ^Pq phase-shifts of the nucleon-nucleon potential shown versus the centre- 
of-mass kinetic energy available in the proton-proton system. The values have been extracted from the 
SAID data base [93] (solution SM97). For higher energies the S- and P-wave phase-shifts arc nearly the 
same. This is because the collision parameter required to yield the angular momentum oilh diminishes 
significantly below 1 fm with increasing energy and consequently the interaction of nucleons — objects 
of about 1 fm size — becomes almost central, (b) The variable tim as a function of the excess energy for 
7r°, r] and 77' mesons produced via pp — > pp Meson reactions, (c) The phase-space volume Vpg (defined 
by equation (21)) versus the excess energy Q. The picture indicates that for the production of "heavy 
mesons" in the nucleon-nucleon interaction at a given Q value there is only a slight difference of the 
available phase-space volume on the produced meson mass, which is larger than that of 7r° production 
by about 30% only. Therefore, for the comparative studies of the production dynamics of different 
mesons, Q is as much a suitable variable as Vps- 

In order to demonstrate the influence of the nucleon-nucleon interaction on the cross section de- 
pendence on the excess energy let us consider the production of neutral pseudoscalar mesons (vr", 77, 
T]') in the collision of protons. Due to the short life-time of these mesons the study of their low-energy 
interaction with nucleons in direct scattering experiments is not feasible. In general for a three-body 
exit channel one expects a dependence of the total cross section on the energy which can be described 
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by the linear combination of partial cross sections of equation (16). Therefore, to be able to extract the 
information about the final state interaction of the outgoing particles from the energy dependence it is 
necessary to know precisely the contribution originating from different partial waves. Yet appropriately 
close-to-threshold there is only one important combination of the angular momenta of emitted parti- 
cles (Ss) and in this region the energy dependence of the total cross section is uniquely determined. 
However, the range of excess energies for the S-wave dominance changes strongly with the mass of the 
produced meson. 

For the production of the 7r° meson in proton-proton collisions the investigations with polarized 
beam and targets [82,94] allowed to deduce that the Ss partial-wave accounts for more than 95 % of the 
total cross section up to r]M ~ 0.4, as can be seen in figure 10c, where the Ss contribution is indicated by 
the dotted line. The Ss contribution was inferred assuming the rj^ and rj^ dependence for Ps and Pp 
partial waves, respectively. These are power-laws taken from proportionality (16), which was derived 
under the assumption of non-interacting particles. Relatively small values of ^Po-wave nucleon-nucleon 
phase-shifts at low energies (compared to ^Sq phase-shifts in figure 11a) and similarly weak low-energy 
interactions of P-wave protons in other spin combinations [95] justify this assumption. 

The measurements with the polarized beam and target allow a model independent determination of 
the contribution from individual partial waves. In case of the pp ppn^ reaction, due to the identical 
particles in the initial state and the rotational equivalences [82, 96] there are only seven independent 
polarization observables, e.g. two beam analyzing powers and five linear combination of spin correlation 
coefficients. 

In particular, as shown in references [82,94] the close to threshold contributions of the Ps and Pp 
partial waves can be determined in the model-free way from the measurements of the spin dependent 
total cross sections only. For example the strength of the Ps final state can be expressed as [82] : 



where atot denotes the total unpolarized cross section and Aut and Acr^ stand for differences between the 
total cross sections measured with anti-parallel and parallel beam and target polarizations. Subscripts 
T and L associate the measurements with the transverse and longitudinal polarizations, respectively. 

In accordance with the phenomenology of Gell-Mann and Watson [4] described in the previous 
section one expects that also in the case of heavier mesons the Ss partial wave combination will constitute 
the overwhelming fraction of the total production cross section for rjM smaller than 0.4. This implies 
— as can be deduced from the relation between tim and Q illustrated in figure lib — that mesons 
heavier than the pion are produced exclusively via the Ss state in a much larger excess energy range 
and hence larger phase-space volume (see figure 11c). In particular in case of rj' this is larger by more 
than one order of magnitude. Thus, whereas in case of tt^ the onset of higher partial waves is observed 
at Q around lOMeV it is expected only above 100 McV and above ~ 40McV for rj' and r] mesons, 
respectively. Figures 12b and 12c present the angular distributions of the created rj meson in proton 
collisions. It is evident that at Q = 15.5 MeV and still at Q = 41 MeV the production of the rj meson 
is completely isotropic within the shown statistical errors. Altliougli at Q = 41 MeV the accuracy of 
the data docs not exclude a few per cent of the contribution originating from higher partial waves, 
the dominance of the s-wave creation is evident. Similarly, the measurements of the differential cross 
section (fig. 12a) for the pp — > pprj' reaction performed at SATURNE [97] at Q = 143.8 MeV are still 
consistent with pure Ss-wave production, though the relatively large error bars would allow for other 
contributions on a few per cent level {^^ 10%). 
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Figure 12: Differential cross section of the pp — > pp Meson reaction as a function of the meson centre-of- 
mass polar angle. Dashed lines indicate the isotropic distribution. Shown are results of measurements 
for the pp — > pprj' reaction taken at Q = 143.8 MeV [97] (a) and for the pp — > ppr] reaction at Q = 
15.5 MeV [66] (b) and Q = 41 MeV [98] (c). In all pictures only statistical errors are plotted, which in 
figure (b) are smaller than the size of symbols. The distribution presented in picture (b) is consistent 
with a measurement performed at an excess energy of Q = 16 MeV by means of the PROMICE/WASA 
detector [99], whereas the data at Q = 37 MeV also from reference [99] deviate significantly from 
isotropy However, data shown in picture (c) have been taken with a detector of much higher angular 
acceptance. 



Table 2: Momentum transfer Ap calculated according 


to equation 


17 and the corresponding distance 


R ^ h/Ap probed by the NN 


NN Boson reaction at the kinematical threshold for different particles 


produced. The table has been 


adapted from [20]. 






particle 


mass [MeV] 


Ap [fm 


^] R [fm] 


7 





0.0 


oo 


TT 


140 


1.9 


0.53 


V 


550 


3.9 


0.26 


p,uj 


780 


4.8 


0.21 


v' 


960 


5.4 


0.19 





1020 


5.6 


0.18 



Let us now consider to what extent the energy dependence of the total cross section in the estimated 
range of the dominance of the Ss partial waves can be understood in terms of the phase-space variation 
and the interaction between the particles participating in the reaction. Watson [18] and Migdal [100] 
have proposed the factorization of the amplitude in the case where the production is of short range and 
the interaction among the outgoing particles is of long range. This requirement is well fulfilled in the 
case of the close-to-threshold meson production due to the large momentum transfer (Ap) between the 
interacting nucleons needed to create the considered mesons (tt, r^, . . . ,0). According to the Heisenberg 
uncertainty relation the large momentum transfer brings about a small space in which the primary 
creation of the meson takes place. In table 2 the distance probed by the NN — ^ A^A^ Meson reaction at 
threshold is listed for particular mesons. It ranges from 0.53 fm for pion production to 0.18 fm for the (j) 
meson, whereas the typical range of the strong nucleon-nucleon interaction at low energies determined 
by the pion exchange may exceed a distance of few Fermi and hence is by one order of magnitude larger 
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than the values hsted in table 2. Thus in analogy to the Watson-Migdal approximation for two-body 
processes [18] the complete transition matrix element of equation (5) may be factorized approximately 
as ^ 

|M„p^|2 ^ \Mfsi? ■ |Mo|2 • Fisi, (30) 

where Mq represents the total short range production amplitude, Mpsi describes the elastic interaction 
among particles in the exit channel and Fjsi denotes the reduction factor accounting for the interaction 
of the colliding protons. Further, in the first order approximation one assumes that the particles 
are produced on their mass shell and that the created meson does not interact with nucleons. This 
assumption implies that the iMpsil"^ term can be substituted by the square of the on-shell amplitude 
of the nucleon-nucleon elastic scattering: 

\Mfsi\'^ = \Mnis[^nn\'^- (31) 

Effects of this rather strong assumption will be considered later, while comparing the estimation with 
the experimental data. 

In the frame of the optical potential model the scattering amplitude is determined by phase-shifts. 
Particularly, the ^Sq proton-proton partial wave — relevant for further consideration — can be expressed 
exphcitly as follows [101]: 



where C denotes the square root of the Coulomb penetration factor defined by equation (28), k stands 
for either proton momentum in the proton-proton rest frame and the phase-shift is indicated by Spp. 
The phase-shifts 5pp(^So) can be extracted from the SAID data base (see fig. 11a) or, alternatively, can 
be calculated according to the modified Cini-Pubini-StangheUini formula including the Wong-Noyes 
Coulomb correction [102-104] 

k ctg(5pp) + 2krj^ h{rj,) = + ^ - (33) 

where h{T]^ = -ln{r]c) - 0.57721 + rj'^ J2n=i niJ+vI) t^^l' '^'^ ^ ™ 

The phenomenological quantities Upp = — 7.83fm and bpp = 2.8 fm denote the scattering length and 
effective range [102], respectively. The parameters Ppp — 0.73 fm^ and Qpp — 3.35 fm^ are related to the 
detailed shape of the nuclear potential and derived from a one-pion-exchange model [102] . Substituting 
equation (33) into equation (32) allows to calculate the low-energy amplitude for the proton-proton 
elastic scattering ^: 



I' = -— — r-. — rrr^^^^TT^ ^- ^^4) 



The result is presented as a solid line in figure 13a and is in good agreement with the values obtained 
from the phase-shifts of the VPI partial wave analysis [93], shown as solid circles and with the phase- 
shifts of the Nijmegen analysis [105], shown as open squares. 

^For a comprehensive discussion of the FSI and ISI issue including a historical overview and a criticism of various 
approaches the reader is referred to [22] . 

^In principle the formula is valid for k < 133 MeV/c [103]. \Mpp^pp\'^ is taken to be constant for larger values of k. 
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Figure 13: (a) Square of the proton-proton scattering amplitude versus k, the proton momentum in 
the proton-proton subsystem from [101,102] (sohd hne), [106] (dashed hne), and [107,108] (dotted 
hnc). The filled circles arc extracted from [93] and the opened squares from [105]. The curves and 
symbols are arbitrarily normalized to be equal at maximum to the result from reference [106], shown 
as the dashed line, (b) Total cross section for the pp — > ppn*^ reaction as a function of the centre- 
of-mass excess energy Q. Data are from refs. [8,81,89,109]. The dashed hne indicates a phase-space 
integral normalized arbitrarily. The phase-space distribution with inclusion of proton-proton strong 
and Coulomb interactions fitted to the data at low excess energies is shown as the solid line. The dotted 
line indicates the parametrization of reference [110] written explicitly in equation (35), with e = 0.3. 



The factor is always less than unity due to the Coulomb repulsion between protons. At higher 
energies, where is close to unity, the nuclear scattering will be predominant and for the very low 
energies the Coulomb and nuclear interactions are competing. The Coulomb scattering dominates 
approximately up to about 0.8 MeV of the proton energy in the rest frame of the other proton, where 

equals to one-half [83]. In the case of the pp — > pp Meson reaction the maximum possible energy of a 
proton seen from another proton is equal to 0.8 MeV already at an excess energy of about Q = 0.4 MeV. 




Figure 14: (a) Total cross section for the pp pprj' reaction as a function of the centre-of-mass excess 
energy Q. Data are from refs. [62,90,97,111,112]. Statistical and systematical errors are separated by 
dashes. The solid line shows the phase-space distribution with inclusion of proton-proton strong and 
Coulomb interactions. The dotted line indicates the parametrization of reference [110] written explicitely 
in equation (35), with e = 0.3 and the dashed line indicates a phase-space integral normalized arbitrarily, 
(b) Stars represent the data of the pp — > pK'^A reaction [113-119]. The dashed line presents the energy 
dependence defined by the phase-space. The sohd line shows the calculation with proton-ir+ Coulomb 
repulsion and a proton-A strong interaction taken into account as described in the text. 
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Therefore, a significant influence of tfie Coulomb repulsion on the energy dependence of the total 
production cross section is expected only at very low excess energies, i.e. conservatively for Q < 2 MeV. 

Assuming that the on-shell proton-proton amplitude exclusively determines the phase-space pop- 
ulation one can obtain the total cross section energy dependence substituting equation (34) into the 
formula (19). Solid lines in figures 13b and 14a represent the determined dependence for the pp ppTi^ 
and pp — > ppt]' reactions, respectively. The absolute scale of the calculations has been fixed by the nor- 
malization to the data. One recognizes the good agreement with the experimental points in the excess 
energy range up to Q 10 MeV for 7r° meson production and up to Q 150 MeV for rj'. This confirms 
the earlier considerations that the domain of Ss partial wave dominance should widen significantly with 
the mass of the created meson and as expected it extends above 100 MeV in case of the rj' meson. 
Comparing the data to the arbitrarily normalized phase-space integral of equation (21) reveals that the 
proton-proton FSI enhanced the total cross section by more than an order of magnitude for low excess 
energies. This is particularly clear in figure 14a where the logarithmic scale is utilized. However, the 
much weaker nucleon-hyperon interaction increases the total cross section for the close-to-threshold 
meson production only by about a factor of two. Again in that case — as shown in figure 14b — the 
factorization of the matrix element into constant on-shell primary production amplitude and the elastic 
scattering between outgoing baryons leads to the good description of the close-to-threshold data. In 
this case the Coulomb interaction between K'^ meson and proton was taken into account multiplying 
the scattering amplitude by the Coulomb penetration factor. 

The dotted line in figures 13b and 14a, which for Q < 50 MeV is practically indistinguishable from 
the sohd fine, presents the excess energy dependence of the total cross section taking into account the 

proton-proton FSI effects according to the model developed by Faldt and Wilkin [110,120]. Representing 
the scattering wave function in terms of bound state wave function the authors derived a closed formula 
which describes the effects of the nucleon-nucleon FSI as a function of the excess energy Q only. This 
approach is specifically useful for the description of the spin-triplet proton-neutron FSI, due to the 
existence of a bound state (deuteron) with the same quantum numbers. Though the bound state of the 
proton-proton system does not exist, the model allows also to express the total cross section energy 
dependence for a pp ^ pp Meson reaction by a simple and easily utilizable formula: 




where the parameter e has to be settled from the data. The fiux factor F and the phase-space volume 
Vps are given by equations (9) and (22), respectively. The normalization can be determined from the 
fit of the data which must be performed for each reaction separately. 
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Figure 15: (a) Total cross section for the pprj reaction as a function of the centre-of-mass 

excess energy Q. Data are from refs. [50,61-66]. Statistical and systematical errors are separated by 
dashes. The dashed hne indicates a phase-space integral normahzed arbitrarily. The phase-space 
distribution with inclusion of proton-proton strong and Coulomb interactions fitted to the data in the 
excess energy range between 15 and 40McV is shown as the solid line. Additional inclusion of the 
proton-r/ interaction is indicated by the dotted line. The scattering length = 0.7fm + i0.4fm 
and the effective range parameter bp^ = — 1.50fm — i0.24fm [125] have been arbitrarily chosen for 
the calculations described in the text, (b) Dalitz plot distribution of the pp pprj reaction at an 
excess energy of Q = 37.6 MeV, corrected for the detection acceptance. Out of the two invariant 
masses corresponding to two p — r) pairs the one being larger is plotted along the x-axis. The figure 
is taken from [63] (A similar spectrum for the ppn^ system [89] reveals the influence of the proton- 
proton interaction, yet again, as for the total cross section energy dependence (fig. 13b), the proton-7r° 
interaction is too weak to affect observably the density distribution of the Dahtz plot.), (c) Dahtz plot 
distribution of the pp — > pprj reaction at Q = 15.5 MeV, corrected for the detection acceptance and the 
proton-proton FSI [66]. The proton-proton FSI enhancement factor was calculated as square of the 
on-shell proton-proton scattering amplitude as written in equation (34). 

In the case of 7] meson production the interaction between nuclcons is evidently not sufficient to 
describe the increase of the total cross section for very low and very high excess energies, as can be 
seen in figure 15a. The solid line in the figure was normalized to the data at an excess energy range 
between 15 McV and 40 MeV. The increase of the total cross section for higher energies can be assigned 
to the outset of higher partial waves. As expected from the previous considerations, this is indeed seen 
at Q f5i 40 MeV where the energy dependence of the total cross section starts to change its shape. On 
the contrary, the close-to-threshold enhancement — being by about factor of two larger than in the 
case of the tt^ and rj' mesons — can be assigned neither to the contribution from other than Ss partial 
waves nor to the variation of the primary production amplitude \Mq\. The latter is expected to change 
at most by a few per cent for excess energies below 20 MeV [121]. Instead, this discrepancy can be 
plausibly explained by the influence of the attractive interaction between the t] meson and the proton. 
Note, that the real part of the scattering length of the rj-pioton potential — depending on the analysis 
method and the studied reaction — is 3 to 10 times [122, 123] larger than the scattering length for 
7r°-proton scattering {ap^^ ?s 0.13 fm) [124]. Hence, the modiflcations of the total cross section energy 
dependence due to vr^-proton interaction are too weak to be observed within the up-to-date accuracy 
of measurements and calculations. However, the influence of the j^-proton interaction presented in 
figure 15 is evident and hereafter it will be considered whether it may serve for the estimation of the 
77-proton scattering parameters. 

The strict quantitative calculation requires, however, the evaluation of the three-body Faddeev 
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equation Here, we will rather present a simple phenomenological treatment which shall lead to 
the qualitative understanding how the mutual interaction among three outgoing particles affects the 
total cross section dependence on excess energy. One of the simplest possibilities based on the naive 
probabilistic interpretation of the incoherent pairwise interaction would be to factorize the overall 
enhancement factor into corresponding pair interactions [34,127]: 
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where |Mjj^jjp denotes the square of the elastic scattering amplitude of particles i and j. The |Mpp^ppp 
can be evaluated according to the formula (34), which for the s-wave 77-proton scattering, after substi- 
tution of = 1 and rjc — Ppp — 0, reduces to 
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denotes the 77 momentum in the proton^ rest frame. The scattering length and effective range b. 



•^prj 



are complex variables with the imaginary part responsible e.g. for the pr] — > pn conversion. 

The factorization of both i) the overall production matrix element (eq. (30)) and ii) the three 
particle final state interactions (eq. (36)) applied to the formula (19) gives the following expression for 
the total cross section of the pp — > ppr] reaction: 
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where the protons are distinguished by subscripts. Exploring formulas (34) and (37) gives the results 
shown as the dotted line in figure 15 a. Evidently, the inclusion of the proton-77 interaction enhances 
the total cross section close-to-threshold by about a factor of 1.5 and leads to a better description 
of the data. The effect of the proton-r/ interaction is also seen in the experimental distribution of 
the differential cross sections a / {d Sp-^Tj dsp^j)) shown in figures 15b and 15c. These distributions 
originating from kinematically complete measurements comprise the whole experimentally available 
information about the interactions of the ppTy-system. In figure 15b one recognizes the increase of the 
distribution density at regions of small invariant masses of the proton-proton and proton-jy subsystems. 
At this excess energy (Q = 37.6 MeV) these regions are quite well separated. However, since this is 
close to the energy where the advent of higher partial waves is awaited, the possible contribution from 
the P-wave proton-proton interaction can not be a priori excluded. Specifically, the latter leads to 
the enhancement at large invariant masses of the proton-proton pair [41] and hence affects the phase- 
space region where the modification from the proton-77 interaction is expected. Figure 15c presents 



^"An exact derivation of the Faddeev equation can be found for example in [126]. 

^^Note that the sign of the term — 1/a from equation (34) was changed because the imaginary part of the proton—?; 
scattering length is positive [125], as we also adopted here, whereas in the majority of works concerning nucleon-nucleon 
interaction, the scattering length is negative [31]. We are grateful for this remark to A. Gasparyan. 
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the high statistics data taken at an excess energy of Q = 15.5 MeV, where the assumption of the total 
dominance of the Ss partial wave is rather save. At this excess energy, corresponding to the small 
relative momentum range k^"^ ~ 105MeV/c, the variations of the proton-77 scattering amplitude are 
quite moderate (fig. 18c) and as expected from simulations presented in the previous section, the increase 
of the phase-space population due to the ?7-meson interaction with one proton is not separated from the 
other one. Interestingly, the density growth with increasing invariant mass of the proton-proton system 
is much faster than expected from the simulations presented in figure 4d, which have been performed 
with a scattering length equal to a^^ = 0.7 fm + zO.Sfm. However, when reducing the proton-proton 
FSl effect to a multiplicative factor, one finds that it depends on the assumed nucleon-nucleon potential 
and on the produced meson mass [128]. 

This issue was recently vigorously investigated e.g. by authors of references [20-22, 128, 129] and 
we shall briefly report this here as well. Up to now we factorized the transition matrix element into 
a primary production of particles and its on-shell rescattering in the exit channel (eqs. (30) (31) (36)). 
Though it is a crude approximation, neglecting the off-shell effects of the production process completely, 
it astoundingly leads to a good description of the energy dependence of the total cross section, as already 
demonstrated in figures 13b, 14a,b and 15a. The off-shell effects, as pointed out by Kleefeld [22], could 
have been safely neglected in case of the electromagnetic transitions in atoms or (5 decays, where the 
excitation energy of the involved nucleons is by many orders of magnitude smaller than their masses 
and the initial and final states go hardly off-shell [22]. However, in case of the NN — > NN Meson 
process the large excitation energy of the colliding nucleons is comparable with the nucleon masses and 
the primary interaction may create the particles significantly far from their physical masses, so that the 
off-shell effects can not be a priori disregarded. 




Figure 16: Diagrammatic representation of the DWBA expressed by equation (39). (a) The primary 
production term, (b) The loop diagram including the nucleon-nucleon FSI. T'(k',k) stands for the 
half-off-shell T matrix with k and k' denoting the on- and off-shell centre-of-mass momentum in the 
nucleon-nucleon system, respectively, q indicates the momentum of the created meson in the reaction 
centre-of-mass and p the momentum of the colliding nucleons. 

Generally, the decomposition of the total production amplitude into the primary production and 
the subsequent nucleon-nucleon interaction visualized in figure 16 is expressed by the formula: 

M = Mo"" + M^"GTnn: (39) 

where the second term of the equation represents the integration over the intermediate (k') momenta 
of the off-shell production amplitude and the half-off-shell nucleon-nucleon T matrix [128]. Assuming 
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that the primary production occurs in such a way that one nucleon emits the meson which then re- 
scatters on the other nucleon and appears as a real particle the authors of reference [128] found that the 
enhancement of the cross section due to the nucleon-nucleon interaction depends strongly on the mass 
of the created meson. This is because with the increasing mass of the produced meson the distance 
probed by the nucleon-nucleon interaction decreases (sec table 2) and hence the relevant range of the 
off-shell momenta becomes larger. The effect for the pp pp Meson reactions is presented in figure 17a, 
where one can see that, when utilizing the Bonn potential model for the nucleon-nucleon T matrix, 
the enhancement in case of the 7r° production is by about a factor of four larger than for the r] or rj' 
mesons. A similar conclusion, but with the absolute values larger by about 40 %, was drawn for the 
Paris NN potential [128]. On the contrary, when applying the Yamaguchi potential into calculations the 
enhancement grows with the increasing mass of the meson, as shown in figure 17b. Thick solid curves 
in figures 17a and 17b show the results of the frequently applied approximation of the nucleon-nucleon 
FSI effects: 

M = Mo"" + M^^^GTnn ~ M°"-(l + GTnn) = M^T\-k) = F^Nik), (40) 

where the overall transition matrix element M is factorized to the primary on-shell production and the 
A^A^ FSI expressed as the inverse of the Jost function J~^(— k) [108]. As can be seen in figures 17a 
and 17b, the variation of the absolute values — of such obtained enhancement factors — with the 
applied potential is severe. 




k[MeV/c] k[MeV/c] b) k[MeV/c] 



Figure 17: The FSI factor for Bonn [131] a) and Yamaguchi [132] b) potentials. The solid, dotted and 
dashed lines correspond to 7r°, 77 and 77' meson production, respectively. Thick solid lines indicate the 
inverse of the squared Jost function (|Fpp(k)p = | J(— k)|~^). c) The FSI factors for 7r° production of 
Paris (dotted curve) [55] , Bonn (dashed curve) and Yamaguchi (solid curve) potentials normalized to be 
equal at maximum to the inverse of the squared Jost function of the Yamaguchi potential (thick solid 
line). The shapes stemming from different potentials are almost indistinguishable. Note that the thick 
solid fine corresponds to the dashed fine in figure 13a. The figures are adapted from reference [128]. 

Since the physical value of the total cross section can not depend on the off-shell features of the 
potential used in calculations, which are in principal not measurable [130], the differences in the mag- 
nitude of the \Fpp\ factor must reflect itself in a corresponding dependence of the primary production 
amplitude on the potential used. Therefore, it is of crucial importance to realize that the values of 
the threshold amplitudes \Mq\ are significant only in the context of the potential they were extracted 
from. Figure 17c demonstrates, however, that the shapes of the enhancement factors, with the meson 
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exchange mechanism assumed for the primary production, are pretty much the same, independently 
of the apphed NN potential and correspond to the form of the Jost function inferred from the Yam- 
aguchi potential. This indicates that the energy dependence of the NN FSI factors is predominantly 
determined by the on-shell A^A^ T matrix. 

In references [20-22] the formula for the transition matrix element with explicit dependence on the 
considered off-shell features for the initial and final state interaction is derived: 

M = |l + i [,7(k)e^^^« - 1] • [1 + P/(p,k)]| Mo |l + ^ [77(p)e^^^(P) - l] • [1 + P.(p,k)]| , 

(41) 

where subscripts i and / indicate the initial and final state, respectively. The functions P(p, k) exhibit 
all the off-shell effects of the NN interaction and the primary production current [20] and S and 1] denote 
the phase-shift and inelasticity, correspondingly. At threshold, the inelasticity in the exit channel is 
equal to unity (^^(k) = 1) due to the small relative momentum of the outgoing nucleons. The last term 
of the formula expresses the influence of the initial state interaction on the production process. Due 
to the large relative momenta of the colliding protons needed to create a meson it is characterized by 
a weak energy dependence in the excess energy range of a few tens of MeV. For example in figure 11a 
one can see that the phase-shift variation of the ^Pq partial wave (having predominant contribution to 
the threshold production of pseudoscalar mesons) in the vicinity of threshold for mesons heavier than 
7r° is indeed very weak. Taking additionally into account that the initial state off-shell function Pi{p, k) 
is small (at least it is the case for meson exchange models [21]) one reduces the influence of the A^A^ 
initial state interaction to the reduction factor Fisi which can be estimated from the phase-shifts and 
inelasticities only: 

Fisi = \ \v{v) e'''^"^ + l\' . (42) 

At the threshold for tt meson production this is close to unity since at this energy the inelasticity is 
still nearly 1 and the ^Pq phase-shift is close to zero (sec figure 11a). However, at the r] threshold, where 
the phase-shift approaches its minimum, the proton-proton ISl diminishes the total cross section already 
by a factor of 0.2 [21]. A similar result was obtained using a meson exchange model for 77 production 
in the pp — > ppr) reaction and calculating the proton-proton distortion from the coupled-channel ttNN 
model [19]. The authors of reference [19] concluded that the initial proton-proton distortion reduces 
the total cross section by about a factor of ~ 0.26, which keeps constant at least in the studied range 
of 100 MeV in kinetic beam energy. Hence, the closed formula (42) disregarding the off-shell effects 
(Pi(p, k)) permits to estimate the cross section reduction due to the initial state distortion with an 
accuracy of about 25 %. The shape of the ^Pq phase-shift shown in figure 11a indicates that the effect 
is at most pronounced close to the rj production threshold, yet for the rj' meson the formula (42) leads 
to a factor Fjsi = 0.33 [133]. The primary production amplitude as well as the off-shell effects of 
the nucleon-nucleon FSI (P/(p,k)) are also weakly energy dependent [20], since they account for the 
short range creation mechanism, which shall be considered in the next section. The accordance of the 
experimental data with the simple factorization represented by solid lines in figures 13b, 14 and 15a 
fully confirms the above considerations which imply that the total cross section energy dependence 
and the occupation density on a Dalitz plot is in the first order determined by the on-shell scattering 
of the outgoing particles. However, since the distortion caused by the nucleons is by some orders of 
magnitude larger than that due to the meson-nucleon interaction, even small fractional inaccuracies in 
the description of nucleon-nucleon effects may obscure the inference on the meson-nucleon interaction. 

The differences between the square of the on-shell proton-proton scattering amplitude and the Jost 
function prescription are presented in figure 13a. To minimize the ambiguities that may result from 



35 



these discrepancies at least for the quantitative estimation of the effects of the unknown meson-nucleon 
interaction one can compare the spectra from the production of one meson to the spectra determined 
for the production of a meson whose interaction with nucleons is estabhshed. For instance, to estimate 
the strength of the rjpp and rj'pp FSI one can compare the appropriate observables to those of the 7r°pp 
system. 

Figures 18a and 18b show the dependence of \Mq\ on the phase-space volume for r] and rj' production 
normalized to \Mq \. The values of |Mo| were extracted from the experimental data by means of 
equation (38) disregarding the proton-meson interaction (|Mp^(^/)^p^(^/)| was set to 1). If the influence 
of the neglected interactions were the same in case of i]{ri') and 7r° the points would be consistent with 
the solid line. This is the case for the pp pprj' reaction visualizing the weakness of the proton-77' 
interaction independently of the prescription used for the proton-proton FSI [37]. In case of the r]' 
meson its low-energy interaction with the nucleons was expected to be very weak since there exists no 
baryonic resonance which would decay into Nrj' [134]. Figure 18a shows — independently of the model 
used for the correction of the proton-proton FSI — the strong effects of the r]pp FSI at low Vpg. 




Figure 18: The ratios of a) \M'^\/\Mf\ and b) \M'^\/\Mf\ extracted from the data by means of 
equation (38), assuming a pp-FSI enhancement factor as depicted by the dotted line in figure 13a and 
neglecting the proton-meson interaction [37]. c) Arbitrarily normalized enhancement factors for pd-, 
dfj- and pr^-FSl. The T^-proton factor is calculated according to equation (37), with apr, = 0.717 fm + 
i 0.263 fm [135] and bp^ — — l.SOfm — i0.24fm [125]. The enhancement factor for 77-deuteron has been 
extracted from the data of panel d) parametrizing the ratio of the cross section to the phase-space 
volume by the expression [136]: Fdr, = 1 + 0.5/(0.5 + (Q/5)^). The proton-deuteron FSI factor is 
calculated according to ref. [70]. The figure is taken from [136]. d) Total cross section of the quasi-free 
pn df] reaction as a function of the excess energy [48, 50] . The curve indicates the energy dependence 
proportional to the function ■ (1 + Q/83.5). The relative s- and p-wave contribution was taken as 
determined for the np — > d7r° reaction [43,44,47]. It was assumed that the p-wave to s-wave ratio is 
the same for the pn — > drj and pn — > dir^ reactions at a corresponding value of rjM- 
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With the up-to-date experimental accuracy, from all Meson NNsystems the rjNN one reveals by 
far the most interesting features. The dynamics of the rjNN system has become recently a subject 
of theoretical investigations in view of the possible existence of quasi-bound or resonant states [137]. 
According to reference [138] — within the present inaccuracy of a^Tv — the existence of quasi-bound 
?7-mesic light nuclei could be possible. A direct measure of the formation — or non-formation — of an 
77-nuclear quasi-bound state is the real part of the 77-nucleon scattering length [139]. The determined 
values of Re(a^iv) range between 0.20 fm [27] and 1.05 fm [144] depending on the analysis method and 
the reaction studied [122], and at present an univocal answer whether the attractive interaction between 
the 77 meson and nucleons is strong enough to form a quasi-bound state is not possible. 



The shape of the energy dependence of the pd — > ^He r) cross section implies that either the real or 
imaginary part of the 1] ^He scattering length has to be very large [140], which may be associated with 
a bound 7] ^He system. Similarly encouraging are results of reference [141], where it is argued that a 
three-body tjNN resonant state, which may be formed close to the rjd threshold, may evolve into a 
quasi-bound state for Re (a^Ar) > 0.733 fm. Also the close-to-threshold enhancement of the total cross 
section of the pp ppr) reaction was interpreted as being either a Borromean (quasi-bound) or a reso- 
nant r]pp state [142] provided that Re (flr/iv) > 0.7 fm. Contrary, recent calculations performed within a 
three-body formalism [137] indicate that a formation of a three-body rjNN resonance state is not pos- 
sible, independently of the rjN scattering parameters. Moreover, the authors of reference [143] exclude 
the possibility of the existence of an rjNN quasi-bound state. Results of both calculations [137, 143], 
although performed within a three-body formalism, have been based on the assumption of a separa- 
bility of the two-body rjN and A^A^ interactions. However, in the three-body system characterized by 
the pairwise attractive interactions, the particles can be pulled together so that their two-body poten- 
tials overlap, which may cause the appearance of qualitatively new features in the r^A^A^ system [137]. 
Specifically interesting is the rjd final state where the pair of nucleons alone is bound by the strong 
interaction. Figure 18d shows the total cross section for the pn — > drj reaction measured close to the 
production threshold. For excess energies below lOMeV the data are enhanced over the energy de- 
pendence determined for the pn — > dTi^ reaction indicated by the solid curve. This is in qualitative 
agreement with the calculations of Ueda [145] for the three-body rjNN-TrNN coupled system, which 
predict the existence of an 77A"A" quasi-bound state with a width of 20 MeV. Ueda pointed out that the 
binding of the r]NN system is due to the Sn r]N and ^Si NN(d) interaction which is characterized by no 
centrifugal repulsion. Such repulsion makes the ttNN system, in spite of the strong P33 nN attraction, 
hard to be bound [146]. Whether the observed cusp at the pn drj threshold is large enough to confirm 
the existence of the rjNN bound state is recently vigorously discussed [147, 148] . The considerations 
concerning the possible existence of the light 77-mesic nuclei will be continued in section 6.2. 



The enhancement factor for the deuteron-?] interaction inferred from the data in figure 18d varies 
much stronger in comparison to the proton-?; one, as demonstrated in figure 18c. This suggests that the 
effects of this interaction should be even more pronounced in the differential distributions of the cross 
section for the pd — > pdr] reaction as those observed in case of the pp — > ppr] process, especially because 
the "screening" from the proton-deuteron interaction is by more than an order of magnitude smaller 
compared to the proton-proton interaction as can be deduced from the comparison of the solid lines 
in figures 18c and 13a. The experimental investigations on that issue [136, 149] as well as the search 
for the ?7-mesic nuclei [150] by measuring proton-deuteron induced reactions in the vicinity of the 77 
production threshold are on the way. 
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4 Single meson production in NN scattering 



4.1 Dynamics of pseudoscalar meson production 

The considerations from the previous section led to the conclusion, that close to the kinematical thresh- 
old, the energy dependence of the total cross section is in the first approximation determined via the 

interaction among the outgoing particles and that the entire production dynamics manifests itself in a 
single constant, which determines the absolute scale of the total cross section. As a first step towards 
the understanding of the creation mechanism underlying the production let us compare the total cross 
sections of neutral pseudoscalar mesons tt^, rj and rj'. Since the masses of these mesons are significantly 
different the influence of the kinematical flux factor F on the total cross section and the suppression 
due to the initial state interaction Fjsi depend substantially on the created meson. Therefore, for the 
comparison of the primary dynamics we will correct for these factors and instead of comparing the total 
cross section we will introduce — according to reference [37] — a dimensionless quantity a ■ F/Fjsi, 
which depends only on the primary production amplitude Mq and on the flnal state interaction among 
the produced particles. 
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Figure 19: (a) Total cross section {a) multiplied by the flux factor F and divided by the initial state 
interaction reduction factor Fjsi versus the available phase-space volume for the reactions pp ppr] 
(squares [50,61-66]), pp ppm^ (diamonds [8,81,89,109,151]), and pp pprj' (circles [62,90,97,111, 
112]). (b) Experimentally determined ratio between pn — > pnr) and pp —>■ ppr] cross sections as a 
function of excess energy Q. The flgure is adapted from reference [49]. 

Close-to-threshold the initial state interaction, which reduces the total cross section, is dominated 
by proton-proton scattering in the ^Pq state which may be estimated in terms of phase-shifts and 
inelasticities by employing equation (42). The Fjsi factor is close to unity for pion production and 
amounts to ~ 0.2 [21] and ~ 0.33 [133] for the i] and rj' meson, respectively, at threshold. 

A comparative study of the production of mesons with signiflcantly different masses encounters 
the difficulty of flnding a proper variable at which the observed yield can be compared. The total 

^^The 7r°, r] and rj' masses amount to 134.98 MeV/c^, 547.30 MeV/c^ and 957.78 MeV/c^, respectively [134]. 
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cross section is defined as the integral over the available phase-space volume of transition probabilities 
— refiecting the dynamics of the process — from the initial to the final state as written explicitly in 
equation 5. Thus, if the dynamics of the production process of two different mesons were exactly the 
same then the above introduced yield would also be strictly the same for both mesons provided it was 
extracted at the same value of Vps- This inference would however not be valid if the production yields 
were compared at the tjm or Q variables. Therefore the volume of the available phase-space for the 
produced particles is the most suited quantity for the regarded comparison [37]. This could be also 
the best choice for the investigation of isospin breaking where the cross sections for the production of 
particles with different masses need to be compared (e.g. n^d pprj and n^d nnr] [152]). 

Figure 19a shows the yield of 7r°, r] and r]' mesons in the proton-proton interaction as a function 
of the available phase-space volume. The onset of higher partial waves is seen for 7r° and r] mesons in 
the Vps range between 10^ and 10^ MeV^, whereas the whole range covered by the r]' data seems to be 
consistent with the pure Ss production. One can also recognize that the data for the Ss final state are 
grouped on parallel lines indicating a dependence according to the power law a -F/Fisi ~ « ■ Vps^^ [1^3] 
and that over the relevant range of Vps the dynamics for r]' meson production is about six times weaker 
than for the tt^ meson, which again is a further factor of six weaker than that of the rj meson. This 
is an interesting observation, since the quark wave functions of rj and rj' comprise a similar amount of 
strangeness (~ 70% [154]) and hence, in the nucleon-nucleon collision one would expect both these 
mesons to be produced much less copiously than the meson tt^ being predominantly built out of up 
and down quarks. On the hadronic level, however, one can qualitatively argue that the 1] meson owes 
its rich creation in the threshold nucleon-nucleon collisions to the existence of the baryonic resonance 
A/"* (1535) whose branching ratio into the Nrj system amounts to 30-55% [134]. There is no such 
established resonance, which may decay into an s-wave rj'N system [134] and the 7r° meson production 
with the formation of the intermediate A(1232) state is strongly suppressed close-to-threshold, because 
of conservation laws. 

Prior to the more comprehensive discussion of the production mechanism for different mesons let 
us present one more interesting observation visualizing that the production dynamics depends strongly 
not only on the structure of the created meson but also on the isospin configuration of the colliding 
nucleons. 

For example, in the case of the t] meson the comparison of the cross section from the production 
in proton-proton and proton-neutron collisions revealed that the yield of the r] meson in the isosinglet 
nucleon configuration exceeds the one in the isospin triplet state by more than an order of magnitude. 
The ratio 

^ a{pn pnrj) 
^ cr{pp pprj) ' 

presented in figure 19b, was determined to be about 6.5 in the excess energy range between 16MeV 
and 109 MeV [49]. This implies that the production of the 7] meson with the total isospin equal to zero 
exceeds the production with isospin / = 1 by a factor of 12, since a{pn — > pnrj) = (cr/=o + (Ji=i)/2 
and a{pp — > pprj) — ai=i and hence ai=o — {2 Rjj — 1) (Ti=\. On the mesonic level, this large difference 
of the total cross section between the isospin channels suggests the dominance of isovector meson (tt 
and p) exchanges (figure 21 graph b) in the creation of t] in nucleon-nucleon collisions [3,49], but 
alternatively it can also be explained on the quark-gluon level assuming the instanton induced flavour- 
dependent quark-gluon interaction [24] (see below flgure 26 graph e). The above examples visualize 
one of the most interesting problems in the investigations of the dynamics of the close-to-threshold 
meson production: Namely the determination of the relevant degrees of freedom for the description of 
the nucleon-nucleon interaction, especially in case when the nucleons are very close together As 

^^These investigations are listed as one of the key issues in hadronic physics [155]. 
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pointed out by Nakayama [20] , the transition region from the hadronic to constituent quark degrees of 
freedom does not have a well defined boundary and at present both approaches should be evaluated in 
order to test their relevance in the description of close-to-threshold meson production in the collision 
of nucleons. 
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Figure 20: A cartoon illustrating in naive 
geometrical terms that forr < 2Y]\f + 2tm 
meson exchange is unlikely to be appro- 
priate to the description of the internu- 
cleon potential. Figure and caption are 
taken from reference [156]. 



Simple geometrical considerations presented in figure 20 indicate that at distances smaller than 2 fm 
the internucleon potential should begin to be free of meson exchange effects and may be dominated 
by the residual colour forces [156]. In the previous section it was shown, that the close-to-threshold 
production of mesons occurs when the colliding nucleons approach distances of about 0.5 fm in case of 
ir^ and of about 0.18 fm in case of production (see table 2). This distance is by about an order of 
magnitude smaller than 2 fm and it is rather difficult to imagine — in coordinate space — an exchange 
of mesons between nucleons as mechanism of the creation process. Such small collision parameters 
imply that the interacting nucleons - objects of about 1 fm - overlap and their internal degrees of 
freedom may be of importance. 
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Figure 21: Diagrams for the pp — » pp Meson reaction near threshold: (a) - 
(nucleonic current) (b) — "rescattering" term (nucleonic current) (c) 
meson-exchange (d) — emission from virtual meson (mesonic current) 
intermediate resonance (nucleon resonance current). 
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Nonetheless, in the last decade, the effective theory based on the meson exchanges, which accounts 
for the size of the participating particles by introduction of the momentum transfer dependent form 
factors, has been extensively utilized for the description of the creation process. Figure 21 represents 
the regarded mechanisms. 

The question which processes — on the mesonic level — are responsible for the production was 
considered already in 1966 by Koltun and Reitan [6]. The authors deemed the direct pion production on 
one of the protons (fig. 21a) to play the most important role in the pp — > ppir^ reaction. This supposition 
was anticipated until the precise experiments performed at lUCF [8,9] and CELSIUS [89] revealed that 
this mechanism can not account for more than 5 % of the total yield, as demonstrated in figure 22. 
Subsequent investigations showed that the inclusion of the processes illustrated in figure 21b, when the 
pion is produced on one of the protons and then scatters on the other one (in the irN s-wave [157] 
and nN p-wave via A(1232) [159]), increases the cross section essentially yet not enough to describe 
the data, even if the off-shell properties of the ttN amplitude are taken into account [12,160]. The 
contribution of the rescattering mechanism is shown as the dashed line in figure 22. The agreement 
with the data can, however, be achieved — as shown by the solid line — if one takes additionally into 
account the mechanisms of heavy-meson-exchange associated with an intermediate virtual nucleon- 
antinucleon state [11,161] illustrated in figure 21c. A thorough discussion concerning tt meson production 
in differently charged channels can be found e.g. in references [13, 162]. 



41 




Figure 22: Total cross section for the reaction 
PP ~^ ppn°. The dotted hne shows the results 
for direct production only. Adding rescattering 
gives the dashed line. Including also the contri- 
bution from heavy meson exchange leads to the 
solid line. The data are from refs. [8,89]. The 
figure is adapted from reference [158]. Please 
note, that similarly as shown in figure 13b, 
the solid line representing the full calculation 
starts to deviate from the data at r]M ~ 0.4 
(Q pa llMeV), yet here the magnitude reflects 
the assumed mechanism and in figure 13 it was 
settled to fit the data. Thus indeed, the primary 
dynamics fixes in practice only the scale. 
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The dominance of the heavy-meson-exchange and the off-shell nN rescattering mechanisms, both 
refiecting the short range nucleon-nucleon interactions, confirms the simple estimation of the reaction 
region based on the momentum transfer between the colliding nucleons. Thus, in case of the 7r° meson 
creation, this approach succeeded to describe the magnitude of the total cross section. However, the 
meson-exchange model at the present stage does not reproduce the angular distributions of the spin 
observables for all pp NNn isospin channels, simultaneously. Specifically, though the only available 
calculations of the spin correlation coefficients based on the Jiilich model [14] are in a good agreement 
with the experimental data for the reactions pp pnn~^ [168] and pp — > dir'^ [169] they reveal se- 
rious discrepancies when confronted with the angular dependences of the spin correlation coefficients 
Ay, = + Ayy and A^z of the pp ppTi^ reaction [82]. Therefore, further investigations concnerning 
the relative contributions of various mechanisms, as well as the different partial waves are still required. 

The role of the mechanisms leading to the tt production have been also considered in context of the 
lowest-order chiral perturbation theory [34] and relativistic impulse approximation [163] . In contrast to 
the meson exchange models the chiral perturbation method results [164-166] in the destructive inter- 
ference between the direct and the pion-rescattering terms, and consequently implies more significant 
contributions originating from the shorter-than-pion-range mechanisms, like the heavy-meson exchanges 
or the creation of pions in the fusion of exchanged mesons [167]. 

In case of the rj and rj' mesons the investigations of the mechanisms underlying the production 
process are even more difficult. It is partly due to the fact that in contrary to the pion case even the 
coupling constants determining the strength of the NNrj^r]') vertex are poorly established. 

The values derived by various methods differ by more than an order of magnitude [154, 170]. For ex- 
ample, the estimation based on the dispersion relation results in gNNrj' smaller than 1 [171], whereas the 
value obtained from the fits to low-energy nucleon-nucleon scattering in one-boson-exchange models 
amounts to 7.3 [172]. This uncertainty prevents from the exact estimation of the contribution from the 
nucleonic current. However, according to the systematic investigations of Nakayama et al. [20,133,173] 
this mechanism is of minor importance. The calculations have been performed with gNN-q' = QNNri = 6.1 
— used in the NN scattering analysis [174] — which is close to the upper limit of the predicted values. 

It is at present rather well estabhshed [19, 33, 121, 173, 175-179] that the 77 meson is produced 
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predominantly via the excitation of the Sn baryonic resonance A^*(1535) which subsequently decays 
into 7] and nucleon and whose creation is induced through the exchange of vr, 77, p, cr and u mesons, 
as shown in figure 21e. Though all the quoted groups reproduce the magnitude of the total cross 
section, their models differ significantly as far as the relative contribution from the vr, rj and p exchange 
mechanisms is concerned. The discrepancies are due to the not well known strength of Meson-NSn 
couplings and the rjN scattering potential. For example, while the dominance of the p meson exchange 
is anticipated by authors of references [33, 121, 175-177], it is rather the exchange of the rj meson 
which dominates the production if one takes into account effects of the off-shell rjN scattering [180], 
or utilizes the multi-channel multi-resonance model [19]. In any case, the hitherto performed studies 
aiming to describe the total cross section show that the close-to-threshold production of r] mesons 
in nucleon-nucleon collisons is dominated by the intermediate virtual Sn nucleon isobar whose width 
overlaps with the threshold. In order to disentangle between various scenarios of the Sn excitation a 
confrontation of the predictions with other observables is needed. The interference between considered 
amplitudes causes a different behaviour — depending on the assumed scenario — e.g. of the rj meson 
angular distributions. These differences, however, are too weak in the close-to-threshold region to judge 
between different models. Also the ratio of the 1] meson production via the reactions pp — >■ pprj and 
pn — > pnr] can be equally well described either by assuming the p meson exchange dominance [177] or by 
taking pseudoscalar and vector mesons for exciting the Sn resonance [173]. In the latter case, shown in 
figure 23, the excitation of the resonance via the p meson exchange was found to be negligible. Yet, the 
predictions of the analyzing power depends crucially on the assumed mechanism [173,177]. This fact 
has triggered already experimental investigations which aim to determine the spin observables [181]. 

1 r.2 



Figure 23: Total cross sections for 
the pp pprj (upper panel) and 
pn — > pnrj (lower panel) reactions 
as a function of excess energy. The 
dashed curves correspond to the nu- 
cleonic current contribution and the 
dash-dotted curves to the mesonic 
current consisting of r/pp, r]ujuj and 
r/aovr contributions. The resonance 
current presented by the dotted line 
consists of the predominant Sii(1535) 
and of the Pii(1440) and Di3(1520) 
resonances excited via exchange of tt, 
7], p and uj mesons. The solid curves 
are the total contribution. The devi- 
ation from the data at low Q in the 
upper panel reflects the rjp FSI which 
was not included in the calculations. 
The data are from refs. [49, 61-65]. 
20 40 60 80 100 The figure is taken from [20]. 

Q (MeV) 

Figure 23 presents also a contribution from the mesonic current where the rj meson is created in 
the fusion of virtual e.g. p ot u mesons emitted from both colliding nucleons. In comparison to the 
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overwhelming strength of the resonance current this mechanism is however by a factor of 30 too weak. 
In contrast, it is enough to explain the magnitude of the close-to-threshold rl meson production in the 
proton-proton interaction, which is just by about a factor of 30 smaller compared to the 77 meson. Among 
the mesonic currents regarded by authors of reference [133] the pprj gives the dominant contribution, 
by a factor of five larger than the arjr]'- and a;a;?7'-exchange. However, the understanding of the r}' 
production on the hadronic level is far from being satisfactory. The magnitude of the total cross section 
was also well reproduced in the frame of a one-boson-exchange model (nucleonic current) where the 
virtual boson {B = it, 1], a, p, uj, ag) created on one of the colliding protons converts to the rj' on the 
other one [35]. Taking into account the off-shell effects of the Bp — > r/'p amplitude it was found that 
the short range a and p meson exchanges dominate the creation process, whereas the tt exchange plays 
a minor role. On the contrary, other authors [182] reproduced the magnitude of the total cross section 
with TT exchange only and concluded that the rj, p and cu exchange currents either play no role or cancel 
each other. However, in both of the above quoted calculations [35, 182] the initial state interaction 
between protons, which reduces the rate by a factor of about 3, was not taken into account and hence 
the obtained results would in any case reproduce only 30 % of the entire magnitude and could be at 
least qualitatively reconciled with the mentioned result of reference [133], where the nucleonic current 
was found to be small. Moreover, the choice of the another prescription for the form factors could 
reduce the one-pion-exchange contribution substantially [133]. However, the picture that the rj' meson 
is predominantly created through the mesonic current, which we would like to advocate, remains at 
present unclear as well. This is because the magnitude of the total cross section could have been also 
described assuming that the production of rj' is resonant [133] . As possible intermediary resonances the 
recently reported [183] Sii(1897) and Pii(1986) have been considered. These resonances were deduced 
from rj' photoproduction data, under the assumption that the close-to-thrcshold enhancement observed 
for the 7p — > rj'p reaction can be utterly assigned to resonance production. Further, as well strong, 
assumptions have been made in the derivation of the gNN*ri' and gNN*ir coupling constants [133]. Hence, 
it is rather fair to state that in the case of the close-to-threshold tt^ and rj meson production in nucleon- 
nucleon coUisons the dynamics is roughly understood on the hadronic level, but the mechanisms leading 
to the rj' creation are still relatively unknown. Unil now it was not possible to estimate satisfactorily 
the relative contributions of the nucleonic, mesonic and resonance current to the production process. 
In fact, model uncertainties allow that each one separately could describe the absolute values of the 
pp — > ppr]' total cross section. This rather pessimistic conclusion calls for further experimental and 
theoretical research. The understanding of the production dynamics of the rj' meson on both hadronic 
and the quark-gluon level is particularly important since its wave function comprises a significant 
gluonic component [184], distinguishing it from other mesons and hence the comprehension of the 
mechanism leading to its creation in collisions of hadrons may help in the determination of its quark- 
gluon structure which is also of importance to investigate possible glueball candidates [184]. Hereafter 
we will briefly report about the gluonic mechanisms which may — additionally to the meson exchange 
processes discussed above — contribute to the 77' production. 

As already mentioned, the close-to-threshold production of 1] and rj' mesons in the nucleon-nucleon 
interaction requires a large momentum transfer between the miclcons and hence can occur only at dis- 
tances of about 0.3 fm (see table 2). This suggests that the quark-gluon degrees of freedom may indeed 
play a significant role in the production dynamics of these mesons. A possibly large glue content of the 
Tj' and the dominant flavour-singlet combination of its quark wave function may cause that the dynam- 
ics of its production process in nucleon-nucleon collisions is significantly different from that responsible 
for the production of other mesons. In particular, the rj' meson can be efficiently created via a "contact 
interaction" from the glue which is excited in the interaction region of the colliding nucleons [185]. A 
gluon-induced contact interaction contributing to the close-to-threshold pp — > pprj' reaction derived in 
the frame of the U(l)-anomaly extended chiral Lagrangian is discussed in references [185-187]. The 
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strength of this contact term is related to the amount of spin carried by polarized gluons in a polar- 
ized proton [187, 188] , thus making the study of the close-to-threshold 77' meson production even more 
interesting. 

Figure 24 depicts possible short-range mechanisms which may lead to the creation of the vl meson 
via a fusion of gluons emitted from the exchanged quarks of the colliding protons [189] or via an exchange 
of a colour-singlet object made up from glue, which then re-scatters and converts into 77' [191]. The 
hadronization of gluons to the v( meson may proceed directly via its gluonic component or through its 
overwhelming flavour-singlet admixture ri\ (see fig. 25). Contrary to the significant meson exchange 
mechanisms and the fusion of gluons of figure 24 graph a), the creation through the colour-singlet object 
proposed by S.D. Bass (graph 24b) is isospin independent, and hence should lead to the same production 
yield of the rj meson in both reactions (pp pprj' and pn pnr]') because gluons do not distinguish 
between flavours. This property should allow to test the relevance of a short range gluonic term [192] 
by the experimental determination of the cross section ratio Rr,' = (^{pn pnri')/a{pp — > pprj'), which 
in that case should be close to unity after correcting for the final and initial state interaction between 
participating baryons. The other extreme scenario — assuming the dominance of the isovector meson 
exchange mechanism — should result in the value of R^' close to 6.5 as was already established in the 
case of the rj meson (see fig. 19b). The experimental investigations of that issue are in preparation [51]. 




Figure 24: Diagrams depicting 
possible quark-gluon dynam- 
ics of the reaction pp pprj'. 
(a) — production via a fusion 
of gluons [189] with rearrange- 
ment of quarks, (b) — produc- 
tion via a rescattering of a "low 
energy pomeron" . 
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Figure 25: Coupling of rj and rj' to two gluons through (a) quark and antiquark triangle loop and 
(b) gluonic admixture. rji denotes the fiavour-singlet quark-antiquark state. This indicates that 
gluons may convert into the rj or rj' meson via a triangle quark loop only by coupling through their 
fiavour singlet part. The figure is taken from reference [190]. 
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In addition to the QCD motivated investigation presented above [185,186], the interesting features of 
the close-to-threshold meson production, in particular the large cross section for the rj meson in proton- 
proton interactions exceeding the one of the pion, or even more surprisingly large cross section of the 
7] production in proton-neutron collisions, encouraged also other authors to seek for the underlying — 
OZI rule violating — creation mechanisms in the frame of microscopic models of QCD [24,25,193]. 
The hitherto regarded processes are presented in figure 26. As indicated in the pictures the structure 
of participating baryons has been modeled as quark-diquark objects with harmonic confinement [194]. 
In the upper graphs the large momentum transfer is shared by the exchanged gluons with a subsequent 
interchange of quarks to provide a colourless object in the final state. The lower graphs depict the two 
examples of instanton induced interactions with a 6-quark-antiquark (fig. 26d) and two-gluon vertex 
(fig. 26e). 



Adjusting the normalization to the cross section of the pp ppir'^ reaction at a single energy point 
the model [25] accounts roughly for close-to-threshold cross sections of other pseudoscalar and vector 
mesons in proton-proton collisions. Though this approach is characterized by significantly smaller 
number of parameters than in the meson exchange models, their uncertainties allow for the description 
of the data equally well either by the gluon exchange or by the instanton induced interactions, at least 
in case of the vr*^, 77' , uj and mesons. Yet for the resonance dominated 7] and meson production in 
proton-proton collisions it was found that the instanton induced interaction presented by graph 26d is 
not sufficient. Similarly, the authors of reference [24] argue that the instanton induced interaction with 
a quark-gluon vertex (graph 26e) should be of no importance for the r] production in proton-proton 
collisions. This arises from the properties of the vertex which lead to the r] production only in case 
of an interaction between quarks of different flavours {ud udgg) and in the quark-diquark model of 
baryons the proton consists of -urf-diquark and -u-quark and correspondingly the neutron is modeled 
as M(i-diquark and d-quark. Although negligible in case of proton-proton interactions, this mechanism 
may contribute significantly to the total cross section of the pn — > pnr] reaction and indeed as shown in 
reference [24] it reproduces the data. 
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Figure 26: Diagrams for the qq — > qq{qq) production operator: (a), (b) Two-gluon exchange and 
rescattering mechanism (c) Correlated, colourless two-gluon exchange. The dashed line indicates the 
excitation of an intermediate NN* system. (d) Instanton induced 6-quark interaction. Figures (a-d) 
according to [25]. (e) The instanton contribution to the rj meson production in the proton-neutron 
interaction. Figure copied from [24]. 

However, the magnitude of the cross section is very sensitive to the size of the instanton in the QCD 
vacuum and — similar to the uncertainty of coupling constants in case of the hadronic approach — 
at present it makes precise predictions impossible. Thus, despite the partial successes, at the present 
stage of developments both approaches — on the quark-gluon and on the hadronic level — do not 
deliver an unambiguous answer to the dynamics of the close-to-threshold meson production in the 
nucleon-nucleon interaction. 



4.2 Vector meson production 

Vector meson production studies in the mass range up to 1 GeV/c^ in hadronic collisions are especially 
motivated by questions concerning the strangeness content of the nucleon. Within the quark model 
the (f) meson is an (almost) ideally mixed ss state, whereas the nucleon is composed only of u- and 
(i-quarks. According to the OZl rule [195] (p production in the nucleon-nucleon interaction would 
be forbidden as a process with disconnected quark lines, if the (p was a pure ss state. 

However, due to a small deviation from ideal mixing (Sy — 3.7° [134]) the (f) meson is allowed to 
couple to the nucleon through its (uu + dd) admixture. Quantitatively, the naive OZI rule results in a 

^^for a review see [196]. 
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suppression of (j) compared to 00 production in hadronic interactions 



after phase-space corrections, where A, B and X denote hadronic systems consisting only of hght 
quarks [197]. Any significant deviation might be interpreted as a hint for a strangeness component in 
the nucleon. A considerable contribution of strange sea quarks to the nucleon wave function is suggested 
by results on the E^jv term in pion nucleon elastic scattering [198], on the nucleon's structure function 
in deep inelastic scattering with polarized muons [199] and, more recently, on charm production in deep 
inelastic neutrino scattering [200,201]. 

Large violations of the OZI expectation (eq. (43)) have been found in pp annihilation experiments by 
the ASTERIX, Crystal Barrel and OBELIX collaborations at LEAR (as reviewed in [202]). However, 
deviations are predominantly found in pp annihilation at rest as compared to higher energies and seem 

to be restricted to S-wave processes only. Furthermore, the results are strongly dependent on the final 
state: While the (j) tt and 07 channels exceed the OZI value by a factor of 20 and 100, respectively, little 
or no effect is seen in the tt tt and 77 final states. 

These data have been interpreted as "shake-out" and "re- arrangement" of a negatively polarized 
ss Fock space component of the proton wave function [202,203], which also accounts for large double 
00 cross sections [204] in pp annihilation measured at the JETSET experiment [205]. However, OZI- 
allowed two-step processes via intermediate KK or AA states have been shown to describe available 
data from pp annihilation without any strangeness in the nucleon [206-213]. 

In nucleon-nucleon scattering, effects of competing two step processes for production via inter- 
mediate KK or UK states are expected to be of minor importance [212,214]. Furthermore, in close- 
to-threshold production the entrance channel must be in a ^Pi state (see table 1) due to parity and 
angular momentum conservation and it is the spin triplet fraction in pp annihilation which is strongly 
correlated to the meson yield [215]. Thus, cross section ratios for and u production in proton-proton 
scattering should clearly indicate possible OZI violations and probe the ss component of the nucleon. 

The exclusive production ratio has been determined in the reaction pp pp(l)/uj at a beam momen- 
tum of 3.67 GeV/c corresponding to excess energies of 82 MeV and 319 MeV to the and uj thresholds, 
respectively, by the DISTO collaboration at SATURNE [216]. After phase-space corrections, the ob- 
served ratio exceeds the OZI expectation by an order of magnitude, while data at excess energies larger 
than 1.1 GeV [217] show an enhancement smaller by a factor of four (see fig. 27). One-boson-cxchange 
models [218-221] underestimate the absolute value of the production cross section measured at the 
DISTO experiment by at least a factor of three. 
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Figure 27: total cross section ra- 
tio in proton-proton collisions versus 
excess energy with respect to the 
production threshold from [216]. The 
recent data point from the DISTO 
collaboration (sohd square [216]) is 
shown together with high energy data 
(solid circles [217]) and calculations 
taking into account the OZI prediction 
and three-body phase-space [216] and 
results obtained within a one-boson- 
exchange model including the proton- 
proton FSl [218]. Data includes both 
statistical and systematical errors. 



Figure 27 includes results obtained within a one-pion-exchange model [218], which provides a rea- 
sonable description of high energy data. Model predictions — taking into account the proton-proton 
final state interaction only — falls short of the recent close-to-threshold value [216]. Using instead 
a parametrization of the energy dependence of the uj production cross section and considering both 
the finite uj width as well as effects of proton-proton FSl [219,222], the (^juj ratio can be determined 
at equal excess energies, thus reducing uncertainties due to the available phase-space, partial wave 
amphtudes and effects of the proton-proton FSL In conclusion, the excess over the OZI prediction 
is decreased to a factor of five [216], i.e. to a modest enhancement as compared to the results from pp 
annihilation experiments 

In the DISTO results, production at an excess energy of 82 MeV appears to be dominated by 
s-wave relative to the nucleon-nucleon system in agreement with theoretical predictions [224] and 
phenomenologically described in section 3. However, the proton-proton angular distribution relative 
to the direction is found to deviate significantly from isotropic emission at this beam energy with 
evidence for a P-wave contribution in the proton-proton subsystem. Consequently, production at 
82 MeV above threshold, although dominated by a ^Pi proton-proton entrance channel, partly proceeds 
via the ^Sq and ^V)i spin singlet partial waves [216]. Thus, the influence of the spin-triplet entrance 
channel, prerequisite to a "re-arrangement" of an ss component in the nucleon wave function [202], 
might be diluted in the presently available data and would be expected to be more prominent in data 
closer to the production threshold. 



^^The w angular distribution relative to the system measured at an excess energy of 319 MeV shows evidence for p- 
and d- waves contributing [216]. 

^^Comparing the DISTO data at Q = 82 MeV [216] and the recent COSY-TOF result at an excess energy of 92 MeV [223] 
leads to a similar enhancement of w 7 compared to the OZI expectation. 
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Figure 28: Vector meson production currents in nucleon-nucleon scattering [NN — > NNv) studied 
within a relativistic meson exchange model in [20,225,226]: 
(a) nucleonic current and (b) mesonic current. 

M — 77,7), p, cu, a, ao denotes the exchanged meson and B the intermediate off mass-shell baryon. 



Recently, within the framework of a relativistic meson exchange model, nucleonic and meson- 
exchange currents have been used explicitly to study vector meson production in nucleon-nucleon scat- 
tering (fig. 28 [20,225,226]). Uncertainties in the model calculation for the nucleonic current arise mainly 
from a large range of values quoted for vector and tensor NNv coupling constants [171,174,212,227-232], 
as well as from the vector and tensor vertex form factors for the coupling of an on-mass shell vector 
meson to off mass-shell nucleons. Note, that a resonance contribution to the nucleonic current (fig. 28a) 
is neglected in the model [225], since there are no isospin-1/2 nucleon resonances experimentally known 
coupling to the vp channel [233]. Due to the strength of both the A'^A'^tt and NNp couplings and the g-rrpw 
coupling constant as compared to other possible choices, only the Trpv vertex is considered in evaluating 
the meson-exchange current. For a; production, contributions from the latter and the nucleonic current 
are found to interfere destructively when added coherently, which leaves two choices for a given value of 
the NNlo tensor-to-vector coupling ratio — when fitting the close-to-threshold data on the reaction 
PP ~^ PP^ taken at SATURNE [236] — for the only remaining free parameter of the model, the cut-off 
parameter at the NNuo form factor. Effects from a destructive interference of mesonic exchange current 
and nucleonic current contribution are also concluded by other authors for uj production in [234] and 
for production in [235]. In Nakayama's analysis the energy dependence of the total cross section is 
slightly influenced by either choice, as shown by the solid and dashed curves in figure 29 together 
with a paramctrization based on onc-pion-exchange (dash-dotted line [218]). In comparison, total cross 
section data taken at the SPES 3 spectrometer up to excess energies of 30 MeV [236] , the DISTO data 
point discussed above [216], high energy data from [238] and recent results obtained at the TOF facility 
at COSY at excess energies of Q = 92 MeV and 173 MeV, respectively [223], are included. 



^''Kaiser's study [234] is based on the calculation of various tree level diagrams in comparison with the threshold 
transition amplitude extracted from data of [236] and follows the approach developed in [34,237] (see section 4.4). The 
dominant nucleonic current contribution arisc;s from exchange. Howcvcir. the combined result of the studied amplitudes 
shows similarly to [225, 226] a dependence on the choice of the tensor-to-vector coupling ratio. 

^^The two curves correspond to the extreme cases found for the tensor coupling values Jmnv = ±0.5 qnnv [226], with 
cut-off masses of 1545 MeV (dashed line) and 1170 MeV (solid line), respectively. 
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Figure 29: Total cross section 
for pp ppuj as a function of 
the excess energy. Data are from 
refs. [216] (open triangle), [223] 
(solid squares), [236] (solid cir- 
cles) and [238] (open squares). 
Calculations include a one-pion- 
exchange based parametrization 
(dash-dotted curve [218]) and re- 
sults obtained within a relativis- 
tic meson-exchange model con- 
sidering nucleonic and mesonic 
current contributions (solid and 
dashed curves [225,226], for de- 
tails see text). The figure is taken 
from ref. [223]. 

The one-pion-exchange parametrization [218] deviates strongly from the close-to-threshold data, 
where FSI effects become important, which have not been included in the calculation. Both predictions 
obtained within the relativistic meson exchange model of [20, 225, 226] overestimate the recent COSY- 
TOF data [223] by at least a factor of five. The solid curve in figure 29 seems to follow the data slightly 
better, which would correspond to a lower cut-off parameter in the nucleonic current and, consequently, 
a dominating meson exchange current contribution. 

Within the model [225,226], angular distributions of the produced mesons are expected to differ 
depending on the dominant production mechanism for both uj and production: While the mesonic 
current leads to an isotropic angular distribution, the nucleonic current results in a cos^ behaviour for 
the vector meson emission angle 0^, in the overall centre-of-mass system Thus, angular distributions 
offer a unique possibility to separate both contributions, or to isolate the nucleonic current contribution 
in order to determine the (pNN coupling directly in a combined analysis of cu and production data 
in pp induced reactions: The experimentally observed isotropic (p angular distribution [216] clearly 
favours a dominant (ppTr meson exchange over the nucleonic current [20,212]. Contrary, the centre-of- 
mass angular distribution of the uj meson obtained at the COSY-TOF facility at an excess energy of 
173 MeV [223] exhibits a strong anisotropy similar to the observation by the DISTO collaboration at 
an excess energy of 319 MeV [216], indicating the presence of higher partial waves and rather favours a 
dominant nucleonic current contribution instead of meson exchange currents suggested by the shape of 
the energy dependence of the total cross section in figure 29 

Further constraints on the relevant vector meson production mechanisms in nuclcon-nuclcon scat- 
tering might be obtained from comparing with complementary meson induced reactions: An analysis 
of ttN — > N(f)/u! data concludes cu production to result from strong interferences between the meson ex- 
change current and both nucleonic and resonance currents, which contribute with similar strengths [240]. 
Near threshold, the Pn A^*(1440) resonance and the Sn A^*(1535) and A/"* (1650) excitations are found 
to contribute, with the Sn contributions canceling to a large extent 

In conclusion, at least for the reaction nN — > Nuj/(f), the nucleonic current contribution appears to 

^^The cos^ 8^ behaviour arises according to [225] from the spin-dependent part of the nucleonic current contribution, 
the "magnetization current" [239]. 

^°It should be noted, that the excess energy of the COSY-TOF data [223] is slightly out of the energy range covered 
by the predictions from Nakayama et al. [225, 226] due to the onset of inelasticities in the final state. 

^^For the resonance contributions in the work of Titov et al. [240], u)NN* couplings and relative phases have been 
taken from [241]. 
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be much larger for cu production relative to the meson exchange current, contrary to 4> production [240]. 
Thus, the scenario from pion induced vector meson production seems to agree — qualitatively — with 

the observations in nucleon-nucleon scattering, and is interesting to note especially in view of the u 
angular distribution from the COSY-TOF data [223]. However, angular distributions at lower excess 
energies would clearly help to constrain the interpretation of the presently available nucleon-nucleon 
data in future. 

In a combined analysis of cu and production within the relativistic meson exchange framework [20] 

a vahic of the coupling to the nucleon has been extracted from the DISTO data and is found to be 
consistent with the OZl expectation gNN<i> = gNNpSinSv ~ — (0.60 ± 0.15) [226]. Yet, as stated by 
Nakayama, the value still has large uncertainties, especially due to the lack of more and precise data in 
the threshold region [20]. Indeed, on the basis of extracted transition matrix elements [219] combining 
ttN and pp induced reactions, 8k(t)/uj production ratio has been derived, which exceeds the OZI estimate 
of eq. (43) by a factor of five in both cases. 

As suggested in [216], future experiments on (j) and oo production in close-to-threshold proton- 
neutron collisions might provide further evidence for the intrinsic strangeness content of the nucleon, 
as predictions for production cross sections based on the "shake-out" and "rearrangement" mecha- 
nism [242] and on meson exchange models [221,224] differ by more than an order of magnitude. 

4.3 Scalar sector 

The 1 GeV / c^ meson mass range is continuously under discussion regarding the nature of the scalar reso- 
nances /o(980) and ao(980), which have been interpreted as exotic two-quark two-antiquark states [243], 
possible glueball candidates [244], excitations of the light quark vacuum condensate {minions [245]), 
conventional qq states [246,247], or molecular like KK bound states [248,249]. Just already the similar- 
ity of the masses of these objects with twice the kaon mass, see figure 30a, suggests a strong correlation 
between these resonances and the K-K system. 

Within the framework of the Jiihch meson exchange model for tttt [251] and ttt^ [252] scattering 
the KK interaction dominated by vector meson exchange gives rise to a KK bound state identified 
with the /o(980) in the isoscalar sector, while the isovector ao(980) is concluded to be a dynamically 
generated threshold effect [252,253]. 

As shown in figure 30b, both shape and absolute scale of tttt — > KK transitions turn out to depend 
crucially on the strength of the KK interaction, which, in turn is prerequisite of a KK bound state 
interpretation of the /o(980). Similar effects might be expected for the elementary kaon-antikaon 
production in proton-proton scattering. 

Although the KK decay mode of /o(980) and ao(980) is rather weak in comparison to the dominant 
TTTT and Tirj decay channels [134], even a new theoretical analysis based on the chiral approach can not 
account for the /o(980) and ao(980) if the KK channel is not introduced additionally to the tttt and 
T^rj interaction [254]. An analysis of the tttt and KK interaction [255] showed that the /o corresponds 
to two poles on the second and third Riemann sheet, respectively, and appears physically as an object 
with a decay width of about 400 MeV and a narrow peak width of about 50 MeV. The same parameters 
of the /o were found by utilizing a unitarized quark model, according to which the /q was interpreted as 
a qq state with a large admixture of a KK virtual state [256]. The origin of the scalar resonances was 
also thoroughly studied by means of a coupled channel analysis considering tttt, KK and a a meson- 
meson scattering [257,258]. Decreasing gradually the interchannel coupling constants it was inferred 
that for some solutions at the limit of the fully uncoupled case the /o corresponds to a KK bound 
state [258,259]. 

Recently notable progress has been made concerning the discussion of this topic due to approaches 
in unitary extensions of chiral perturbation theory where the resonances under consideration qualify as 
dynamically generated resonances from the multiple scattering with the lowest order chiral Lagrangian, 
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as outlined in the review [260] . 




Figure 30: (a) Weighted mass average for the /o(980) scalar resonance from [250]. (b) Predictions of 
the Jiilich model for tttt — > KK transitions using the full model and neglecting any KK interaction, 
taken from ref. [252]. 



In high energy experiments, the /o meson is observed as a resonance in the system of two pions 
produced in a variety of hadro-production reactions [261-264] and in the hadronic decays of either 
heavier mesons [265-267] or the Z° boson [268,269] created in e^e~ coUisions. These experiments study 
the invariant masses of the created neutral (vr'^vr^) [261,263,267] and charged {'k'^'k") [262,264-266,268] 
pion pairs. Similarly, charged [270,271] and neutral [272] mesons were observed as a clear signal in 
an invariant mass spectrum of the rjn system. 

Complementary to these approaches, studying the interaction of tttt, KK and 7r?7 meson pairs, the 
COSY-11 collaboration presently investigates the possible manifestation of the mesons /o and/or oq 
as doorway states leading to meson production in proton-proton collisions, namely pp pp/o(oo) 
pp Mesons. By measuring the missing mass of the pp-system the overlapping f 0-0,0 resonances are 
studied as a genuine particle produced directly at the reaction. The production of these resonances 
is carried out at a mass range a few tens of MeV below the KK threshold (where the resonances can 
only decay into non-strange mesons) and above the KK thresholds, where the equivalent part of these 
broad resonances can be excited. 

It is obvious that for the production of a broad resonance the phrase "close-to-threshold" is not well 
defined and implies here that the beam momentum is such that masses just in the range of the resonance 
can be excited. Different mass ranges of such objects arc produced with different excess energies making 
the interpretation of the experimental data difficult but still possible [278]. It is worth noting, that 
recently first measurements of the fo meson production relatively close to its threshold but still above 
the KK threshold were performed by observing the (f) meson decays into the /07 channel [266,267], where 
the (f) is only about 40MeV/c^ heavier than the fo- In a recent publication [77] data on the close-to- 
threshold K^K^ production following the proton-proton interaction at an excess energy of Q = 17 MeV 
are presented. The obtained distribution of the missing mass to the pp system is shown in figure 31 
and demonstrates that the non-resonant K^K~ production (solid line) is hardly distinguishable from 
the resonant pp — > pp/o(980) ppK'^K' reaction sequence (dashed line) [273]. 
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Figure 31: Experimental spectrum of the 
K^K~ invariant mass measured for the 
reaction pp — > ppK^K~ at a beam mo- 
mentum of 3.356 GeV/c corresponding to 
an excess energy of Q = 17MeV with 
respect to the K^K~ threshold (data 
points). The width of the bins cor- 
responds to the experimental resolution 
of the mass determination (FWHM ^ 
2 Me V) . The sohd and dashed hues show 
Monte-Carlo simulations assuming the 
direct and resonant production, corre- 
spondingly [273]. 
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Since the statistics of the data were not sufficient to favour one of the two processes, the cross section 
was extracted for both and resulted in nearly identical values of: (J non-resonant — (1-80 ± 0.27 ^oJf) nb 
and aresonant = (1-84 ± 0.29]'^q;33) nb including statistical and systematical errors, respectively. The is- 
sue, whether there is a chance to distinguish between KK pairs originating from the decay of a genuine 
/o/oo resonance and those produced by a strong tttt-KK correlation is at present under theoretical 
investigation [253]. Still, calculations based upon the one-pion exchange and a Breit-Wigner presenta- 
tion of the /o resonance indicate [274] that no /o signal might be extracted from the K'^K~ invariant 
mass spectrum in pp — > ppK^K~ reaction at near-threshold energies, due to the large contribution 
from other reaction channels [79] as arising from pion and kaon exchange. 

Considering the ao (980) as a conventional qq meson, predictions on the total and differential cross 
sections for ao production via the reaction pp —>■ dag are given for the COSY energy range in [275]: 
Within a two-step model, where the Oq is found to be produced predominantly from u-channel tt 
exchange with a subsequent fusion of the nucleons to a deuteron bound state, a rather large total cross 
section is concluded. An experimental verification at the ANKE facility [276] might give an additional 
hint on the qq content of the isovcctor scalar resonance. The measurements also aim at determining 
the rather poorly known gq branching ratios to the KK and rjTf final states (see discussion in [275] and 
references therein). 

The role of final state interactions in the reactions pp —>■ dK^K^ and pp — > dn^r] is explicitly studied 
in [277]: Considering primary production amplitudes within meson-baryon chiral perturbation theory, 
K^d final state interactions — due to the proximity of the A (1405) to the K^n system — arc expected 
to change the cross section by up to an order of magnitude and to significantly infiucncc invariant mass 
distributions. The K^K^ and Tr+r^ final state systems are investigated within a coupled-channel chiral 
unitary approach, in which the oq (980) is generated dynamically as a meson-meson resonance. Once the 
two remaining independent model parameters — related to the couplings of the two pseudoscalar meson 
systems — are fixed by sensitive experimental data, i.e. K^K^ and K^d invariant mass distributions, 
the approach will provide absolute predictions for the r]7r system. Thus, the underlying dynamics of 
the coupled-channel system giving rise to the ao (980) should be probed directly by the experimental 
approach of [276]. 
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Further experimental efforts on studying the iso-singlet /o and the iso-triplct Uq scalar resonances 
are continuing at the COSY-11 [278,279] and ANKE [280] facihties at COSY using different approaches. 



4.4 Associated strangeness production 

In elementary hadronic interactions with no strange valence quark in the initial state the associated 
strangeness production provides a powerful tool to study reaction dynamics by introducing an "im- 
purity" to hadronic matter. Thus, quark model descriptions might be related to mesonic or baryonic 
degrees of freedom, with the onset of quark degrees of freedom expected for kinematical situations with 
large enough transverse momentum transfer. 

A large fraction of the nucleon spin is carried by the intrinsic orbital angular momentum of quarks 
which - in the condensate - are pairwise in a ^Pq state. The challenge of the strange quark, with its 
intermediate mass between light and heavy quarks, is to explore the effective degrees of freedom via the 
associated strangeness production in the low energy structure of QCD. The qq (ss) vacuum condensate 
is the link between current quark and gluon of QCD and valence quark and potentials of the quark 
model, structures which can give the massless QCD quarks their constituent masses. 

In close-to-threshold kaon-hyperon production, effects of the low energy hyperon-nucleon interac- 
tion are inherent to the observables, and allow to constrain hyperon nucleon interaction models. 

As the strange quark is heavy compared to up and down quarks (m„, m^ << m^ ~ ^qcd) it is not a 
priori clear, whether the strange quark is appropriately treated as a light quark on the hadronic scale in 
chiral perturbative approaches with an expansion parameter of mx/47rF7r = 0.43 (with the kaon mass 
rrix and the pion decay constant F^^) compared to mT^/AnFj^ = 0.12. In exploratory calculations, SU(3) 
baryon chiral perturbation theory only has been applied to kaon photo- and electroproduction [281], 
for a more detailed discussion see [282]). 

With a significant contribution of strange sea quarks to the nucleon wave function being suggested 
by complementary experimental approaches during the past years (see section 4.2 and [198-201]) the 
creation of an ss pair in the nucleon-nucleon interaction should provide a further piece of valuable 
information on the structure of the nucleon [193]. Especially polarization observables are expected to 

act as very sensitive tools [283]. 

At the same time, both the kaon and antikaon (see section 5.4) production cross section in nucleon- 
nucleon scattering are important input parameters within the framework of transport model calculations 
used to describe the strangeness flow in heavy ion collisions: Due to strangeness conservation, kaons 
will not be absorbed in the nuclear medium. Thus, kaon yields in subthreshold production are expected 
to probe the early stage of the collision, i.e. hot, dense matter, the nuclear equation of state at high 
densities [284] and the in-medium kaon potential [285-287], see also section 5.4), which relates to the 
partial restoration of chiral symmetry [288,289] (for a review see [290]). Elementary cross sections are 
also of interest in view of the recent inclusive subthreshold measurements in nucleus-nucleus collisions, 
resulting in comparable K'^ and K~ yields at the same energy per nucleon below threshold [291-293]. 

Exclusive data close-to-threshold on the elementary kaon production via pp — > pK^A have only be- 
come available recently with the advent of dedicated experiments at COSY. Figure 32 shows the available 
total cross section data obtained at the COSY-11 [113-115] and COSY-TOF facihty [116-118] together 
with an earlier result from the BNL Cosmotron [119] and theoretical calculations and parametrizations 
as a function of the excess energy Q. 
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Figure 32: Total cross section of 
the reaction pp — > pK'^K. Data 
obtained at the COSY-11 facil- 
ity are marked by filled [113,114] 
and open squares [115], COSY- 
TOF data by filled circles [116] 
and triangles [117,118] and the 
BNL result [119] by an open 
cross, respectively. Experimen- 
tal data are the same as in fig- 
ure 14b. The arrow denotes the 
opening of the pK^TP thresh- 
old. For comparison, theoreti- 
cal descriptions [294-300] are in- 
cluded (for details see text). 



excess energy Q [MeV] 

Assuming the production amplitude to be constant and both the ph. system and the relative to 
the subsystem to be in s-wave leads to a production cross section proportional to t]'^ (section 2.1, 
eq. (16)). This dependence has been adopted by Schiirmann and Zwermann [294] for a quartic (in terms 
of the maximum meson momentum) parametrization of the cross section (black solid line in fig. 32), 
leading to a reasonable fit [295] of data at excess energies above 150 MeV. Although quite commonly 
employed in the description of nucleon-nucleus and nucleus-nucleus collisions in the past [297], the 
quartic parametrization underestimates the close-to-threshold data by an order of magnitude, as does 
a linear parametrization once suggested by Randrup and Ko (grey solid line [296]), overestimating the 
presently available data at small excess energies by more than two orders of magnitude. It should 
be noted that both descriptions have been estimated from high energy data, before any of the near- 
threshold measurements were done. 

Close-to-threshold, deviations from the energy dependence result from the strong s-wave final 
state interactions in the three two-body pK^, pA and A'+A subsystems. As the strong interaction 
in the pK~^ and AT+A systems appears to be weaker by an order of magnitude as compared to the 
proton-A case [301,302], an accurate description of the energy dependence is expected to be achieved 
by considering the strong attractive pA-FSI and repulsive Coulomb corrections in the proton-kaon 
subsystem. Following the discussion in section 3, a reasonable fit of the energy dependence is obtained 
by factorizing the reaction amplitude Mpp_^pK+A according to eq. (30): 

~ Ffsi (Ppa) Fcoui {PpK+) ■ \Mof ■ Fisi (44) 

Here, Mq is the short-range production amplitude, which is assumed to be constant in the energy range 
close-to-threshold, Mpsi and Fisi describe the elastic interaction in the final state and the reduction 
due to the pp initial state interaction (see section 3). According to the above discussion the interaction 
in the three-body final state Mpsi is approximated by Fcoui{VpK+) ^^^1 FFsi{l>ph)-, which denote the 
corrections for the Coulomb interaction and the pA FSI, depending on the relative momenta Pp^+ and 
Ppy^ in the proton-kaon and proton-A centre-of-mass systems, respectively. The Coulomb correction is 
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essentially given by the Coulomb penetration factor of eq. (28) 



Fcoul {PpK+) — 

with the relativistic Coulomb parameter 

'nc,pK+ = a 



(45) 



SpK+ - ml- ml.+ 



m 



where a is the fine structure constant, SpK+ denotes the total energy in the pK'^ subsystem and the 
triangle function A is defined in eq. (7). 

The modification of the energy dependence close-to-threshold due to the pA FSI 



Ffsi (Ppa) 



PpA (1 + '^m) p^^ + (-l/a,A + r,A p^^ / 2) 



(46) 



is determined by the relative momentum in the pK centre-of-mass system and the pK S-wave phase- 
shift 5pA, which has been expressed in terms of an effective range expansion in eq. (46), with a^A and 
fpA as spin-averaged values of pA scattering length and effective range Following the result for the 
ppr] final state in eq. (38) in section 3 and using relation (21) from section 2.1, the total cross section 
for the reaction pp — > pK'^A writes as: 



cr (pp — > pK'^A) — 



Fisi 



Mn 



dVr 



ps 



PpA+ (-V«pA + rpAPpA/2)' 



(47) 



where F is the flux factor of the colliding nucleons from deflned in eq. (9). Using the energy dependence 
given in eq. (47) an accurate description close-to-threshold is achieved, as shown in figure 14b in 
section 3 (see also [113]). 

Within meson-exchange models, contributions to the reaction dynamics of the elementary strange- 
ness dissociation are expected from direct production (fig. 33a) and from the exchange of nonstrange 
(tt, ?7, p, cr, lS) mesons including N* resonance excitation and (X, i<'*)-exchange in the strange sector 
(fig. 33b). 
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Figure 33: Possible production mechanisms for kaon-hypcron production in nuclcon-nuclcon scattering, 
a) — Direct kaon emission, b) — Nonstrange and strange meson exchange diagrams, c) — Quark- 
gluon exchange mechanism. For the latter see the discussion in section 4.1. 

^^As the pA system can be in a spin singlet ^Sq and spin triplet ^Si state, the amplitudes should be added incoherently. 
However, model calculations [303-305] and experimental data [306,307] suggest that the scattering parameters are rather 
similar for both spin states (fig. 34). 
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Pion- and kaon — exchange were suggested as mechanisms of the associated strangeness production in 
the nucleon-nucleon interaction already by Ferrari [308]. As noted by Laget [310], the contribution from 
direct production (fig. 33a) to the total cross section is small but might be of relevance for differential 
observables as discussed below. 

Modern parameterizations considering tt- and i^T-exchange reproduce the magnitude of the total 
cross section close-to-threshold within a factor of two or three (grey dashed [297] and grey dotted [298] 
hues in fig. 32). In refs. [297,298] the elementary t:N KY and KN — > KN transitions are included 
phenomenologically and both exchange diagrams arc found to be important in the case of the reaction 
pp — > pK~^A, while E production appears to be dominated by vr-exchange However, interference 
terms of the amplitudes as well as the effects of final state interactions have been neglected, although 
the latter have been shown to be of relevance in earlier studies [310]. 

The same holds for a resonance model (dashed hue in fig. 32 [299, 315]) considering tt-, r)- and 
p-exchange exciting baryon resonances up to 2GeV/c^ coupling to the kaon-hyperon channel (fig. 33b 
left), i.e. the Sn N*{1Q50), Pn iV*(1710), P13 A^*(1720) and P33 A(1920) resonances. The fully rel- 
ativistic model essentially uses parameters fixed by the previous study on tt induced strangeness pro- 
duction [311] and experimental branching ratios, interferences of the amplitudes are neglected. Non- 
resonant exchange diagrams, i.e. especially kaon exchange, are not considered from arguments based 
on the author's studies on nN — > KY reactions [311]. As a result, tt exchange and excitation of the 
A^*(1650), coupling strongly to the nN channel, is found to give the dominant contribution close-to- 
threshold and is only exceeded by p exchange at excess energies beyond 1 GeV. However, as FSl effects 
are neglected — according to the authors — the approach is neither expected nor intended to describe 
the close-to-threshold regime [315]. 

Within a factor of two, the total cross section data is described by an alternative approach [300] con- 
sidering TT exchange and excitation of the second Sn resonance, the A^*(1650) (dotted lines in fig. 32) — 
in analogy to 77 production arising from tt exchange and excitation of the first Sn resonance (A^*(1535)). 
Assuming the iV*(1535) and the TV* (1650) to dominate the rjN and K^A systems, respectively, leads to 
identical forms of the production operators and the spin-angular momentum algebra. Thus, after nor- 
malizing the scale of the pp — > pK^A excitation function to low-energy data on the ppr] final state, the 
model gives an absolute prediction for A production, while the energy dependence closc-to-threshold is 
essentially determined by three-body phase-space modified by the pA final state interaction. The two 
corresponding curves shown in figure 32 (dotted lines) differ in the choice of low energy pA scattering 
parameters [304-306] determining the proton-hyperon FSI (see eq. (46)). 

The situation is summarized by Kaiser [237]: Separating the on-shell S-wave final state interaction 
and the invariant T matrix at threshold, the information from the closc-to-threshold total cross section 
can be "condensed" into essentially one number, the threshold transition amplitude parametrizing the 
T matrix, with the energy dependence being accurately reproduced by three-body phase-space and 
FSI (see also our discussion in section 4.1). In an analysis based on tree-level meson exchange diagrams 
pointhke u exchange alone or the total vector p°, u and K* exchange, are found to describe the transition 
amplitude However, also pseudoscalar amplitudes, i.e. pointlike exchange alone, or — in line 
with [300,316] — tt exchange followed by an excitation of the Sn A?"* (1650) resonance when added to 

^■^For TT induced strangeness production, either experimental values [309] (as in the work of Lagct [310]) or paramctriza- 
tions (cross sections determined from the resonance model approach [311] used by Li and Ko [297], and Sibirtsev [298]), 
respectively, are employed. KN scattering data are available experimentally [312] and in terms of phase-shift analy- 
ses [313]. 

^''The relative strength of the K exchange mechanism in A and TP production reflects the smallness of the KNTi 
couphng constant gKNT, ^ compared to Zknk (see e.g. [305,314]) 

^^The transition amplitudes extracted in [237] based on the approach developed in [34] arc related to the production 
amplitude \Mq\ used in this work by simple kinematical factors, i.e. ) Mq [^ = ) T]^ • (Smpm^ • \{s,xn^,vcLp)) / s, with m^ 
denoting the mass of the second baryon B = p,Y m the final state and ) T ) as threshold transition amplitude, i.e. ) A \ 
and \ C\ for TT and 77 production in [34] and ] K ]/V3 in [237]. 
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the total pseudoscalar and vector exchange contributions are in good agreement with the experiment. 
It is concluded [234] that the unpolarized pp — > pK'^A total cross section close-to-threshold does not 
provide enough information to determine the underlying production mechanism. 

Instead, further information should be provided in the future by angular distributions and polariza- 
tion observables, as reported at higher excess energies by the DISTO collaboration [317]: The impor- 
tance of K exchange for A production, which might already be expected from elementary amplitudes 
has been supported experimentally by data on the normal spin transfer coefficient D^vat in exclusive 
data at an excess energy of Q = 430 MeV. Djviv is found to be large and the A polarization correlated 
with the beam polarization is oriented opposite to the beam spin for forward A production. Based on 
the model described in [310], calculations by Laget show Da? at to discriminate between the tt and K 
exchange mechanism from a different spin coupling at the vertices between the polarized proton in the 
initial and the hyperon in the final state While from angular momentum and parity conservation tt 
exchange leads to Dnn = +1 for forward production, kaon emission results in a spin flip (Datat = — 1) 
and consequently, the latter is clearly favoured by the experimental result [317] Recently, the nega- 
tive sign of the spin transfer coefficient, i.e. the dominance of K exchange for A production, has been 
conffimed for two different beam momenta corresponding to excess energies of 319 MeV and 200 MeV, 
respectively [321]. However, unexpectedly large negative values are also seen for backward production, 
where the hyperon is preferentially associated with the unpolarized target proton — an observation, 
which is so far not reproduced within the meson exchange picture. 

The ffist close-to-threshold measurement of the A recoil polarization at the COSY-TOF facility has 
shown a negative polarization for transverse momentum transfers pj > 0.3 GeV/c [116] which implies 
evidence for S- and P-wave interference terms at the respective excess energy of 138 MeV according 
to [237]. 

Differential observables, like the occupation in a Dalitz plot, or in projection invariant mass distri- 
butions of two-particle subsystems, might also provide important information to identify the relevant 
production mechanism. A resonant production should influence for example the AK~^ invariant mass 
distribution significantly at excess energies in the range of 100 MeV, as calculated within the resonance 
model approach of [299,315] by Sibirtsev et al. [323] First data has already been pubhshed [116] 
and higher statistics distributions at several excess energies will be available from data taken at the 
COSY-TOF facility in the near future. 

A production on a neutron target, i.e. via the reactions pn — > p{n)K^{K~^)A, might turn out to 
be sensitive on the relative strength of tt, p and r] exchange, as it has been discussed in case of rj 
production [49,300]. 

An important aspect in the theoretical interpretation of data close to the production threshold is 
the interplay between the type of hyperon-nucleon interaction employed and the result of the model 
calculation. As already discussed by Laget [310], with experimental observables being significantly 
influenced by the strong final state interaction at low transverse momenta, exclusive close-to-threshold 
data should be expected to eventually constrain strongly YN low-energy scattering amplitudes. Within 
a meson-exchange approach, i.e. considering direct kaon emission as well as tt and K exchange (figs. 33a- 
b) and using an off-shell T matrix description for the hyperon-nucleon interaction the sensitivity on 
the microscopic hyperon-nucleon potential used has been studied in detail by Kelkar and Jain [325] 
for A production in comparison to the inclusive SATURNE data from [326]. The authors conclude, 
that — in order to separate different channels contributing at the SA^ threshold — i) exclusive 

^^The total irp KY cross section is in the range of Imb compared to 10 mb for K^p scattering (see fig. 10a). 
^''For a more general discussion on the interpretation of spin transfer measurements see [318]. 

^^The DISTO results show an opposite sign c;oniparcd to inclusive Djvjv high energy (Q = 16.9 GeV) data [319], which 
have been accounted for within quark model calculations [320] . 

^^A negative A polarization is also observed for pj > 0.4GeV/c in proton-proton scattering at higher energies [322]. 
^•^The calculations on differential distributions have recently been extended to the reaction — > pK^YP in [324]. 
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data is needed (which shall be discussed below) to constrain hyperon-nucleon interaction models and 
that ii) an improved knowledge on the hyperon-nucleon interaction is prerequisite to further detailed 
investigations on the relevant production mechanisms. We shall use these conclusions as guideline for 
the discussions during the next paragraphs. 

Due to the short hyperon lifetimes direct YN scattering experiments are difficult to perform and, con- 
sequently, low-energy scattering parameters arc still rather unknown experimentally. In the case of the 
A-proton system, experimental results have only been available at excess energies above 3.8 MeV [306, 
307]. In contrast, a study in close-to-threshold production via the strong final state interaction covers 
an excess energy range down to threshold, otherwise inaccessible for elastic scattering experiments. 




Figure 34: Ap S-wave singlet (open symbols) and triplet (closed symbols) scattering length a^A and 

effective range Tpj^. Experimental data are from Ap elastic scattering (stars [306] and triangles [307]), K~ 
capture on a deuteron (diamond [327]) and from the FSI approach in threshold production (cross [328], 
spin-averaged value). Results of YN interaction models are included by squares [303] and circles and 
double circles (solutions A and B of [304,305]), respectively. Errors have only been determined for the 
K" capture experiment [327]. 



Assuming the energy dependence of the reaction pp pK~^A close-to-threshold being determined 
by eq. (47), i.e. by three-body phase-space modified by the strong pA FSI and the Coulomb repulsion 
in the pK'^ subsystem, constraints on the spin-averaged low energy Ap scattering parameters have 
been extracted in [328] from the differential occupation of Dalitz plot distributions obtained from 
measurements at the COSY-11 facility [113]. However, in the analysis scattering length and effective 
range appeared strongly correlated, as inherent to the parametrization used in eqs. (46) and (47). 
Thus, only a combined fit using the complementary information of low-energy Ap elastic scattering 
experiments [306,307] constrained by the Dalitz plot analysis allowed to derive separate spin-averaged 
values [328] . The result — new spin averaged values of Ap S-wave scattering length and effective range 
{ap\,rp\) — (— 2.0fm, 1.0 fm) — is shown in figure 34 together with other experimental and theoretical 
values, where the extracted value of the Ap scattering length is clearly fixed by the old data. 

Although the FSI approach to low-energy A-proton scattering proves to be a powerful tool, it should 
be noted that the above values from [328] have been derived assuming a cross section ratio according 
to the number of possible magnetic spin quantum numbers, i.e. a ratio of three to one between the 
triplet and singlet Ap final states, which was confirmed in the early bubble chamber experiments. 
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Consequently, the analysis implies that the underlying reaction mechanism in proton-proton scattering 
does not favour either the singlet or triplet state, different to the assumed ratio of spin states. Future 
threshold production experiments employing both polarized beam and target should dispel this concern 
and allow to extract singlet and triplet hyperon-nucleon scattering parameters separately. 

Eventually referring to the first conclusion of Kelkar and Jain [325] cited above, further constraints 
on the present understanding of close-to-thrcshold hyperon production have arisen from data on the 
reaction pp — > pi^+E^ taken at the COSY-11 installation [114]: At equal excess energies — up to 13 MeV 
with respect to the hyperon production thresholds — the total cross sections have been determined for 
both A [113,114] and S° [114] production. A fit to the data including the hyperon-nucleon FSI is shown 
in figure 35 and has been performed using the parametrization of the energy dependence developed by 
Faldt and Wilkin [110, 120,300] as given by eq. (35), yet neglecting the weak variation of the cross 
section due to the variation of the flux factor F. It should be noted that Gasparyan [329] concludes that 
the proton-hyperon FSI is less important for the T^N rather than for the AN channel. A similar result 
gives the preliminary analysis of new data taken at COSY-11 for both reaction channels up to excess 
energies of Q = 60 MeV [115]. 
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Figure 35: Total cross sections of the reactions pp — > pK^K/TP [113, 114]. In comparison, fits of the 
energy dependence for both hyperon production channels using a pure phase-space calculation (dotted 
hues) and phase-space modified by the proton-hyperon FSI (solid hues, eq. (35) [110, 120, 300]) are 
shown. 

The most striking feature of the data is the observed TP suppression in the energy range close-to- 
threshold with 

_ a{pp^pK+A) _ Q 
^^/"-a(pp-pir+E0)-2^-^' (48) 

while at excess energies > 300 MeV this ratio is known to be about 2.5 [238]. Considering only tt 
exchange, data on n induced hyperon production via nN — >• KA (S°) result in a ratio of 7^a/s ~ 
0.9 [330], clearly underestimating the experimental value of relation (48). Kaon exchange essentially 
relates the ratio TZ\fs to the ratio of coupling constants squared SNAK/sifEK- Although there is some 
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uncertainty in the literature on their values [304,305,314,331,332], a A/S° production ratio of 27 
follows from the suitable choice of the SU(6) prediction [232,333], in good agreement with experiment. 
However, effects of final state interaction as well as the importance of tt exchange especially for E° 
production are completely neglected by this simple estimate. 

Inclusive production data taken at the SPES 4 facility at SATURNE in proton-proton scattering 
at an excess energy of Q = 252 MeV with respect to the pK'^K threshold show enhancements in the 
invariant mass distribution at the Ap and TiN thresholds of similar magnitude [326]. However, with 
only the kaon being detected, it is not clear whether the enhancement at the SA^ threshold is due 
to S production. Qualitatively, a strong EA^ Ap final state conversion might account for both 
the inclusive results from SPES 4 as well as the depletion in the COSY-11 data. In production 
experiments, evidence for EA^ Ap conversion effects is known experimentally from exclusive hyperon 
data via K~d — > 7r~Ap [327,334,335] (for a theoretical interpretation see [336,337]). In hypernuclear 
physics, the strong EA^ AN mode appears as the dominating decay channel of E-hypernuclei [338], 
and the EA^ interaction reflects in both shift and broadening of the last X-rays observed from E~ 
atoms [339]. Theoretically, modern hyperon-nucleon interaction models show a significant cusp effect 
for the AA^ cross section at the EA^ threshold [304,314]. 

In exploratory calculations performed within the framework of the Jiilich meson exchange model [330] , 
taking into account both tt and K exchange diagrams and rigorously including FSI effects in a coupled 
channel approach, a final state conversion is rather excluded as dominant origin of the observed E*^ 
suppression. While A production is found to be dominated by kaon exchange — in agreement with 
the DISTO results at higher excess energies discussed above — both n and K exchange turn out to 
contribute to the E° channel with similar strength. It is concluded by Gasparian et al. [330], that 
a destructive interference of tt and K exchange diagrams might explain the E'^ suppression close-to- 
threshold and a good agreement with the COSY-11 data is obtained after including the overall reduction 
Fjsi in eq. (47) due to the pp interaction in the initial state [19]. 

An experimental study of E production in different isospin configurations should provide a crucial 
test of the above interpretation: For the reaction pp — > nir+E+ an opposite interference pattern is 
found, i.e. the nK'^T:'^ channel is enhanced for a destructive pattern as compared to the construc- 
tive interference of K and tt exchange. For the choice of a destructive interference favoured by the 
p/T+E^ channel, the corresponding ratio o" (ni^+E"*") /a{pK^TP) ^ 3 is in good agreement with high 
energy data, whereas the alternative choice is rather difficult to reconcile with existing data at higher 
energies [330]. In contrast, for the reaction pp — > pK^T.'^ the same interference pattern occurs as for 
the pK^YP channel, resulting in a ratio cr (pir°E+) /a{pK^YP) ^ 3.3 close-to-threshold for either 
choice [340], while a ratio close to unity is found in the literature for higher energies [238]. 

As stated by the authors of ref. [330], contributions from direct production (fig. 33a) as well as 
heavy meson exchanges (fig. 33b) have been neglected so far in these calculations but might infiuence 
the A/TP production ratio (see also [237,299,315]). In fact, employing both a tt and K exchange based 
meson exchange model without any interference of the amplitudes following [298] and the resonance 
model of [299,315] the close-to-threshold A/Tf production ratio has been studied in [341]. Within 
a factor of two of the experimental error bars in [114] both models — with parameters fixed at data 
from excess energies > 1 GeV — are in reasonable agreement with the data, i.e. — according to the 
authors — the total cross section data close-to-threshold are not sensitive on the details of the model 
used [341] 

The latter study has been critically discussed in [342], where the strangeness production is modeled 
in an effective Lagrangian approach following [107, 316, 343] and strangeness production proceeds via 

^^It should be noted, that unhke [330] tt and K exchange amphtudes are concluded by Sibirtsev et al. [341] to contribute 
to both A and E° production with similar magnitude. As mentioned in [341], an experimental determination of the kaon 
exchange contribution close-to-threshold — e.g. from polarization observables as in [317] — should be crucial to identify 
the dominant reaction mechanisms. 
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TT, p, u) and a exchange exciting the nucleon S-wave resonance A'^*(1650) and the P-wave resonances 
A^*(1710) and A^*(1720). For both the A and E° production channel, tt exchange followed by an 
excitation of the A^*(1650) is found to dominate the total cross section close-to-threshold (see also [300]), 
though an on-shell coupling of the A^*(1650) to the K^Tf channel is suppressed with the channel 
opening at 1686MeV/c^. In contrast, the influence of the iV*(1650) on TP production and possible 
interference effects of resonance contributions have been neglected in [341] At larger excess energies 
(Q > 300 Me V) the A'"*(1710) is concluded in the work of Shyam et al. [342] to dominate both reaction 
channels, which would be expected when considering coupling constants only 



It is worth noting, that the OBE calculations performed by Laget [310,346], although not using the 
interference of pion and kaon exchange graphs explicitely and just choosing the relative sign to maximize 
the cross section, not only describe the recent close-to-thrcshold data on the h./TP ratio within a factor 
of two, as well as the polarization transfer results of the DISTO experiment [317], but also give an 
accurate description of the YN invariant mass distributions obtained in the inclusive measurements at 
SATURNE [326] cited above: In particular, the FSI contribution from the direct Ap — > Kp amphtude 
is found to significantly enhance the kaon spectrum just above the h.p threshold, while the strong 
EA^ — > Ap coupling "induces only a small dip close to the SA^ threshold" in the tt and K exchange 
amplitudes [310]. On the other hand, it is this strong channel coupling and the resulting rapid variation 
of the ^Si AA^ amplitude near the SA^ threshold, which appears responsible for the sharp rise of the 
invariant mass distribution obtained at SPES 4 at the EA^ threshold via the "kaon direct emission 
amplitude" , i.e. the coupling of one proton to a K^Y pair and a following interaction with the second 
proton (fig. 33a). 



^^It has been remarked by Shyam et al. in [342] that the final state interaction employed in [341] is "at variance" 
with [108] and other approaches in the literature [33,316,344]. 

^^Accordingly, close-to-threshold production in the resonance model of [315] proceeds preferrably via tt and r] 
exchange, the latter being enhanced for the TP compared to the A channel due to the large branching ratio of the 
N*{\1\Q) to the Nt] channel. However, as pointed out in [324] the r] exchange contribution is tainted with rather large 
uncertainties due to different values of the t}NN coupling constant extracted from experimental data (see [345] and 
section 3). 
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Figure 36: A/E'^ production ratio in proton-proton collisions as a function of the excess energy. Ex- 
perimental data within the shaded area, which corresponds to the range described by relation (48), are 
taken from [114], data at higher excess energies from [115]. Calculations [347] are performed within the 
Jiilich meson exchange model, assuming a destructive interference of K and tt exchange and employing 
the microscopic YN interaction models Nijmegen NSC89 (solid hne [314]) and the new Jiilich model 
(dotted line [348]), respectively. 



Recently, measurements of the A/S° production ratio in proton-proton collisions have been extended 
up to excess energies of Q = 60 MeV at the COSY-11 installation [115], in order to study the transition 
between the low-energy enhancement of the ratio reported in [114] and data at excess energies larger 
than 300 MeV. In comparison to the experimental data, in figure 36 calculations are included obtained 
within the approach of Gasparian et al. [330] assuming a destructive interference of tt and K exchange 
with different choices of the microscopic hyperon nucleon model [347]. 

As it was already emphasized in [330], the result depends on the details — especially the off-shell 
properties — of the microscopic hyperon-nucleon interaction employed, although the actual choice does 
not alter the general result in [330] of only the destructive interference of K and tt exchange explaining 
the experimentally observed suppression of the S° signal close-to-threshold. At the present stage of 
theoretical investigations, both the good agreement found for Jiilich model A [304] with the close-to- 
threshold result for TZ\f£ in relation 48 and for the Nijmegen model (solid line in fig. 36 [314]) with the 
energy dependence of the cross section ratio in figure 36 should rather be regarded as accidental. In 
the latter case an SU(2) breaking in the ^Si UN channel had to be introduced [314]. Consequently, the 
relation between the E^'p amplitude with the S+p and E~p channels becomes ambiguous. Only one of 
the choices leads to the good agreement with the data shown in figure 36 (solid line), whereas the other 
one results in a completely different prediction [349]. 

Calculations using the new Jiihch model (dotted hne in figure 36 [348]) do not reproduce the tendency 
of the experimental data. It is suggested in [347] that neglecting the energy dependence of the elementary 
amplitudes and assuming S-waves in the final state might no longer be justified beyond excess energies 
of 20 MeV. 

Conclusions from figure 36, as well as considerations by Kelkar and Jain [325] for the pp — > pK~^A 
reaction clearly demonstrate, that — once the reaction mechanism for close-to-threshold hyperon pro- 
duction is understood — exclusive hyperon production data should provide a strong constraint on the 
details of hyperon-nucleon interaction models. 

Finally, in future high resolution hyperon production experiments close to threshold may give further 



64 



information on the possible existence of low-lying strange dibaryons. Narrow, stable six-quark objects 
Dg and Df with a gr^ (g) cluster substructure and strangeness S = — 1 coupling to P-wave nucleon- 
hyperon systems have been postulated from an extended MIT bag-model [350] below and above the 
TiN thresholds with widths in the range of 10 keV to 10 McV. Experimental signatures consistent 
with a strange dibaryon hypothesis have been found both in strangeness transfer reactions K~d — > 
n-YN [327,351,352] and in the inclusive pp K+X data taken at SATURNE [326] at excess energies 
of 175 and 306 MeV with respect to the TPp threshold. However, the statistical accuracy in these 
experiments turned out to be too small for unequivocal conclusions. For a high resolution study of the 
inclusive pp — > K^X reaction channel new experimental efforts arc on their way at the BIG KARL 
spectrometer at COSY, with a resolution improved by a factor of five and a higher statistical accuracy 
compared to the SATURNE measurements [353]. 

5 Meson pair production 

5.1 Double pion production in NN scattering 

Single and double pion production in elementary vrA^, and A^A^ reactions is an important source 
to gain information on nucleon-nucleon, nucleon-pion and pion-pion interactions. Apart from these 
fundamental aspects, the pion production on nucleons also enables the study of the properties of nu- 
cleonic resonances, which might be excited in the reaction processes. The reaction — > Ntttt for 
instance has been found to be influenced by the formation of the A^*(1520) resonance [354] with a sub- 
sequent decay into Avr. Therefore, studies on this reaction channel allow to gain information about the 
A^*(1520) — A(1232)7r decay amplitudes. Additionally, contributions of the A^*(1440) resonance, fol- 
lowed by its decay A"* (1440) — > N{7ni)Jz^^-^^, play a non-negligible role in the near-threshold production 
of pion-pairs via the reaction channels 7A^ — > Ntttt and ttN — > Ntttt [355] . 

In case of the pion-pair production in the nucleon-nucleon scattering similar contributions of nu- 
cleonic resonances are under discussion. Recently, detailed calculations on the two-pion production 
in nucleon-nucleon interactions have been performed [356], considering non-resonant terms as well as 
contributions from nucleonic resonances with subsequent decays into a nucleon and pions. One of the 
main results of these studies is the prediction of a dominant effect of the A^*(1440) Roper resonance, 
followed by its decay A^*(1440) — > ^{'^'^)s-wave^ region of low excess energies. According to this, 

the absolute scale of the low energy total cross section of reaction channels, where the pions can be in 
an isospin zero state, e.g. pp — > ppTT'^7T~ and pp — > ppn^n^, should be determined by the excitation of 
this resonance. In this region of excess energies contributions of non-resonant terms are discussed to 
be negligible for the description of the excitation function of both these reaction channels. However, 
at higher energies (Q > 200 MeV) the excitation of A resonances has been evaluated to dominate the 
total cross sections. 

A different situation is given in the case of the pp pmr^TT^ reaction channel. Here contributions 
of the Roper resonance with the subsequent decay into N {tt'k)'^~^^^^ are forbidden by conservation laws 
and contributions of remaining A^* production diagrams are suppressed at low excess energies. Since the 
non-resonant terms are also negligible in this case, the dominating processes are given by A resonance 
excitations. 

It is worth noting that the discussed model calculations yield predictions for the excitation functions 
for different two-pion production reactions from threshold up to high excess energies (Q ~ 300 MeV) and 
all arc based on the same Feynman diagrams. However, especially in the close-to-threshold region the 
data situation for most double pion reaction channels is far from being optimal. Recently, three isospin- 
independent reaction channels on the two-pion production in proton-proton interactions, pp — > ppTT'^TT~, 
pp — > ppTT^TT^ and pp — > priTT'^TT^, have been studied in high statistics precision measurements at the 
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PROMICE/WASA facility at CELSIUS and resulted in the first total and differential cross sections in 
the near-threshold region [357, 358] . 

In fig. 37a the total cross sections for the 7r"'"7r~ production in proton-proton scattering are presented 
up to an excess energy of Q ^ 300 MeV. The solid and the dashed lines represent the above discussed 
predictions of ref. [356] for this reaction channel and correspond to two possible solutions in the applied 
model. Obviously these calculations are in good agreement with the experimental data both close-to- 
threshold and in the region of higher excess energies, while the model seems to underestimate the data 
at intermediate energies. It should be noted that the data points from Dakhno et al. [362] have been 
obtained using a deuterium filled bubble chamber and selecting events from the quasi-free two-pion 
production pp{n) — > pp{n)7T~^n~ . However, in the analysis of the data the Fermi motion of the nucleons 
in the deuteron has been neglected. In later calculations an adequate consideration of the deuteron wave 
function has been found to shift the effective excess energy by ~ 20 MeV towards higher energies [365], 
corresponding to a shift of the data towards the region expected by the discussed model calculation. 
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Figure 37: Total cross section data for the reactions: 

a) pp — >• pp7l~^TT~ , 

b) pp — > ppTT^TT*^ and 

c) pp — > pmr'^Tr'^. 

The solid and the dashed curves represent model 
calculations from [356]. 
The data were taken from [357-364]. 
A version of the article including a figure of better qual- 
ity can be found at the COSY-11 homepage: 
http : / /ikpellOl.ikp.kfa — juelich.de/cosy — 11 
/pub/ List-of -Publications. html^papers 
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Figure 37b displays the experimental situation for the reaction channel pp — > pp7r°7r°. The new data 
from PROMICE/WASA disfavour the solution of the model calculation indicated by the solid hne but 
are in very good agreement with the prediction shown by the dashed line. However, the high energy 
data from ref. [363] are in disagreement with both curves and are clearly overestimated. 

Finally, the close-to-threshold total cross section data of the reaction channel pp — > pnTr+vr^ are 
presented in fig. 37c. 

Different to the previous situation now the data is underestimated by both solutions of the model 
calculations, especially in the near-threshold region. Since now a production via the A'"*(1440) 
^(^^)sXave is forbidden, this observation might indicate contributions from processes like NN — > 
A^A^* —>■ NAtt NNtttt larger than predicted by theory. 

It should be mentioned that contributions of higher partial waves and final state interactions have 
been neglected in the model calculations, indicated by the presented curves. However, the effect of 
the pp FSI has been investigated and is reported to increase the close-to-threshold cross sections by 
nearly one order of magnitude, while at Q ~ 120 MeV this effect is reduced to an increase of a factor 
of ~ 2 [366]. According to this, considering the interaction of the final state baryons is expected to 
reduce the apparent discrepancy between the n'^n^ total cross section data and the corresponding model 
calculations. Surprisingly, in case of the tt+tt" and tt^tt^ pair production, the near-threshold data are 
described reasonable by these calculations neglecting the pp FSI. 



Furthermore, signals for the contribution of higher partial waves are reported for the near-threshold 
data [357], which also have been neglected in the discussed model calculations. This observation might 
be interpreted as a signal of a heavy meson exchange (ex, p) between the interacting protons. 



5.2 Double pion production in pd scattering 



It is well known and cited in the literature as the ABC-anomaly [367, 368] that the missing mass 
enhancement at around 310 MeV with a width of 50 MeV observed in the inclusive measurements of the 
pd ^HeX^ reaction at a proton beam energy of Tp > 745 McV (corresponding to an excess energy 
of Q > 190 McV for the tttt production) is not excited in the ^"H final state. Detailed investigations 
confirm that this observation is to be associated with the isospin-zero s-wave tttt double-pion state 
but additionally stress that the missing mass peak and its width vary with the beam momentum [369]. 
Further, since the isoscalar tttt scattering length is small [370] the observed anomaly is supposed to be 
of kinematical origin and associated with the intermediate excitation of two A resonances [371]. 
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Figure 38: Differential cross sections for tlie pd ^He tt+tt^ reaction at Tp = 546 MeV as a function 
of four kinematical variables. Dashed curves represent phase-space predictions normalized to the data 
whereas the solid curves are calculations assuming the tt'^t:' system to emerge in a relative p-wave, 
see [372]. 

Recently the MOMO (Monitor-of-Mesonic-Observables) collaboration at COSY investigated the 
TTTT production closer to threshold in an exclusive measurement of the pd "'iJe tt+tt" reaction at a 
beam energy of Tp = 546 MeV (excess energy Q = 70 MeV for the tttt production), where the ABC 
enhancement is expected to show up near the centre of the available TTTT-invariant mass distribution [372]. 

In figure 38a-d differential cross sections in terms of four kinematical variables are presented: 

a) The tt'^tt' excitation energy T^tt = T^nn — 2m^, where m^Tr is the two-pion invariant mass. 

b) The Ti-^He excitation energy. 

c) The angle (^^) between the proton and one of the tt's in the overall centre-of-mass (cm.) system. 

d) The angle (^^ttttp) between the relative 7r"'"7r~ momentum and the beam axis in the overall cm. system. 

The Ttttt distribution at this near-threshold data is pushed closer to the maximum values of excitation 
energy and is in contrast to the original ABC result [367,368] which showed an enhancement over phase- 
space at Ttttt ~ 30 MeV. The distribution of the tt- ^He excitation energy is rather consistent with the 
phase-space prediction. The significant anisotropy of confirms the importance of higher partial waves 
involved in the reaction process. The MOMO collaboration concludes [372] that the agreement between 
the Ttt-w data of fig. 38a with the solid line (phase-space multiplied by Ttttt) and the fair linearity of 
sin^ Ot^t^p indicates that the 7r"'"7r~ system is essentially produced with an internal angular momentum 
Z = 1. 

On the other hand, as outlined in section 5.1, the Valencia group [356,373] developed a model for 
the two pion production resulting in a dynamical origin for the small invariant masses of the low energy 
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part, where the two pions can be in an isospin zero state. Here the reaction dynamics is dominated by 
the intermediate excitation of the A'"*(1440) Roper resonance followed by the N* — > Att — > N {7r7r)Jz^ave 
decay. 



The experimental situation seems to be the following: s-wave tttt production has been observed in the 
pd ^He TT+TT^ reaction very close-to-threshold [374] and at high excess energies (Q ~ 200 MeV) [367, 
368] whereas at an intermediate excess energy range a p-wave n'^n^ production process dominates. 
Similar results have been obtained in np — > dX [375,376] and in the ir'^d — > Tr'^Tr'^nn {ir^n'pp) [377,378] 
reactions. 



5.3 Double K— meson production in pd scattering 



After the remarkable results of p-wave production in the tt+tt system the MOMO collaboration con- 
tinued equivalent measurements for the pd — > ^He K^K~ reaction at excess energies of Q = 56 MeV, 
40 MeV and 35 MeV. 



The K^K~ invariant mass spectra in units of the K^K^ relative energy 1k+k- ai'c shown in 
figure 39. A clear signal of the 0-meson production is seen above the continuous non-resonant K^K~ 
spectrum. The dashed lines in all three parts of the figure correspond to simple phase-space calculations, 
where the Q-dependences of the cross sections {ax+K- and a^) were taken to scale as ax+K- oc 
for the three-body and o",^ oc Q^^^ for the two-body final state, respectively. The relative angular 
momentum between the two kaons was assumed to be / = 0. No significant deviation of the data from 
these summed phase-space distributions is observed and thus it is interesting to note, that no indication 
of p-wave in the K'^K~ system is present in the mass spectra and angular distributions, in contrast to 
the above discussed data for the two pion production via the reaction pd — > ^HeTr'^Tr" . 
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Figure 39: K+K' 
invariant mass spec- 
tra from the reaction 
pd ^HeK+K- at 

g = 56 MeV, 40 MeV and 
35 MeV above threshold 
plotted in units of the 
K'^K~ relative energy 
Tk+k-, see [379]. 



5.4 Strangeness dissociation into KK— pairs 

Measurements on the KK meson pair production are of considerable interest in the context of the 
continuing discussion on the nature of the scalar resonances /o(980) and ao(980), as already discussed 
in section 4.3. 

Available close-to-threshold data [77, 78] on the elementary antikaon production channel in proton- 
proton scattering, pp ppK^K^ , are shown in figure 40. In comparison, parametrizations of the energy 
dependence of the total cross section based on a four-body phase-space behaviour and for a resonant 
production via the /o are included. Obviously, within the experimental error bars and neglecting the 
possible influence of higher partial waves the data is consistent with both a non-resonant and 
resonant production, as already concluded from figure 31 in section 4.3, i.e. from the K^K" invariant 
mass distribution obtained at an excess energy of Q = 17 MeV [273]. 



^^Preliminary data by the DISTO collaboration [380] indicate that K^K production at an excess energy of Q = 
114 MeV is still consistent with isotropic emission. 



70 




s m 10 



Figure 40: Total cross sections for the free K'^K~ pair production in proton-proton collisions [77, 
78]. The sohd hnc indicates calculations on the basis of the non-resonant four-body phase-space 
expectations fitted to the data points. The dashed lines correspond to three-body phase-space for the 
resonant K^K~ production via the /o, normalized to the COSY-11 data point). The pp FSI and the 
Coulomb interaction have been taken into account. The effect of the large uncertainty about the width 
of the /o(980) resonance (F = 40 - 100 MeV [134]) is indicated by the dashed area. 
A version of the article including a figvirc of better quality can be found at the COSY-11 homepage: 
http : / /ikpellOl.ikp.kfa — juelich.de/cosy — ll/pub/List_of-PublicaUons.htmli^papers 



Furthermore, exclusive K~ production data are of special interest with respect to subthreshold 
kaon production experiments in nucleus-nucleus interactions, which are expected to probe the antikaon 
properties at high baryon density. Recent inclusive subthreshold measurements [291-293] resulted in 
comparable K'^ and K" yields at the same energy per nucleon below the production thresholds for the 
elementary reactions pp — > K'^X, which has been discussed as an indication for a possible restoration 
of chiral symmetry (as reviewed in [290]). 

Available inclusive elementary and total cross sections in nucleon-nucleon interactions are 
shown in figure 41. In contrast to the results obtained in subthreshold measurements, in proton-proton 
scattering the inclusive cross section exceeds K~ production by more than an order of magnitude in 
the near threshold region. To explain this observation, different models [74,286-289,381,382] consider 
kaons {K^,K^) and antikaons {K^,K'^) to be subject of repulsive and attractive forces within the 
nuclear medium, respectively. Recent theoretical studies using a chiral unitary approach to describe the 
K~N S-wave interaction allow for a self-consistent microscopic implementation of medium effects [286, 
287] . The K~ self-energies obtained are consistent with kaonic atom data, as shown in [383] . 
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Figure 41: World data on the 
inclusive total K'^ and K~ pro- 
duction cross sections in proton- 
proton collisions at excess ener- 
gies below Q = 2 GeV. The data 
are taken from [77,113-118,238, 
317,384]. 
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6 Meson production on light nuclei 

Ncar-thrcshold meson production experiments in nuclcon-nucleon interactions might help to obtain 
fundamental information on the elementary production processes and final state interactions. Extending 
those measurements from the binary nucleon scattering to reactions with three or four participating 
nucleons, i.e. meson production in the proton-deuteron or deuteron-deuteron scattering, allows to study 
the role of further participating nucleons. This facilitates a gain in information on possible multi-step 
production processes, which are naturally absent in the nucleon-nucleon case. 



6.1 Dynamics in three and four nucleon systems 

Of considerable interest is the production of rj mesons in the pd — > ^Herj reaction. Studied at the SAT- 
URNE accelerator in the near threshold region [385, 386] this channel is reported to expose remarkable 
features. In spite of the much higher momentum transfer, the observed threshold amplitude for the 
77 production, defined by the unpolarized centre-of-mass cross section 

^ {pd ^ 'Her]) = 4 \fnipd - 'Her])\\ (49) 

was found to be comparable to that for the pd 'Hen^ reaction at its respective threshold [385]. In 
addition, while being consistent with s-wave production, |/t,(p^)|^ decreases by a factor of three from 
threshold up to a centre-of-mass momentum of emitted r] mesons of p* = 70McV/c (=Q ~ 7MeV). 

This observed strong decrease in the near-threshold region can be addressed to a strong rj 'He final 
state interaction and is therefore of special interest in view of the existence of ?7-mesic nuclei. Although 
being the topic of several theoretical investigations [138, 140,387-391], the possibility for the formation 
of quasi-bound 'He rj states is still an open question. 
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Figure 42: Total cross section a and square of the production amplitude / for the pd ^Herj reaction. 
The lines are explained in the text. The data are from references [386, 392, 393] and show only the 
statistical uncertainties. 

A version of the article including a figure of better quality can be found at the COSY-11 homepage: 
http : / /ikpellOl.ikp.kfa — juelich.de/cosy — 11 / pub / List_of -Publications. htmlij^papers 

The experimental situation for the pd ^Herj reaction is presented in fig. 42, displaying total cross 
section data from threshold up to an excess energy of Q = 100 MeV. 

The SATURNE data [386] (circles) expose an energy dependence, which clearly differs from pure 
phase-space expectations: According to naive two-body phase-space considerations, neglecting higher 
partial waves, an excitation function according to o" oc Q^^^ is expected. One promising ansatz to 
describe the reaction process was given by Kilian and Nann [394], taking into account both nucleons 
from the target and considering double-scattering diagrams. A sketch of such a double-scattering 
diagram is illustrated in the left part of fig. 43. 
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Figure 43: Production diagrams for the pd ^He rj reaction. 

Calculations by Faldt and Wilkin [395] are in line with this approach and describe the pd — > ^He rj 
reaction as an initial pp dir'^ {pn — > di:^) production followed by a final reaction according to 
TT'^n — > r]p (tt^p — > r]p). These investigations lead to a prediction of the spin averaged threshold 
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amplitude squared of |/(p* = 0)p = 0.69/xb/sr, underestimating the measured results by a factor of 
2-3. Since none of the considered sub-reactions exposes a significant variation of the cross section in 
the relevant energy region, the observed energy dependence of the production amplitude has been 
interpreted to be a signal of a strong final state interaction of the ?7-nucleus system. 
Introducing the ^Herj final state interaction by an s-wave scattering length formula 



i — I a.'<Heri Pr, 

with fs as the threshold amplitude, the best description of the observed energy dependence of the SAT- 
URNE total cross section data is achieved using a scattering length of asj^^gr? = ±(3.6 ±0.5) fm + i(2.0± 
0.9) fm, as can be seen in fig. 42 (solid line) with the absolute scale of the curve being scaled to fit the 
data. Such a large scattering length might be associated with the existence of a quasi-bound ^He rj state. 

In the right part of fig. 42 the extracted production amplitude squared |/(p*,)P is presented as 
function of the centre-of-mass momentum of the ejectiles. It should be noted that in the close vicinity 
of the production threshold and with a given cross section, the production amplitude squared |/(pJ^)P 
depends crucially on the momentum p* and, therefore, on the exact masses of the ejectiles. Here a mass 
of the ^He-midens of msj^e = 2.808392 GeV/c^ [396] and a meson mass of m^ = 0.54730 GeV/c^ [134] 
have been applied. While a constant amplitude would imply a behaviour according to pure s-wave 
phase-space considerations, the observed shape exposes a strong momentum dependence and can be 
reproduced well by the discussed formalism up to excess energies of Q ~ 7MeV, corresponding to 
p; fa 0.35 fm~\ 

However, recent measurements at higher excess energies from the GEM collaboration at COSY [392] 
and from the WASA/PROMICE collaboration at CELSIUS [393] clearly deviate from the predicted 
energy dependence of the cross section (fig. 42). Moreover, different to the SATURNE results, the 
centre-of-mass angular distribution of the emitted rj mesons strongly differs from an isotropic emission 
and exposes a peaking in forward direction. Motivated by results on the rj meson production in 7p 
and pp interactions, another approach to describe the production process has been evaluated [392] 
where the f] meson creation is described by the excitation of the Sn A^*(1535) resonance as shown in 
figure 43 (right). Neglecting effects of other resonances, the production amplitude can be expressed by 
a Breit-Wigner distribution [397]: 
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(51) 



with Er as the mass of the Sn resonance and as the corresponding width. Since the production 
threshold is located in the vicinity of the resonance mass Eji, the energy dependence of the width r(E) 
has to be considered: 

r(E) ^rR{b,^ + K^ + K,]. (52) 

Pr/R PnR 



The parameters bi indicate the branching ratios for different decay modes of the resonance while the 
centre-of-mass momenta of the rj and tt mesons arc represented by p*. A corresponding calculation, 
scaled to fit the GEM data point, is given in figure 42 by the dash-dotted line. Obviously, the prediction 
of this approach fails to describe the shape of the SATURNE data. This observation might indicate an 
influence of a strong ^He rj flnal state interaction only in the close vicinity of the production threshold. 
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Alternatively, this effect might be caused by a transition between different production processes. 

Concluding, the pd ^Herj reaction appears to be an interesting and important channel to study 
both the meson-nucleus final state interaction as well as reaction processes involving three nucleons. 
However, the body of data is far from being optimal, especially in the near-threshold region. New data 
from PROMICE/WASA and COSY-11 is currently under evaluation and will improve the experimental 
situation. 



6.2 Quasi— bound states 

Prom the theoretical point of view the existence of mesic nuclei or quasi-bound meson-nucleus systems 
is not excluded, however, up to now a compelling experimental proof for the formation of such a 
state is still missing. It is obvious that the possibility to create such states crucially depends on 
both the properties of the meson and the nucleus as well as on the sign and the strength of the 
meson-nucleus interaction. Since the vr-nucleon and i^T"*" -nucleon interaction have been found to be 
repulsive, it is unlikely that they will form quasi-bound states. Contrary, there are evidences for 
an attractive X^-nucleon interaction, however, corresponding experiments would suffer from the low 
production cross sections for K~ mesons. Furthermore, the Coulomb interaction is expected to screen 
corresponding physical observables and might lead to the formation of mesonic atoms bound by the 
Coulomb interactions. Different to this, the r] meson is uncharged and the observation of an attractive 
77-nucleon interaction led to speculations concerning the existence of 77-nuclear quasi-bound states. In 
the absence of 77-meson beams such states bound by the strong interaction would offer a new possibility 
to study the 77-nucleon interaction since the meson would be trapped for a relatively large time in 
nuclear medium. Predicted by Haider and Liu [387,388], detailed studies on the formation of such 
states have been the topic of recent investigations [138,140,389-391,398], resulting in different lower 
hmits of the atomic number A for which bound states should exist. While e.g. in [387, 388] the lower 
mass number is estimated in the framework of an optical model approach to be A = 12, there are 
considerations that even for nuclei of lower masses (d, t, ^He, ^He) bound systems might exist [138, 140], 
see also section 3. 

Experimental evidences for bound states have been found in the near-threshold production of 'q 
mesons in the reaction channel pd ^Herj. The unexpected large production amplitude as well as its 
rapid decrease with increasing energy (figure 42) is attributed to a strong s-wave final state interaction 
associated with a large r]~^He scattering length [140]. 

Such a large scattering length (see fit in figure 42) might imply the existence of a quasi-bound state. 
Assuming the ^Herj to form a quasi-bound state, it is therefore interesting to compare pd — > ^Herj 
production data with the ones of the reaction dd '^Herj, where according to the higher mass number, 
the formation of a bound state should be even more probable. A comparison of corresponding near- 
threshold data is presented in figure 44. 

A fit to the data using the scattering length formula (50) provides information on the scale of both 
the real and imaginary part of the scattering length of the meson-nucleus system but is unsensitive to 
the sign of the real part of auer]- However, information on this sign can be obtained using a lowest-order 
optical potential [140]: 



with the 77-nucleon reduced mass m.^, the mass number of the nucleus A and the 77-nucleon scattering 
length a^M ~ (0.52 + i0.25)fm. Using Gaussian nuclear densities with rms radii of p(r) = 1.63 fm 




(53) 
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(1.9 fm) for the rj^He (rj^He) system, meson-nucleus scattering lengths of 



ai^Heri) ~ (-2.3 + i 3.2) fm 
a{*Her}) ~ (-2.2 + i 1.1) fm 

with negative real parts for both systems have been obtained. 



(54) 
(55) 
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Figure 44: Averaged squared pro- 
duction amplitudes of the reactions 
pd ^Herj and dd '^Herj as 
function of the ccntrc-of-mass rj 
momentum. The curves are dis- 
cussed in the text. The data are 
taken from [386,399,400]. 



Fits to existing close-to-threshold data for both reaction channels using these scattering lengths 
and eq. (50) are presented in figure 44. Obviously, the data are reproduced well by this formalism. 
These results imply the existence of quasi-bound meson-nucleus states in the systems 7] ^He and 7] ^iJe, 
however, unambiguous evidence is still missing. Therefore, detailed investigations on both reaction 
channels as well as on reactions using heavier targets are highly recommended. 



6.3 Test of invariances 

Conservation laws restrict the quantum numbers for the production of mesons in the nucleon-nucleon 
interaction. Here especially the pion production is of considerable interest since it exhausts all inelas- 
ticity up to about 2GeV/c beam momentum in the nucleon-nucleon collision. 

As outlined in references [4,5,401,402] the conservation of isospin allows the total cross section for pion 
production to be expressed as an incoherent sum of three independent partial cross sections given by: 
(Jii, (Jio and (Toi, where the first (second) index indicates the isospin of the initial state (final state) 
two-nucleon system and (Tqq is ruled out due to isospin conservation arguments. 
The three cross sections can be determined via the reactions: 

pure (Til channels: 
pure (Jio channels: 
mixed Uiq -\- (Th channels: 
mixed dio + ctqi channels: 
mixed cth -|- ctqi channels: 

Charge symmetry relates those reactions denoted by a * to their mirror systems given as the first 
one in the same row. 



a['pp — ^ ppTT ), a{nn — > nnn )* 

a{pp ^ dn^), a{np—*d'K^), cr(nn — > dTr"")* 

a{pp — > pn7r+), a{nn npir')* 
a{np — > pmr^), 

a(np — > ppir"), a(np — > nmr'^). 
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The two reactions pd — > ^if 7r+ and pd — > ^Hen^, which permit the study of isospin breaking effects 
have been studied systematically [403] . Figure 45 compares the differential cross sections for these two 
reactions at two proton momenta close to the production threshold. At these momenta no significant 
differences in the differential cross sections were observed. In fact, interpolating the cross sections 
to the same four-momentum-transfer yields an expected average ratio R = a{pd — > TT~^)/a{pd — > 
^Hen^) 2.0 within the experimental uncertainties. However, the COSY-GEM collaboration seems to 
have indications that a clear effect of isospin breaking is observed at a momentum of about 1.57GeV/c, 
close to the 77-production threshold [404] . 




Figure 45: Differential 
cross sections for the two 
isospin symmetric reac- 
tions pd — > ^Hen'^ and 
pd — > ^Htt~^ at beam 
momenta of 800MeV/c 
and 750MeV/c, respec- 
tively. Data are from the 
GEM collaboration [403], 
the error bars are smaller 
than the symbols. 



cos(0) 



As pointed out in reference [405,406], in context of QCD isospin and charge symmetry are broken 
due to both i) the mass difference between up and down quarks and ii) the electromagnetic interaction. 
The observation of mixing between mesons belonging to the same SU{3) multiplets offers a unique 
possibility to study charge symmetry breaking effects. As suggested by Wilkin [407] a possible isospin 
symmetry breaking can manifest itself e.g. via external 77*^-77 mixing. Searching for such effects the 
authors of reference [405] are investigating the cross section ratio for the pd ^He 7r° and pd — > tt"^ 
reactions in particular for large relative pion-proton angles and at beam energies of the corresponding 
1] production threshold. The idea is based on the assumption that — besides the direct reaction — 
an intermediately formed fj meson undergoes a transition to the final Tr'^-meson for the pd — > ^He tt" 
reaction (pd — > ^Hefj — > ^Hen^) with a yield governed by the strength of the mixing angle while 
such an intermediate state can not occur for the pd — > 7r+ case. The authors [405] predict angular 
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dependent cross section ratios as: Rej,^=i8o° = 2.4 and -R6>p,r=o° = 2.03. The predicted symmetry breaking 
magnitude appears to be in reach of the experiments planned by the GEM collaboration at COSY. 

Another direct evidence for the isospin invariance breaking would be the observation of the non-zero 
yield for the isospin forbidden dd ^He 7r° reaction. The upgrade of the ANKE detection system at 
COSY will allow to search for the effect [408]. 
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7 Experimental facilities 



Studies of meson production at threshold require crucially high quality and precise quantitative knowl- 
edge of the accelerator beam quality especially due to the strong cross section dependence on the excess 
energy. The rapid growth of the meson production cross section in the vicinity of the kinematical 
threshold is evidently seen in figure 2. It has been shown in the previous sections that new precise 
data are available for different reaction channels. Here we present a very brief overview of technical 
features of such experimental facilities as: accelerators, targets and detector arrangements where only 
a representative selection can be given. 



7.1 Accelerators for hadronic physics at medium energies 

Initiated by the rather high costs for experimenting with antiprotons and aiming in a most efficient 
usage, storage rings have been developed. The possibility of controlling the emittance of such beams 
by electron as well as stochastic cooling, combined with some unique benefits of using internal targets 
lead to the construction of medium energy cooler rings as the low energy antiproton ring (LEAR) which 
might be regarded as the father for many facilities built. 



accelerators for hadron physics at medium energies 
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Figure 46: Accelerators 
medium energy physics. 
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The maximum available beam momentum limits the producible mass. In figure 46 a comparison is 
presented between proton and electron accelerators in the range of producing in the proton-proton and 
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7-proton collision mesons up to the 0-meson at 1.02 GeV/c^. This self-explaining figure demonstrates 
e.g. that the Uppsala accelerator CELSIUS can reach beam momenta to produce the 77-meson in the 
proton-proton scattering which is about the same mass production range as that of MAMI via the 

7-proton interaction. 

Complementary to the investigations with charged hadronic beams the high precision secondary 
neutron beam of the TRIUMF facihty [46] permits the close-to-threshold meson production in neutron- 
proton collisions. 

In Table 7.1 some accelerators used for hadron physics at medium energies are listed, which are or 
were important for this field. 



Table 7.1 Hadron accelerators for hadron physics at medium energy, 
LEAR and SATURNE stopped operation in 1997 and 1999, respectively. 



accelerator 


typ 


momentum 
MeV/c 


energy 
MeV 


beam 


cooling 


TRIUMF 


cyclotron 


1090 


500 


n,p 




lUCF 


synchrotron 


1090 


500 


A = 1, 2, 3 


e" 


SIN-PSI 


cyclotron 


1220 


600 


P 




LAMPF 


linear accelerator 


1460 


800 


P 




LEAR 


synchrotron 


2000 


1270 


P 


e' 


CELSIUS 


synchrotron 


2100 


1360 


A=l,2,... 


e~ 


COSY 


synchrotron 


3500 


2620 


A = 1,2 


e~ plus stochasitc 


SATURNE 


synchrotron 


3770 


3000 


A - 1...40 





7.2 Targets 

The choice of a certain target material for accelerator experiments as well as its chemical and mechan- 
ical properties is strongly correlated with the given experimental conditions and, certainly, with the 
reactions which should be investigated. Since in this article we concentrate on unpolarized meson pro- 
duction experiments off nucleons and deuterons, here we will restrict ourselves on corresponding targets. 

Proton and deuteron targets are commonly provided using hydrogen and deuterium as raw material. 
In the absence of free neutron targets deuterons can be used as substitute, considering the weakly bound 
deuterons as effective neutron targets and restricting on reactions with the bound protons acting only 
as spectator particles. For the sake of completeness we want to point out that in this case the Fermi 
motion of the nucleons has to be considered, especially in the regime of low excess energies. 

Near threshold production experiments are commonly connected with low cross sections in the 
order of picobarns or nanobarns. To obtain sufficient counting rates during experimental runs, high 
luminosities in the order of ~ 10'^*^ cm~^s~^ are necessary and can be achieved using targets of high areal 
density and intense accelerator beams. Neglecting technical limitations hke space charge effects, from 
the experimental point of view a combination of high intense ion beams and thin targets are preferable 
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since in this case the quality of the cooled accelerator beam is nearly uninfluenced during the passage 
of the target, granting well defined kinematical conditions as well as single particle interactions. In 
practice one has to distinguish between internal and external beam experiments. In the latter case a 
well prepared and extracted accelerator beam is shot onto a target and beam particles not interacting 
with the target are lost. Therefore, liquid or solid state targets of sufficient high areal density have 
commonly to be used in order to obtain the required high luminosities. Contrary, in internal beam 
experiments carried out at storage rings the circulating beam particles pass the target material 10^ to 
10^ times per second and projectiles not interacting with the target remain in the beam. Therefore, 
to achieve similar luminosities compared to external beam experiments only thin targets arc necessary. 
Moreover, to preserve the quality of the circulating beam and to obtain life times of the accelerated 
beam up to hours, effective areal target densities below ~ 10^^ atoms/cm^ are desirable and can be 
achieved using gas and cluster jets, pellets and thin fibers. Considering hydrogen or deuterium as raw 
material such targets can be realized in different ways. 

7.2.1 Fiber Targets 

Fiber targets with diameters of only a few micrometers are of special interest for storage ring experiments 
with the request of high target densities in combination with a point-like target and UHV conditions. 
Due to the relatively simple target installation, Att detection systems are possible. Using carbon for 
example as target material, fiber diameters down to 5//m are accessible, corresponding to ~ 2 • 10^^ 
atoms/cm^ [409] and allowing to compensate beam heating effects by beam coohng devices. Fiber 
targets for interactions off protons and deuterons naturally consist of chemical compounds like propylene 
such as the target for the internal beam experiment EDDA at COSY, where propylene fibers with 
diameters of ~ 5/xm are used. However, apart from those advantages one has to keep in mind that the 
mass of such CH2 fibers is mostly given by "carrier" material carbon, which infiuences the performance 
of the data aquistion system and, moreover, causes a physical background in the data which has to be 
considered. 

7.2.2 Liquid Targets 

Liquid hydrogen and deuterium targets with their high volume densities {plh2 ~ 0.07 g/cm^) are 
reasonable for the use at external beam experiments. Depending on the experimental requirements 
large target volumes even up to several liters are accessible. However, precision experiments on the 
near threshold meson production force to use only small target volumes in the interaction region in 
order to minimize energy losses and effects of multiple scattering in the target material. An example of 
an experimental facility using a thin liquid hydrogen target is the COSY-TOF installation [410,411]. 
Here hydrogen is liquefyed using a commercial cryogenic cold head, providing a compact target volume 
of several millimeters in diameter. Since such targets naturally have to be encased in special vessels, 
the choice of appropriate entrance windows is of great importance to minimize the physical background. 
At COSY-TOF for example thin windows of only ~ 1 //m Mylar foil with a diameter of ~ 4 mm are 
routinely in use, granting good signal to background conditions. 

7.2.3 Pellet Targets 

Pellet target installations can be used to prepare streams of frozen gas micro-spheres as targets for 
internal beam experiments. Similar to fiber targets, pellet targets can be associated to the class of the 
solid state targets and allow to provide high areal densities in combination with point-like interaction 
regions, well suited for An detector geometries. Moreover, contrary to fiber targets these installations 
allow to produce pure hydrogen or deuterium pellets without any carrier material, however, at the cost 
of an increased experimental effort and vacuum conditions of ~ 10~^ mbar in the scattering chamber. 
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At the CELSIUS ring in Uppsala such a pellet target was build up for the WASA experiment [413-415]. 
The pellets are produced by pressing liquid hydrogen at triple point conditions (14 K, 72 mbar) through 
a glass nozzle of 20 /im diameter. By means of a piezo-electric transductor connected to the nozzle the 
stream of hydrogen is broken into uniformly sized micro-spheres (pellets) as shown in figure 47. This 
device was designed to provide a stream of hydrogen pellets with diameters of 50 /xm at a repetition 
rate of ~ 68 kHz. While the pellets itself are locally very thick (~ 10^° atoms/cm^), the effective areal 
density in the scattering chamber is in the order of 10^^ atoms/cm^ due to the minimized geometrical 
beam-target overlap. 



t I 




Figure 47: Formation of deuterium pellets with diameters of 

~ 50 /um, the photograph was taken at the pellet target of 
the CELSIUS/WASA detector at CELSIUS in Uppsala. 




7.2.4 Gas Jet Targets 

Supersonic gas jet beams can be used as targets for both internal and external beam experiments. 
Produced by feeding Laval-nozzles with gases like hydrogen or deuterium, high areal target densities 
(~ 10^^ atoms/cm^) can easily be achieved close to the nozzle exit [416]. Similar to pellet targets, proton 
or deuteron targets without effects of windows or carrier materials can be provided by operating the 
installation with gases of high purity. However, since the volume density of a gas jet beam drastically 
decreases both in longitudinal and transversal direction, the nozzle has to be placed very close to the 
interaction point in the scattering chamber. Therefore, the advantage of high target densities of such 
installations is commonly connected with higher gas loads in the scattering chamber and a limited 
spatial resolution of the interaction point if no vertex detector is used. The great advantage of such 
targets is the possibility to use almost all gases as target material and to adjust the areal density by 
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orders of magnitude by simply changing the gas input pressure. 
7.2.5 Cluster Targets 

Similar to conventional gas jet targets, cluster beams [417] can be produced by using fine nozzles with 
a convergent-divergent shape (Laval-nozzle) and, therefore, almost all gases can be used. Operated at 
appropriate low temperatures, a partially condensation of the used gas is possible and clusters with 
sizes of 10^ to 10^ can grow. Separated from the surrounding conventional gas beam by a set of 
collimators (see fig. 48), cluster beams with a homogeneous volume density, a sharp boundary, a small 
angular divergence and well defined beam dimensions can be prepared for interaction with circulating 
accelerator beams [412]. 



Cluster beam installations are frequently used at cooler rings. Using Laval-nozzles with diameters 
of 11-16 yum, low nozzle temperatures of ~25 K and gas input pressures of ~18 bar, hydrogen and 
deuterium cluster beams with areal densities in the order of ~ 10^^ atoms/cm^ in combination with 
UHV conditions are accessible in the interaction region. Furthermore, similar to gas beam targets the 
areal density can easily be adjusted over orders of magnitude by changing the gas input pressure or the 
nozzle temperature [412]. 

7.3 Typical detector arrangements 

Special efforts are made for an optimal detector design to match the required conditions for kinemat- 
ically completeness and undisturbed determination of the four-momentum vector of the ejectiles. A 
kinematically complete event is characterized by knowing all four-momentum vectors of the participat- 
ing particles. The experimental information is limited by: i) statistical fluctuations, ii) the precision 
for the determination of the components of the four-momentum vectors, iii) contributions from any 
kind of background produced by means of the apparatus or due to physical reasons, iv) the efficiency of 
the detector covering not 100 % of the phase-space due to size or construction limitations and finally 
v) the reconstruction efficiency due to unavoidable inadequate performance whatsoever. Here all col- 
laborations have to meet the optimum between wishful thinking, technical possibilities and affordable 
costs. 




saturated gas 
[Z3 cluster-beam 



Figure 48: Schematic drawing of the prin- 
ciple for the formation of clusters. 
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7.3.1 Example of a 47r internal-detector 

The high luminosity and iir detector CELSIUS/WASA (Wide Angle Shower Apparatus) facility [418] 
is primarily 




Figure 49: Top: CAD-view of the WASA detector facility 

Bottom: Cross section of the central and forward part of the WASA detector facihty. 
A version of the article including a photo of better quality can be found at the COSY-11 homepage: 
http : / /ikpellOl.ikp.kfa — juelich.de/cosy — 11 / pub / List_of -Publications. htmlij^papers 

designed for both measuring rare decay processes of light mesons and investigating meson produc- 
tion studies in the threshold region. The target and beam particles tend to be scattered at forward 
angles while the products of the meson decay will be distributed rather isotropically. The detector is 
constructed following these constrains and in full operation an annual production rate of about 10 77's 
is expected in the pp interaction. 

The essential parts of the CELSIUS/WASA facility are: 

• The pellet target with: 

Pellets diameter and frequency of 25 - 50 jivci and a few kHz, respectively, 
where the beam diameter is about 3 mm, 5 x 10 atoms/cm^. 
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• The forward detector with: 

A segmented forward window counter (FWC) for fast timing in the first trigger level. 

A forward proportional chamber (FPC) out of straw chambers to measure the scattering angle. 

A forward trigger hodoscope (FTH) built out of three layers of circular scintillators used for: 

fast trigger decisions, rough track position determination and particle identification. 
A forward range hodoscope (FRH) as a plastic scintillator calorimeter for measurements of 

particle energy and identification. 
A forward range intermediate hodoscope (FRI) adding a position sensitivity to the FRH. 

• The zero degree spectrometer to filter out low rigidity reaction products employing HPGe- 
or telescopes of Si-/x-strip detectors, scintillation counters and CsI(Tl) crystals. 

• The central detector with: 

A scintillator electromagnetic calorimeter (SEC) consisting out of 1012 CsI(Na) 

between the superconducting solenoid (SCS) and the iron yoke. 
A mini drift chamber (MDC) with 1738 straw tubes for measuring tracks of charged particles. 

The CELSIUS/WASA detector is completed and regular data taking has started in late 2001 on some 
TT-production reactions where the vr^'s are detected by their 77 decays. The tagging of the ?7-production 
and the detection of the jy's by their 77 and Sir^ decays has been proven to work. In other words there 
is a high luminosity hadron facility ready for exciting physics. 



7.3.2 Example of a centre— of— mass 47r external— detector 

The COSY-TOF detection facility at the external beam line of COSY is a multi purpose modular 
arrangement of segmented barrel and circular planar scintillation detectors for measuring the time-of- 
flight and the energy loss of charged reaction products. 

In Fig. 50 two CAD-drawings of a short and an extended version as well as a photograph of the 
outer barrel in it's longest assembly is shown. The modularity serves several purposes but essentially 
helps to adjust the system to the optimal dimensions according to the reaction under investigation. 

Three scintillators placed 1 m, 0.5 m and 0.1 m in front of the target are used to veto a beam halo. 
Further the following components are employed: 

• A few mm^ large liquid hydrogen/deuterium target with very thin windows [410,411] 

• Different granulated start counter systems optimized for the reaction type under investigation. 

• Up to three barrels equiped with granulated cylindrical detectors for timing, spatial resolution 
and energy loss determinations. 

• Granulated three layer stop counters, divided into inner- and outer circular scintillation detector 
hodoscopes [419]. 

• A scintillation calorimeter for energy determination of charged and neutral hadrons, placed behind 
the stop detector, see the short version of the facility, left in fig. 50. 

• In preparation a straw drift chamber hodoscope, to be installed behind the start detector systems. 
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Figure 50: Schematic view of the COSY-TOF detector facihty in two different set-ups, the short and 
the longer version. Bottom: A photo from the very long version of the COSY-TOF detector. 
A version of the article including a photo of better quality can be found at the COSY-11 homepage: 
http : / /ikpellOl.ikp.kfa —juelich.de/cosy — 11/pub/List-of -Publications. html^papers 

7.3.3 Example of a dedicated internal threshold detector 

The COSY-11 installation is an internal facility at COSY designed for threshold meson production 
studies. Close-to-threshold, due to the small transvers momenta, all ejectiles arc limited to a narrow 
forward cone in the laboratory system and therefore a high geometrical acceptance is attainable with a 
comparatively small detector set-up. The essential components of the COSY-11 installation are: 

• A cluster target located in front of a regular C-shaped COSY dipole producing H2 or D2-clusters 
of up to 10 ^ molecules and densities of about 10 atoms/cm^. 

• An exit window (187 x 7.6 cm^) of 30 /xm Al and 2 x 150 /xm carbon fibers to separate the 
ultra high beam line vacuum from atmospheric pressure. It provides high acceptance compared 
to conventional magnetic spectrometers especially in the horizontal plane. 

• Two sets of drift chambers stacks (Dl and D2 in fig. 51) with six and eight planes, respectively 
and a position resolution of about 200 /im in each plane, for particle track reconstruction and 
momentum analysis through the magnetic field for positively charged ejectiles. 
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Figure 51: Top: Top view of the COSY-11 zero degree facility. 

Bottom: A photographic view of the COSY-11 detector at the COSY ring. 
A version of the article including a photo of better quality can be found at the COSY-11 homepage: 
http : / /ikpellOl.ikp.kfa — juelich.de/cosy — 11/pub/List-of -Publications. html^papers 



• Two scintillator hodoscopes (SI and S2 in fig. 51) for fast level trigger determination. 

• A large area scintillation wall [420] (2.2 x 1 m^) S3 in fig. 51 placed 9.1 m downstream, for fast 
trigger, time-of-fiight and consecutively invariant mass determination. 

• A detection system for negatively charged mesons consisting out of a scintillator detector and an 
array of silicon pad detectors, both mounted inside the dipole gap. 

• A monitor system for elastically scattered particles, utilized for luminosity determination. 

• In preparation for an extension: two wire chambers with a hexagonal drift structures to be 
positioned at the side - and at the front exit windows and a threshold Cerenkov counter behind 
the SI hodoscope for pion background subtraction. 
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8 Outlook 



Storage rings with the flexibihty of different projectiles, variable beam momenta, extremely good phase- 
space cooling and with the unrivaled options of internal and external detection systems provide an 
excellent tool for investigations on meson production at threshold. Especially the high quality of 
relative and absolute energy resolutions available at modern accelerators are crucial for such studies. 

In the present report, discussing the close-to-threshold meson production in hadronic interactions, 
the natural restriction to investigations using only hadronic probes was made though interesting and 
informative results were certainly obtained with electromagnetic tools as well. 

The theory of strong interactions (QCD) has a convincing predictive power at high energies. Per- 
turbative expansions in the coupling constant no longer converge for medium energy scales and the 
structure of QCD becomes increasingly complicated. The threshold meson production in nucleon- 
nucleon scattering offers valuable aspects for studying production, decay and interaction of hadrons 
and investigating the reaction mechanisms, the dynamics of the system and appropriately chosen sub- 
systems as well as the structure of hadrons themselves. Such experiments probe the sensitivity of 
nucleon-nucleon interactions at distances of < 1 fm, while the long-range part of the interaction is 
suppressed. Meson production in pp — > ppX reactions (with X — 7r°,77, is a process of strong 
inelasticity in proton-proton scattering. 

8.1 Polarized beams and/or targets 

In the present paper discussions utilizing polarized beam and/or targets are omitted. At the upgraded 
lUCF cooler ring a polarized proton beam in combination with a windowless polarized hydrogen target 
now permits measurements of all spin correlation coefficients for — > NNtt reactions. Very recently 
initial results for spin correlations were reported [421]. The goal of such experiments is to qualify the 
energy dependence and the importance of higher partial waves by determining analyzing powers and 
spin correlations. Here the hadron physics community certainly faces the challenge of new experimental 
techniques with necessary improvements of theoretical developments, since polarization observables 
are sensitive indicators of reaction mechanisms, relative contribution of partial waves and transition 
amplitudes. 

The production of rj mesons near threshold has been measured by different groups [61-63,65]. Very 
close-to-thrcshold only S-wavcs are supposed to contribute to the production process and in fact the 
S-wave production is large due to the presence of the N*(1535)Sii resonance which seems to act as an 
intermediate state. 

Both TOP [98] and COSY-11 [66] took recently rather high statistics data on the rj production 
at Q ~ 15.5 MeV with fairly flat angular distributions as expected for an s-wave production. Still 
- whether some data suggest or don't suggest - differential cross section measurements might be too 
insensitive to reveal e.g. a small ^rf-contribution. Data using a polarized beam, however, might indicate 
effects arising from an interference term of the kind: 

2Re{AssA*sa} cos2e; (56) 

where Ass is the amplitude for the production of a final pp pair with angular momentum Lpp — and 
/ = while for Asd the r] meson has an angular momentum of = 2. 

The importance of d-waves in r] production is well known from both ir'p rjn [422] and 7p rjp [397] 
reactions and the fact that this interference term is negative suggests that the r^-production in proton- 
proton scattering is governed mainly by p-exchange [99]. It is the aim of the COSY-11 group, followed 
by a recent feasibility study [181] to investigate the presence of this and other higher partial waves 
through the study of the analyzing power Ay. 
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Though the asymmetry should also be studied with respect to the final proton directions, the s-d 
interference contributes significantly to the rj analyzing power defined as: 



Tr(MMt) ^^'^ 

with the transition matrix element Ai and the Pauli matrix ay. 

To lowest order in the rj d-wave amplitude, the analyzing power is proportional to the imaginary 
part of the s~d interference [423] , 

A. « 2 (58) 

where (p* = corresponds to the plane of the beam (perpendicular to the polarization direction). On 
angular momentum grounds, this signal is expected to vary like ip^y, i.e. roughly like Q. 

The Pii(1440) resonance could also be of importance for the r] production [424], leading to p- 
waves in the final state. The interference between the Ps and Pp amplitudes gives an analyzing power 
proportional to: 



\^Ss\ 

and such a signal has been observed in the pp ppix^ reaction at low energies [38] . The presence of the 
three p-waves in the interference might suppress such a term close-to-threshold for 77 production and 
no unambiguous signal for them has been seen in the unpolarized differential cross section. Nevertheless 
one must be prepared for their existence. The Ss-Sd interference gives a maximal Aj, signal for ^* ?a 45° 
and vanishes at 90°, whereas the Ps-Pp signal should be largest at this point. 

Both the differential cross section and the analyzing power should be studied in order to detect ef- 
fects associated with the onset of higher partial waves. It should be noted that the effects of incomplete 
acceptance are much less of a problem for Ay than for the da/dft measurement provided that the beam 
characteristics of the two polarization states arc similar. Finally, the results should help to settle the 
on-going discussion of whether the rj production is dominated by p [175,176,425,426], lu [178] or r] [19] 
exchange. 

First data on the rj meson production with a polarized proton beam in the reaction have been taken 
at the internal experiment facility COSY-11 [181]. The measurement was performed with a proton 

— * 

beam momcutum Pp = 2.096 GeV/c corresponding to an excess energy of Q = 40 MeV. 
In Figure 52 a comparison of the data (solid triangles) with theoretical predictions for Q = 37 MeV from 
Faldt and Wilkin [177] (dotted line) and Nakayama [20] (solid line) is shown. In addition, the theoretical 
calculations for Q = 10 MeV are plotted, both based on the one meson exchange model. While Faldt 
and Wilkin predict a dominance of p exchange (dashed line) , Nakayama concludes a dominant tt and rj 
exchange (sohd line). 



89 



Figure 52: Theoretical predictions for the 
analyzing power according to [177] (dotted 
line) and [20] (solid line) for Q = 10 and 
37MeV. The solid triangles represent the 
measured values [181] at Q = 40 MeV. 



30 60 90 120 150 180 

Though the data are consistent with Ay = within the error bars given, there seems to be a slight 
increase of Ay towards backward scattering angles. First results of the analyzing power - deduced 
from the very first feasibility study run - suffer from statistics but tend to differ from both theoretical 
calculations, although the predictions of [177] show a quite good match. For a more precise statement 
further measurements are needed which could benefit from higher polarization now available at COSY. 

8.2 Final state interaction 

An important physics aspect of the modern cooler ring facilities is the interaction of particles in the final 
state. A detailed study of the Dalitz plot of, e.g. the reaction pp — >■ pK~^A measured at COSY-11 [328], 
made possible the extraction of valuable information about the A-N interaction. 

The pn — > drj reaction, measured at TSL in Uppsala [47-50], hmited significantly the range of values 
for the r)-N scattering length. This quantity is the relevant input for calculations studying the existence 
of 77-nucleus bound states. 

8.3 Baryon resonances 

Throughout the article the importance of an intermediate excitation of baryon resonances has been 
pointed out. 

At excess energies of 1-2 GeV, baryon resonances begin to play an important role in nucleon-nuclcon 
collisions. Most of our understanding about these resonances comes from pion and photon induced 
reactions on the nucleon. Current quark models predict many more resonances than can be experi- 
mentally extracted from the data. A large program at both JLAB and ELSA focuses on this so-called 
missing resonance problem and the search for yet unknown baryonic states in electro- nuclear reactions. 
However, besides the possibility that resonances thus far simply went unnoticed because the quality 
of the data was insufficient, another explanation might be that these resonances couple only weakly 
to photons or pions. An accelerator like COSY, built for nucleon-nucleon as well as nucleon-nucleus 
collisions will provide for the necessary experiments to fill that gap in the future. 

8.4 Symmetries 

The investigation of symmetries and their breaking has always led to significant insights into the inner 
workings of nature. Already a large program at COSY [405,406] is studying charge-symmetry breaking 



Q=10MeV 



0.25 



Q=37MeV 
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- an effect closely linked to the difference in the current masses of the light quarks. 

Another interesting aspect of charge-symmetry breaking is its high sensitivity: Once, for instance, the 

/o-flo mixing amplitude is isolated, it will provide us with insights about the nature of scalar resonances. 

Studies about the structure of these resonances are under investigations at COSY [278-280]. 

At the CELSIUS ring in Uppsala the CELSIUS/WASA detection system is in operation now and exciting 

experiments especially on r;-meson decay studies will start. 

8.5 Meson production in the three— nucleon system 

Naively, one would think that near-threshold meson production in pd collisions can be fully understood 
in terms of a short-ranged A^A^ NNn process, with the extra nucleon acting as a spectator. However, 
backward tt production in the pd ^Hen^ cannot be explained in terms of a spectator model. In the 
case of meson production in pd collisions the number of participating amplitudes is small as well, even 
though the transferred momentum is large. Hereby a unique environment is available for studying the 
interaction of more than just two nucleons, e.g. three-nucleon forces. It is likely that eventually the 
research on three-nucleon forces will focus on this process. Again, the role of polarization observables 
is important. 

8.6 Apparative aspects 

For the experiments, the development of suitable experimental configurations, comprising polarized 
targets and associated detector systems that facilitate measurements of polarization observables is 

necessary. 

Present and future experiments become possible only through the outstanding performance of the 
cooler accelerators. With the new injector for COSY the situation will be further improved. It will 
be possible to fill the machine up to the space-charge limit with polarized protons and deuterons, 
whereby electron-cooled beams at injection energy and stochastically-cooled beams at high energy 
of unprecedented brilliance will be available to the users. Efficient correction schemes to preserve 
the polarization of protons during acceleration have already been established, while for deuterons no 
difficulties are anticipated, because no depolarizing resonances occur in the energy range of COSY. New 
polarization experiments in the future will certainly call for further refinement of the spin-manipulation 
systems in COSY, e.g. longitudinally polarized internal and external beams, for which suitable spin 
rotators have to be developed. 

8.7 Epilogue 

The above mentioned examples - which will be the subject of a Summer School and Workshop on COSY 
Physics [427] - demonstrate the richness of physics in meson production reactions. The drawback of 
this wealth comes at the expense that, apart from rare cases, it is difficult to extract a particular piece 
of information from the data, e.g. both resonances and final state interactions modify the invariant 
mass plot in a coherent superposition. 

Fortunately, the use of polarization in the initial state acts as a spin filter and different contributions 
to the interaction can be singled out, because they show up in the angular distributions of different 
spin combinations. It is important to stress, that as long as only the lowest partial waves contribute 
significantly to the process under study, the number of independent observables from experiments with 
polarized beam and target is sufficient to allow for a model independent partial wave analysis. 
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